Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,082 A→C 5.2% intergenic (+480/+134) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,270,0820AC 5.2% 239.0 / 2.7 92intergenic (+480/+134)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (43/40);  new base C (2/4);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGA  >  NC_000913/3269936‑3270227
                                                                                                                                                  |                                                                                                                                                 
tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGCAGCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGc                                                                                                                                                <  2:1457993/150‑1 (MQ=255)
      tGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGttt                                                                                                                                          <  2:543701/150‑1 (MQ=255)
      tGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGttt                                                                                                                                          <  2:543742/150‑1 (MQ=255)
       gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtttt                                                                                                                                         >  1:1662219/1‑150 (MQ=37)
        ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                                                                                                                                        <  2:370656/150‑1 (MQ=37)
        ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                                                                                                                                        <  2:1652644/150‑1 (MQ=37)
         cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                                       <  2:1494300/150‑1 (MQ=37)
             cccccaGTACACACCCGCCGCAGCGGGTTTTTACGTCCCACAGGTGCCCCTGCCCGATCAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                                                                                                   <  2:1583039/145‑1 (MQ=11)
              cAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCCGTCAGTCATCTGTCTGCCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGGTGTTACGTac                                                                                                                                  >  1:893415/1‑149 (MQ=11)
               aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                                                                                 >  1:205840/1‑150 (MQ=37)
                acaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAc                                                                                                                                <  2:1590590/148‑1 (MQ=38)
                aaaGTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAc                                                                                                                                <  2:1092367/150‑1 (MQ=38)
                aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAc                                                                                                                                <  1:956127/150‑1 (MQ=38)
                 aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAcc                                                                                                                               <  1:260889/150‑1 (MQ=38)
                  aGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCa                                                                                                                              <  1:717725/150‑1 (MQ=39)
                    tAAAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGACATTAATATTTATGCCCCTCCAAAAACAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                                                                                                            >  1:1367374/1‑150 (MQ=11)
                    tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTTCGGAAAACAgg                                                                                                                            >  1:151654/1‑150 (MQ=255)
                     aaaaaCCCGCCGAAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGt                                                                                                                           <  1:18296/150‑1 (MQ=37)
                     aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGt                                                                                                                           <  1:18297/150‑1 (MQ=37)
                      aaaaCCCGCTTCGGCGGGTTTTTACTTAAAACAGGTGAAACTGACCGATAAGCCGATATCTTTTGGGTATCGTGTCGTGGACAGTCATTCATCTTTCAGCCCCTCCAAAAGCCAAAACCCGCCGCAGCGGGGTGTTACGGACACCCgggg                                                                                                                          >  2:1243076/1‑148 (MQ=11)
                         aCCCGCCGAAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaaa                                                                                                                       <  1:729677/150‑1 (MQ=37)
                         aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaaa                                                                                                                       <  1:1323397/150‑1 (MQ=37)
                         aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAACCAAAAACCCGCCGAAGCGCGGTGTTAAGTAAACCAGGTGaaa                                                                                                                       >  1:364468/1‑150 (MQ=11)
                         aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTAGAGTGTCGTGGAAAGTCATTCATCTTTAGGCCCCGACAAAAGCAAAAACCCGCCGAAGCGGGTTGTTACGTAAAACAggggcaa                                                                                                                       >  1:1301345/1‑145 (MQ=11)
                         aCCCGCCGAAGCGGGATTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaaa                                                                                                                       <  1:1466061/150‑1 (MQ=37)
                           ccgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACt                                                                                                                     >  1:1267509/1‑150 (MQ=37)
                           ccgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACt                                                                                                                     >  1:607811/1‑150 (MQ=37)
                           ccgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCAGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACt                                                                                                                     >  2:1051181/1‑150 (MQ=21)
                             gccgCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGGGGCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTACGTAAACCAGGTGAAACTGAc                                                                                                                  <  1:278241/150‑1 (MQ=255)
                              ccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAc                                                                                                                  <  1:1565960/150‑1 (MQ=37)
                              ccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAc                                                                                                                  <  1:595775/150‑1 (MQ=37)
                                       gTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGcc                                                                                                         >  1:933780/1‑150 (MQ=21)
                                       gTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGcc                                                                                                         >  1:1026655/1‑150 (MQ=21)
                                       gTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGcc                                                                                                         >  2:93282/1‑150 (MQ=21)
                                        tttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCg                                                                                                        >  1:1571561/1‑150 (MQ=37)
                                               gTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGcggtctg                                                                                                 >  2:828361/1‑144 (MQ=255)
                                                 acatcAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                               <  2:196802/146‑1 (MQ=14)
                                                 aaaaCAGTTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCTCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                               <  2:196801/150‑1 (MQ=255)
                                                  aatcAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                              <  2:1008520/147‑1 (MQ=17)
                                                  aaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                              >  1:240223/1‑150 (MQ=21)
                                                  aaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAAAGTCATTCATCTTTCTGCCCCTACAAAAGCAAAAACCCGCCGAAGCGGGTTGTTACGTAAACCAGGTGAAACTGAACGATAAGACGCTTTCTTTTg                                                                                              >  2:680388/1‑150 (MQ=14)
                                                   atcAGGTGACACTGACCGATAAGCCGCGTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAGCTGTCTGCCCCTCCACAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                             <  2:362852/148‑1 (MQ=11)
                                                   aaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                             <  2:663990/150‑1 (MQ=38)
                                                    aCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg                                                                                            <  1:345409/150‑1 (MQ=38)
                                                    aCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCAACAAAATCAAAAAACCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTTTTggg                                                                                            >  1:492207/1‑150 (MQ=11)
                                                        gTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGGCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTTCGTAAACCAGGTGAAACAGTCCGATAAGCCGCTTTCTGTggggggga                                                                                        >  1:802032/1‑146 (MQ=11)
                                                             aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtc                                                                                   <  1:649801/150‑1 (MQ=37)
                                                               tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcgt                                                                                 <  2:1660512/150‑1 (MQ=37)
                                                               tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcgt                                                                                 <  2:1660567/150‑1 (MQ=37)
                                                               tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCATTCTTTTGGGTATAGcgtcgt                                                                                 >  2:660846/1‑150 (MQ=37)
                                                                    gATAAGCCTCTTTCTTTTGGGTATAGTGTCGTTGAAAGTCATTCATATTTCTGCCCCTCCAAAAGAAAAAACCCGCCGTATCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACa                                                                            <  2:989766/150‑1 (MQ=11)
                                                                    gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCCCAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACa                                                                            <  1:1623379/150‑1 (MQ=255)
                                                                    gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACa                                                                            <  2:243106/150‑1 (MQ=37)
                                                                       aaGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                                                         <  2:205840/150‑1 (MQ=37)
                                                                            gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattca                                                                    <  2:1216309/150‑1 (MQ=37)
                                                                            gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattca                                                                    <  2:1326453/150‑1 (MQ=37)
                                                                                 ctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCttt                                                               <  2:1662219/150‑1 (MQ=37)
                                                                                       tttttaGTGTCGTGGACAGTCATTGTTATTTATGGCACTCCACAATCAAAAACCCGCCGTAGCGTTTTTTTACGTACACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGccc                                                         <  1:481296/146‑1 (MQ=11)
                                                                                       gggctaGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGccc                                                         <  1:481515/146‑1 (MQ=37)
                                                                                              tgtCGTGGACAGTCATTCATCTTTATGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACAGACCGATAAGCCGCTTTCTTTTGGGGATAGCGTCGTGGACAGTCATTCATAATTAAGCCCCACCaaa                                                  >  1:1517046/1‑150 (MQ=11)
                                                                                               gtCGTGGACAGTCTTTCATCTTTCGGCGCTTGCACGAGCTAAAACCCGCCGATGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                 <  2:210779/150‑1 (MQ=11)
                                                                                               gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGTAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                 <  2:569929/150‑1 (MQ=37)
                                                                                                tCGTGGACAGTGAGTCATCTTTCTGGCCCTCCAAAAGCATAAACCCGCCGAATCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAg                                                <  2:63440/150‑1 (MQ=11)
                                                                                                tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGTTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAg                                                <  2:981064/150‑1 (MQ=37)
                                                                                                tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAg                                                >  1:716846/1‑150 (MQ=37)
                                                                                                tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAg                                                >  1:212728/1‑150 (MQ=37)
                                                                                                tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAg                                                >  1:212732/1‑150 (MQ=37)
                                                                                                tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAATCATCATTCTGCCCCTCCAAAAg                                                >  2:190216/1‑150 (MQ=255)
                                                                                                    ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAATCATCTTAATGCCCCGACAAAACCaaa                                            >  2:445774/1‑150 (MQ=11)
                                                                                                        aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc                                        >  2:1011216/1‑150 (MQ=37)
                                                                                                        aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAAcc                                        >  1:742712/1‑150 (MQ=37)
                                                                                                         gTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAccc                                       >  1:766815/1‑150 (MQ=37)
                                                                                                          tcattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGTGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccg                                      <  1:1473936/150‑1 (MQ=37)
                                                                                                             ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCGCCAAAAGAAAAAACccgccg                                   >  2:449912/1‑150 (MQ=255)
                                                                                                             ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGAGTCTATGAAAGTCATTCAACTTACATCACCACCAAAAGCAAAAACccgcag                                   >  2:1017551/1‑148 (MQ=11)
                                                                                                                atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                                >  1:156627/1‑150 (MQ=37)
                                                                                                                atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAACCGCCGAAg                                >  1:82961/1‑150 (MQ=37)
                                                                                                                atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCAACTTTATGCCCCACCAAAAGCAAAAACCCGCCGAAg                                >  2:1223443/1‑150 (MQ=255)
                                                                                                                atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGCATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAATAACCCGCCGAAg                                >  1:509443/1‑150 (MQ=255)
                                                                                                                   tttCTGCCCCTCAAAAAGCTAAAACCCGCCGAAGCGGGTTATTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCATACAAAAGCAAAAACCCGCCGAAGCgg                             <  1:671000/150‑1 (MQ=11)
                                                                                                                         ccccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                       >  2:1137944/1‑150 (MQ=37)
                                                                                                                            cTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                    >  2:337279/1‑150 (MQ=37)
                                                                                                                            cTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                    >  2:337276/1‑150 (MQ=37)
                                                                                                                              ccAAAAGCAAAATCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATTGCGTCGTGGACAGTCATTCATCTTTCTGCCCCGCCAAAAGCAATAACCCGCCGAAGCGGGTTTTTACGtt                  >  2:1272459/1‑149 (MQ=255)
                                                                                                                                aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAAGCAAAAACCCGCCCAAGCGGGTTTTTACGTaaa                >  2:595788/1‑150 (MQ=255)
                                                                                                                                   cgtaaaaaCCCGCTTCGGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGGTTTACGGAAATCa             >  2:949723/4‑150 (MQ=11)
                                                                                                                                     cAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCCAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg           >  1:313197/1‑150 (MQ=38)
                                                                                                                                     aaaaaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATACGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTTCGTAAATCAgg           <  1:1334344/149‑1 (MQ=255)
                                                                                                                                      aaaaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGt          >  1:1114734/1‑150 (MQ=37)
                                                                                                                                          aCCCGCCGAAGCGGGTTTTTATGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGaaa      <  1:13897/150‑1 (MQ=37)
                                                                                                                                            ccgccgCAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACt    <  1:63957/150‑1 (MQ=37)
                                                                                                                                              gccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGa  >  2:1705076/1‑150 (MQ=37)
                                                                                                                                                  |                                                                                                                                                 
TCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGA  >  NC_000913/3269936‑3270227

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: