Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 741,501 G→T 7.5% intergenic (+98/‑402) PP_t12 → / → PP_0635 tRNA‑Thr/group II intron‑encoding maturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947741,5010GT 7.5% 203.8 / 5.5 67intergenic (+98/‑402)PP_t12/PP_0635tRNA‑Thr/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base G (34/28);  new base T (2/3);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

CGTAACGCGAAGGTCGCAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTGACC  >  NC_002947/741361‑741613
                                                                                                                                            |                                                                                                                
cGTAACGCGAAGGTCGCAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCaa                                                                                                         <  2:713329/150‑1 (MQ=255)
 gTAACGCGAAGGTCGCAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCTTACGCGCAAc                                                                                                        <  2:1004365/150‑1 (MQ=255)
 gTAACGCGAAGGTCGCAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAAc                                                                                                        >  2:401733/1‑150 (MQ=255)
  tAACGCGAAGGTCGCAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAAcc                                                                                                       >  1:570668/1‑150 (MQ=255)
      gcgAAGGTCGCAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCTTACGCGCaa                                                                                                         <  2:183775/144‑1 (MQ=255)
      gcgAAGGTCGCAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCTTACGCGCaa                                                                                                         >  1:183775/1‑144 (MQ=255)
             tcacAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTtgt                                                                                            <  2:1120111/147‑1 (MQ=255)
               gCAGGTTCGATTGCAGACTCTGGCACCAGAATCAGTTAGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAATCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCAAAGCATTGCGCGTTGCGCGTAAGCGCGACCCGCTTGGTTgtgg                                                                                          <  2:236322/150‑1 (MQ=255)
               gCAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCACCCCGCTTGGTTgtgg                                                                                          >  1:236322/1‑150 (MQ=255)
               gCAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTgtgg                                                                                          >  1:497927/1‑150 (MQ=255)
                 aGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTgtggtg                                                                                        <  2:569236/150‑1 (MQ=255)
                   gTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGc                                                                                      >  2:1119312/1‑150 (MQ=255)
                       gATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAg                                                                                  >  1:894922/1‑150 (MQ=255)
                          tCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                               >  1:990089/1‑150 (MQ=255)
                            cTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                             <  2:497927/150‑1 (MQ=255)
                            cTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                             <  2:1038127/150‑1 (MQ=255)
                            cTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCACTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                             >  1:1212003/1‑150 (MQ=255)
                                   cGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTgg                                                                      >  1:384951/1‑150 (MQ=255)
                                     gCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                                    <  1:401733/150‑1 (MQ=255)
                                       aCCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGt                                                                  <  1:191278/150‑1 (MQ=255)
                                        ccAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGgtggtg                                                                 >  1:979529/1‑150 (MQ=255)
                                         cAGAATCAGTTTGGCTATGCCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGa                                                                <  2:790640/150‑1 (MQ=255)
                                         cAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTgg                                                                <  1:1119312/150‑2 (MQ=255)
                                         cAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGa                                                                >  1:422585/1‑150 (MQ=255)
                                           gAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                              >  2:1234373/1‑150 (MQ=255)
                                           gAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                              <  1:706900/150‑1 (MQ=255)
                                              tCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGGCTTGTGCGGGCTTTTTTGTGTGCGACAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTc                                                           >  1:413930/1‑150 (MQ=255)
                                                 gTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGCAAGTcccc                                                        <  2:844413/150‑3 (MQ=255)
                                                 gTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                                        <  1:500635/150‑1 (MQ=255)
                                                     gTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTaca                                                    >  1:1339034/1‑150 (MQ=255)
                                                       cTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacaca                                                  >  1:334745/1‑150 (MQ=255)
                                                        tATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacacac                                                 <  1:1215381/150‑1 (MQ=255)
                                                          tGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                               <  2:1075122/118‑1 (MQ=255)
                                                          tGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                               >  1:1075122/1‑118 (MQ=255)
                                                            gCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGt                                                                                 >  2:677024/1‑114 (MQ=255)
                                                            gCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGt                                                                                 <  1:677024/114‑1 (MQ=255)
                                                            gCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCgg                                             >  1:1012898/1‑150 (MQ=255)
                                                                    tcGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCTTACGCGCAAcgc                                                                                                      <  1:745090/85‑3 (MQ=25)
                                                                    tcGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCTTACGCGCAAcgc                                                                                                      >  2:745090/1‑83 (MQ=25)
                                                                    tcGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGa                                     >  1:1372681/1‑150 (MQ=255)
                                                                      gCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAg                                   <  1:1234373/150‑1 (MQ=255)
                                                                      gCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAg                                   >  1:507109/1‑150 (MQ=255)
                                                                        tACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgtg                                 >  2:209042/1‑150 (MQ=255)
                                                                         aCCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgtgt                                >  2:83650/1‑150 (MQ=255)
                                                                                aCCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCgg                         >  2:459899/1‑150 (MQ=255)
                                                                                aCCAATAGAAGCCCGCCTTGTGCGGGCTTTCTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCgg                         >  2:426345/1‑150 (MQ=255)
                                                                                  cAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg                       <  1:459899/150‑1 (MQ=255)
                                                                                  cAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGCGTTAGCCGGAg                       <  2:422585/150‑1 (MQ=255)
                                                                                   aaTAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAgg                      >  1:1229422/1‑150 (MQ=255)
                                                                                   aaTAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAgg                      <  2:507109/150‑1 (MQ=255)
                                                                                   aaTAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAgg                      <  2:990089/150‑1 (MQ=255)
                                                                                    aTAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGc                     >  2:1388157/1‑150 (MQ=255)
                                                                                     tAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCa                    >  2:259377/1‑150 (MQ=255)
                                                                                       gAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAg                  <  1:916973/150‑1 (MQ=255)
                                                                                       gAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAACCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAg                  <  1:209042/150‑1 (MQ=255)
                                                                                        aaGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGTGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAgg                 <  2:979529/150‑1 (MQ=255)
                                                                                          gCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGT               >  2:1315358‑M1/1‑150 (MQ=255)
                                                                                           cccGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTG              >  1:890211‑M1/1‑150 (MQ=255)
                                                                                            ccGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGTTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGgtgt             >  2:343909/1‑150 (MQ=38)
                                                                                                 ttGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCgg        <  2:334745/150‑1 (MQ=34)
                                                                                                 ttGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCgg        >  2:271384/1‑150 (MQ=34)
                                                                                                   gtgCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGt      <  1:363913/150‑1 (MQ=32)
                                                                                                      cGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCct                                                         <  2:569252/96‑1 (MQ=32)
                                                                                                      cGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCct                                                         >  1:569252/1‑96 (MQ=32)
                                                                                                       gggCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTGTACTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTGAcc  >  1:375039/1‑150 (MQ=12)
                                                                                                                                      ttGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGT               >  2:705356‑M1/1‑106 (MQ=255)
                                                                                                                                      ttGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGT               <  1:705356‑M1/106‑1 (MQ=255)
                                                                                                                                            |                                                                                                                
CGTAACGCGAAGGTCGCAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTGACC  >  NC_002947/741361‑741613

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