Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,194,650 G→A 5.2% P101L (CCG→CTG)  rluE ← 23S rRNA pseudouridine(2457) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,194,6500GA 5.2% 632.3 / 15.6 191P101L (CCG→CTG) rluE23S rRNA pseudouridine(2457) synthase
Reads supporting (aligned to +/- strand):  ref base G (96/84);  new base A (7/3);  total (104/87)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 192 total aligned reads displayed.

AACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGC  >  NC_000913/1194503‑1194796
                                                                                                                                                    |                                                                                                                                                  
 aCCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGc                                                                                                                                                  >  4:935166/1‑149 (MQ=255)
  ccACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCt                                                                                                                                                 >  1:1018458/1‑149 (MQ=255)
    aCGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTgg                                                                                                                                               <  1:1564779/149‑1 (MQ=255)
     cgcgGGTTCGTCAACCAGCTCCGCGCCGGCGGGCGGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCTGCTggg                                                                                                                                              >  1:825136/1‑149 (MQ=255)
      gcgGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGt                                                                                                                                             <  2:93851/149‑1 (MQ=255)
           ttCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAAc                                                                                                                                        <  3:1356049/149‑1 (MQ=255)
              gTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAAcgcg                                                                                                                                     >  2:1108610/1‑149 (MQ=255)
               tCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTAACCCAGCCGGG‑TAAAcgcac                                                                                                                                    <  2:1027651/149‑3 (MQ=255)
               tCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACgcgc                                                                                                                                    >  1:663926/1‑149 (MQ=255)
                   ccAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGc                                                                                                                                <  4:794725/149‑1 (MQ=255)
                    cAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCa                                                                                                                               >  2:1464507/1‑149 (MQ=255)
                     aGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCaa                                                                                                                              >  3:961501/1‑149 (MQ=255)
                        tCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAAcgc                                                                                                                           >  4:6737/1‑149 (MQ=255)
                               cggcggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTa                                                                                                                    >  1:923970/1‑149 (MQ=255)
                                 gcggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTAtt                                                                                                                  <  4:1275573/149‑1 (MQ=255)
                                 gcggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTAtt                                                                                                                  <  4:829384/149‑1 (MQ=255)
                                  cggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTg                                                                                                                 <  1:1193060/149‑1 (MQ=255)
                                     gCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTc                                                                                                              <  3:151259/149‑1 (MQ=255)
                                      cAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCa                                                                                                             <  3:1284138/149‑1 (MQ=255)
                                               gCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCa                                                                                                    >  1:1370288/1‑149 (MQ=255)
                                                ccATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTc                                                                                                              >  4:277504/1‑138 (MQ=255)
                                                ccATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTc                                                                                                              >  2:179896/1‑138 (MQ=255)
                                                ccATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTc                                                                                                              <  3:277504/138‑1 (MQ=255)
                                                  atcatTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAAcc                                                                                                 >  2:585974/1‑149 (MQ=255)
                                                     atTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCtt                                                                                              >  4:1106563/1‑149 (MQ=255)
                                                       ttAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCg                                                                                            >  2:43015/1‑149 (MQ=255)
                                                           ggTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTAt                                                                                        <  1:585974/149‑1 (MQ=255)
                                                            gTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTc                                                                                             <  2:1618637/143‑1 (MQ=255)
                                                               aCGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCGGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCCATcgcg                                                                                    <  2:1217775/149‑1 (MQ=255)
                                                               aCGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATcgcg                                                                                    <  4:294804/149‑1 (MQ=255)
                                                                cGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCgg                                                                                   <  2:1056902/149‑1 (MQ=255)
                                                                 gCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAAcc                                                                                                 <  4:762080/134‑1 (MQ=255)
                                                                 gCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGt                                                                                  >  3:755096/1‑149 (MQ=255)
                                                                 gCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGt                                                                                  >  3:1092427/1‑149 (MQ=255)
                                                                     ttGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAg                                                                              >  3:1671859/1‑149 (MQ=255)
                                                                      tGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAgg                                                                             <  2:79351/149‑1 (MQ=255)
                                                                         gcAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGa                                                                          >  3:84362/1‑149 (MQ=255)
                                                                           aaGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGAcc                                                                        >  4:140095/1‑149 (MQ=255)
                                                                                ttCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgctg                                                                   >  2:327843/1‑149 (MQ=255)
                                                                                    aGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATa                                                               <  4:1515524/149‑1 (MQ=255)
                                                                                     gTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATaa                                                              >  4:44607/1‑149 (MQ=255)
                                                                                     gTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATaa                                                              <  1:799728/149‑1 (MQ=255)
                                                                                      tGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATaaa                                                             <  1:1127171/149‑1 (MQ=255)
                                                                                      tGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATaaa                                                             >  4:1462472/1‑149 (MQ=255)
                                                                                      tGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATaaa                                                             <  3:1310713/149‑1 (MQ=255)
                                                                                         gTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgc                                                                     <  3:450002/138‑1 (MQ=255)
                                                                                         gTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAacac                                                          >  1:1518617/1‑149 (MQ=255)
                                                                                          tCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTGATTGGTCAGCACCAGCTACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACAcc                                                         >  2:1077619/1‑149 (MQ=255)
                                                                                                   ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCgg                                                >  3:420189/1‑149 (MQ=255)
                                                                                                   ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCgg                                                >  2:190633/1‑149 (MQ=255)
                                                                                                    gAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                               >  3:355306/1‑149 (MQ=255)
                                                                                                    gAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                               >  1:1247378/1‑149 (MQ=255)
                                                                                                    gAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                               >  1:1002522/1‑149 (MQ=255)
                                                                                                     aaTACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGa                                              >  3:956175/1‑149 (MQ=255)
                                                                                                     aaTACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGa                                              <  4:84362/149‑1 (MQ=255)
                                                                                                     aaTACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGa                                              <  4:520/149‑1 (MQ=255)
                                                                                                      aTACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGAt                                             <  1:440933/149‑1 (MQ=255)
                                                                                                      aTACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGAt                                             <  4:1671859/149‑1 (MQ=255)
                                                                                                          cTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAg                                                                              >  1:495109/1‑112 (MQ=255)
                                                                                                          cTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAg                                                                              <  2:495109/112‑1 (MQ=255)
                                                                                                          cTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAAt                                         >  3:553912/1‑149 (MQ=255)
                                                                                                           ttCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGaa                                                    <  1:165099/137‑1 (MQ=255)
                                                                                                           ttCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAAtt                                        >  1:604036/1‑149 (MQ=255)
                                                                                                           ttCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAAtt                                        >  2:849646/1‑149 (MQ=255)
                                                                                                              cACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCtt                                     <  2:663926/149‑1 (MQ=255)
                                                                                                               aCCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCttt                                    <  1:1714734/149‑1 (MQ=255)
                                                                                                               aCCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCttt                                    <  3:216770/149‑1 (MQ=255)
                                                                                                               aCCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCttt                                    <  2:132342/149‑1 (MQ=255)
                                                                                                                 cTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTc                                       >  3:812208/1‑144 (MQ=255)
                                                                                                                  tGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTaa                                 >  1:406881/1‑149 (MQ=255)
                                                                                                                   gCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAAt                                >  1:685420/1‑149 (MQ=255)
                                                                                                                     acaTAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGt                              >  4:1639475/1‑149 (MQ=255)
                                                                                                                       ataataGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGttt                            <  2:923970/149‑1 (MQ=255)
                                                                                                                        taataGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAacac                                                          <  1:1468250/118‑1 (MQ=255)
                                                                                                                        taataGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGtttt                           <  1:1080660/149‑1 (MQ=255)
                                                                                                                        taataGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGtttt                           >  3:1631067/1‑149 (MQ=255)
                                                                                                                        taataGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGtttt                           <  1:190633/149‑1 (MQ=255)
                                                                                                                           taGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCg                        >  3:1016269/1‑149 (MQ=255)
                                                                                                                           taGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCg                        >  2:717694/1‑149 (MQ=255)
                                                                                                                           taGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCg                        >  2:1717568/1‑149 (MQ=255)
                                                                                                                                  ttCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCt                 <  4:1299895/149‑1 (MQ=255)
                                                                                                                                  ttCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCt                 <  2:1247378/149‑1 (MQ=255)
                                                                                                                                       gcgCGTTTAACCCAGCCGGG‑TAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATc            >  1:746447/3‑149 (MQ=255)
                                                                                                                                       gcgCGTTTAACCCAGCCGGG‑TAAACGCACCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATc            >  1:475586/3‑149 (MQ=255)
                                                                                                                                       gTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGtttt                           <  3:1161854/134‑1 (MQ=255)
                                                                                                                                       gTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATc            >  4:1340859/1‑149 (MQ=255)
                                                                                                                                       gTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATc            >  2:564082/1‑149 (MQ=255)
                                                                                                                                           gtttaaCCCAGCCGGG‑TAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAAt                                >  3:506031/6‑125 (MQ=255)
                                                                                                                                           gtttaaCCCAGCCGGG‑TAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAAt                                <  4:506031/120‑1 (MQ=255)
                                                                                                                                           gTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTg        >  3:711637/1‑149 (MQ=255)
                                                                                                                                            tttaaCCCAGCCGGG‑TAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCttt                                    >  4:95507/5‑120 (MQ=255)
                                                                                                                                            ttt‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTCTCGCGGTCAAGGCGACctgc                                                                     >  4:1167308/1‑87 (MQ=255)
                                                                                                                                            ttt‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAac                                                            <  1:1720321/96‑1 (MQ=255)
                                                                                                                                                aCCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACt    >  3:1516521/1‑149 (MQ=255)
                                                                                                                                                aCCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACt    >  1:1438277/1‑149 (MQ=255)
                                                                                                                                                 cccGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTg   >  3:308139/1‑149 (MQ=255)
                                                                                                                                                  ccGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCAcc                                                                                                        >  3:1454974/1‑47 (MQ=255)
                                                                                                                                                  ccGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGc  >  2:449424/1‑149 (MQ=255)
                                                                                                                                                  ccGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGc  <  2:608120/149‑1 (MQ=255)
                                                                                                                                                  ccGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGc  >  4:541316/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGC  >  NC_000913/1194503‑1194796

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑