Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,507,424 G→C pseudogene (228/255 nt) insO → IS911B regulator fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,507,4240GC97.1% 218.5 / ‑0.2 68pseudogene (228/255 nt)insOIS911B regulator fragment
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (30/36);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

AACACAGGTGCTCCAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTGTTCTCATTTGTTCCCCGCTGCCAAGGACTCTGAGGA  >  NC_000913/4507182‑4507670
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
aacaCAGGTGCCCCAATGAAAAAAAGAAATTTCAGCTCAGAGTTTACACGCGGATCCGCTCAACTGGTTGTCGACCAGAAAAACTCGGTGGCAGATGCCTCGAAAGCTATGGATGTTGGTCTTTCCACAATAACAAGATGGGTCAAACAACTGCGTTATGAGCGTCAGGGCAAAACAACAAAAGCCACTCCGATAACACCTGAAGAAATCGAAATACGTAAGCTGAGGAAAAAGCTATAACGCATTGaaac                                                                                                                                                                                                                                                <  2:99045/251‑2 (MQ=255)
  cacaggatgtCCAAAGAAAAAAAGAAATGTCAGCGTAGAGTTTAAACGCGAATCCGCTCAACGGGTTGGTGACCAGAAAGACTCGGTGGCAGATGCCGCGAAAGCTATGGATGGTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATgg                                                                                                                                                                                                                                              <  2:522005/242‑1 (MQ=255)
        ggCTCCAAGGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGTTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTATCCACAATGACAAGATGGGTCAAACAGCTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATg                                                                                                                                                                                                                                        <  2:326371/250‑1 (MQ=255)
         gCTCCAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGa                                                                                                                                                                                                                                       <  2:450831/251‑1 (MQ=255)
          cTCCAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGga                                                                                                                                                                                                                                             <  1:120343/244‑1 (MQ=255)
          cTCCAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGga                                                                                                                                                                                                                                             >  2:120343/1‑244 (MQ=255)
             cAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAat                                                                                                                                                                                                                                    >  1:183173/1‑250 (MQ=255)
             cAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCACGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAata                                                                                                                                                                                                                                   <  1:217260/251‑1 (MQ=255)
             cAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAAGCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACTCCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAata                                                                                                                                                                                                                                   <  2:55744/251‑1 (MQ=255)
                 gAAAAAAAGAAATTTCAGCGCAGAGGTTAAACGCGAACTCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGAAGCCTCGAAAGCTATGGATGTTGGCCTTTCCAAAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAatatga                                                                                                                                                                                                                                <  2:636286/250‑3 (MQ=255)
                  aaaaaaaGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGAGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTaa                                                                                                                                                                                                                               >  2:171214/1‑250 (MQ=255)
                   aaaaaaGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTaaaa                                                                                                                                                                                                                             <  2:596615/251‑1 (MQ=255)
                    aaaaaGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTaaaa                                                                                                                                                                                                                             >  1:340009/1‑250 (MQ=255)
                         gAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGACAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGCACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACCACGCATTGAAATGGAGAATGAAATATTAAAAAGGCt                                                                                                                                                                                                                        >  1:639306/1‑250 (MQ=255)
                          aaaTTTCAGCGCAGAGTTTAAACTCGAATCCGCCTAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTAc                                                                                                                                                                                                                      <  2:540077/251‑1 (MQ=255)
                          aaaTTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTAc                                                                                                                                                                                                                      <  1:509040/251‑1 (MQ=255)
                                 aGCGCAGAGTTTAAACGCGAATCCGCTAAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctc                                                                                                                                                                                                               <  1:472682‑M2/251‑10 (MQ=255)
                                           ttAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtca                                                                                                                                                                                                      >  2:206881‑M2/1‑232 (MQ=255)
                                             aaaCGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCAGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagac                                                                                                                                                                                                   <  1:491102‑M2/251‑22 (MQ=255)
                                              aaCGCGAATCCGCGCAAATGGTTGTTTAAAAGAAATACACGGTGGCAGATGCTGAGAATGCTGGGGATGTTGGCCTTTCCTCAATGACAAGATGGTTCAAACAACTGCTTGATGAGCGTCAGGGCAAAACACAAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATGGAAATGGAGAATGAAATATTAAAAAGTCTtccgcgctcttgatgtcagact                                                                                                                                                                                                  >  2:379160‑M2/1‑229 (MQ=255)
                                                   gAATCCGCTAAACTGGTTGTTGACCAGAAATACACGGAGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAAGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGATGAGGAAAAAGCTAAAACGCATTGAAATGGAGAATGAAATATTAAAAAGGATaccgcgctctagaagtaagactccctg                                                                                                                                                                                             <  2:338777‑M2/251‑28 (MQ=255)
                                                          cTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagtt                                                                                                                                                                                      <  1:118‑M2/251‑35 (MQ=255)
                                                            cAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTGCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttc                                                                                                                                                                                     >  1:91114‑M2/1‑215 (MQ=255)
                                                             aaCTTGTTGTTGACCAGAAATACCCGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCGTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagacgcctggaacagttcta                                                                                                                                                                                   >  2:246475‑M2/1‑214 (MQ=255)
                                                                            aGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaacc                                                                                                                                                                    <  2:260615‑M2/251‑53 (MQ=255)
                                                                             gAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactc                                                                                                                                                                   <  1:518094‑M2/251‑54 (MQ=255)
                                                                                aTACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcaga                                                                                                                                                                >  1:171214‑M2/1‑195 (MQ=255)
                                                                                        tagcaGATGCGGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACAGCGTGATGAGCGACAGGGCAAAACACCAAAAGCCTCTCCGAAAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAACTACAACGCATTGAAATGGAGACTGAAATATTAAAAAGGCTaccgcgctcttgatggcagactccctgaccagttctcgataatcgggaaactcacagagcatt                                                                                                                                                         <  2:314525‑M2/248‑64 (MQ=255)
                                                                                              aTGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgt                                                                                                                                                  <  1:245‑M2/251‑71 (MQ=255)
                                                                                                       aaGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacact                                                                                                                                         <  2:192610‑M2/251‑80 (MQ=255)
                                                                                                          cTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctg                                                                                                                                      <  1:70426‑M2/251‑83 (MQ=255)
                                                                                                                      gCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcagtatcctgtggtcacactctgccatgtgttcgg                                                                                                                          >  1:314525‑M2/1‑157 (MQ=255)
                                                                                                                         tttCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggt                                                                                                                       <  2:670982‑M2/251‑98 (MQ=255)
                                                                                                                          ttCCACAATGACAAGATGCGTCAAACAACTGCTTGATGATCGTCAGGACAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTTCAACGCATTGAAATGGATAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagtctccctgtacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggt                                                                                                                       >  2:183173‑M2/1‑153 (MQ=255)
                                                                                                                           tCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggtt                                                                                                                      >  2:113640‑M2/1‑152 (MQ=255)
                                                                                                                             cacaATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATTAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcat                                                                                                                   >  1:598867‑M2/1‑150 (MQ=255)
                                                                                                                              acaATGACAAGATGGGTCAAACAACTGCGGGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatc                                                                                                                  >  2:157206‑M2/1‑149 (MQ=255)
                                                                                                                                    aCAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGCGGAAAAAGCTACAACGCATTGAAATGGCGAATGAAATATTAAAAAGGCTaccgcgctcttgctgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagc                                                                                                             <  1:85815‑M2/250‑108 (MQ=255)
                                                                                                                                    aCAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagc                                                                                                             <  1:493600‑M2/250‑108 (MQ=255)
                                                                                                                                     cAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagc                                                                                                             <  1:382304‑M2/249‑108 (MQ=255)
                                                                                                                                           ggggCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctaca                                                                                                      >  2:351215‑M2/1‑136 (MQ=255)
                                                                                                                                             gTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagat                                                                                                   >  1:35991‑M2/1‑134 (MQ=255)
                                                                                                                                                 aacaacTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactg                                                                                               >  1:587390‑M2/1‑130 (MQ=255)
                                                                                                                                                 aacaacTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAAACGAAATACGTAAGCGGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactg                                                                                               >  2:493600‑M2/1‑130 (MQ=255)
                                                                                                                                                  acaccTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgagggcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactgg                                                                                              >  2:350598‑M2/1‑129 (MQ=255)
                                                                                                                                                     acTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaa                                                                                            <  1:520742‑M2/250‑125 (MQ=255)
                                                                                                                                                        gCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaa                                                                                          >  1:670982‑M2/1‑123 (MQ=255)
                                                                                                                                                                     cAGGGCAAAACACCAAAAGCCTCTCCGATTACACCCGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGCAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgctgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacaatctgcattgtgttcggggatcatcgcagcagctacagttactggaaaaaccgtcatgaaaacc                                                                           <  2:98645‑M2/251‑142 (MQ=255)
                                                                                                                                                                          cAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagac                                                                      <  1:2382‑M2/251‑147 (MQ=255)
                                                                                                                                                                             aaCACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaacccagacggc                                                                   <  2:160079‑M2/251‑150 (MQ=255)
                                                                                                                                                                                 ccAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagac                                                               >  2:340009‑M2/1‑98 (MQ=255)
                                                                                                                                                                                  cAAAAGCCTCTCCTATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacg                                                              >  2:2382‑M2/1‑97 (MQ=255)
                                                                                                                                                                                   aaaaGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgg                                                             >  1:334266‑M2/1‑96 (MQ=255)
                                                                                                                                                                                      aGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggct                                                          >  1:208152‑M2/1‑93 (MQ=255)
                                                                                                                                                                                           ctcCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtatt                                                     <  1:411672‑M2/251‑164 (MQ=255)
                                                                                                                                                                                             ccGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtatt                                                     >  1:260615‑M2/1‑86 (MQ=255)
                                                                                                                                                                                                aTAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTACAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgc                                                 <  1:227835‑M2/250‑168 (MQ=255)
                                                                                                                                                                                                aTAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattccgc                                                 <  1:628708‑M2/250‑168 (MQ=255)
                                                                                                                                                                                                       cAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggt                                         >  1:272061‑M2/1‑76 (MQ=255)
                                                                                                                                                                                                                  gAAATAAATAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggacacactctgccatgtgttcggggttcatcaccgcagctacagatactggaaaacccatcctgaaaaaccagacggcagacgggctgtattacgcagtcagtgaattgaactc                               <  2:631024‑M2/250‑186 (MQ=255)
                                                                                                                                                                                                                    aaTACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggtacttgagctacat                            >  2:382304‑M2/1‑63 (MQ=255)
                                                                                                                                                                                                                           aaGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggtacttgagctacatggcatca                     >  1:160079‑M2/1‑56 (MQ=255)
                                                                                                                                                                                                                           aaGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggtacttgagctacatggcatc                      <  1:5366‑M2/250‑195 (MQ=255)
                                                                                                                                                                                                                                      aaaGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggtacttgagctacatggcatcagccacggttcg          <  1:198662‑M2/251‑207 (MQ=255)
                                                                                                                                                                                                                                        aGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTATTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTGTTCTCATTTGTTCCCCGCTGCCACGGActc        >  2:215236/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgacatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggaagacgggctgtattacgcagtccggtacttgagctacatggcctcagccccggttcggc        <  2:495355‑M2/251‑209 (MQ=255)
                                                                                                                                                                                                                                             cAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggtacttgagctacatggcatcagccacggttcggccggag   >  1:224230‑M2/1‑38 (MQ=255)
                                                                                                                                                                                                                                               aCGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggtacttgagctacatggcatcagccacggttcggccggagc  <  1:286888‑M2/250‑215 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
AACACAGGTGCTCCAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTGTTCTCATTTGTTCCCCGCTGCCAAGGACTCTGAGGA  >  NC_000913/4507182‑4507670

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: