Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 493,804 G→T 5.4% E79* (GAA→TAA) ‡ PP_0406 → nucleotidyltransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947493,8040GT 5.4% 235.4 / 2.8 74E79* (GAA→TAA) ‡PP_0406nucleotidyltransferase family protein
Reads supporting (aligned to +/- strand):  ref base G (35/35);  new base T (2/2);  total (37/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 = 5.73e-01

GCCAGCCGCTGATCGAGTACCACTTGCGCGCCCTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGTTTGCAGGCTCCCCTGCAAGGCCTGGCTCACCTGGT  >  NC_002947/493658‑493952
                                                                                                                                                  |                                                                                                                                                    
gCCAGCCGCTGATCGAGTACCACTTGCGCGCCCTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAg                                                                                                                                                   <  1:1138495/150‑1 (MQ=255)
  cAGCCGCTGATCGAGTACCACTTGAAATCCGGGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGt                                                                                                                                                 >  1:1506693/1‑150 (MQ=255)
     ccGCTGATCGAGTACCACTTGCGCGCCCTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGaa                                                                                                                                              >  2:148289/1‑150 (MQ=255)
     ccGCTGATCGAGTACCACTTGCGCGCCCTGGCGGCAGCGGGCGTCTGCCTAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGaa                                                                                                                                              <  1:700463/150‑1 (MQ=255)
      cGCTGATCGAGTACCACTTGCGCGCCCTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAAc                                                                                                                                             >  1:1178503/1‑150 (MQ=255)
       gCTGATCGAGTACCACTTGCGCGCCCTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAAcc                                                                                                                                            >  2:1767121/1‑150 (MQ=255)
           aTCGAGTACCACTTGCGCGCCCTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTg                                                                                                                                        <  1:1328611/150‑1 (MQ=255)
           aTCGAGTACCACTTGCGCGCCCTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTg                                                                                                                                        >  2:253119/1‑150 (MQ=255)
                   ccACTTGCGCGCCCTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACcgg                                                                                                                                <  1:531798/150‑1 (MQ=255)
                   ccACTTGCGCGCCCTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACcgg                                                                                                                                >  2:1757485/1‑150 (MQ=255)
                    cACTTGCGCGCCCTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAAcc                                                                                                                                            <  1:1320262/137‑1 (MQ=255)
                    cACTTGCGCGCCCTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAAcc                                                                                                                                            >  2:1320262/1‑137 (MQ=255)
                           gcgcCCTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCa                                                                                                                        >  2:31681/1‑150 (MQ=255)
                             gcCCTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATc                                                                                                                      <  1:657961/150‑1 (MQ=255)
                                cTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAgg                                                                                                                                                  >  2:368870/1‑119 (MQ=255)
                                cTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAgg                                                                                                                                                  <  1:368870/119‑1 (MQ=255)
                                 tGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCa                                                                                                                  >  1:1636892/1‑150 (MQ=255)
                                  ggcggcAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGTCGACGGTAGCCGCTTCGGACTGAGCATTCGATATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCaa                                                                                                                 <  1:787916/150‑1 (MQ=255)
                                   gcggcAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAg                                                                                                                <  2:815815/150‑1 (MQ=255)
                                          cGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTGGCTGTTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGc                                                                                                         <  2:1036105/150‑1 (MQ=255)
                                          cGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGc                                                                                                         <  1:1428497/150‑1 (MQ=255)
                                           gggCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGc                                                                                                         >  2:66633/1‑149 (MQ=255)
                                           gggCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGc                                                                                                         <  1:66633/149‑1 (MQ=255)
                                             gCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTTCAGGTGaatag                                                                                                                                           >  1:1327697/1‑110 (MQ=255)
                                             gCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTTCAGGTGaatac                                                                                                                                           <  2:1327697/113‑4 (MQ=255)
                                                   ccGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTgg                                                                                                >  1:498326/1‑150 (MQ=255)
                                                     gAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGt                                                                                              <  1:874699/150‑1 (MQ=255)
                                                       aGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGAACATGTGGGCGACGGTAGCCGCGTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGa                                                                                            <  1:417556/150‑1 (MQ=255)
                                                          ggtgATCAACCATGCCTGGCTTGGACAGCAGATCGAAGACCATCTGGACACCGGTAACTGCTTCGGCCGGAGCATTCACTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATgc                                                                                         >  2:709611/1‑150 (MQ=255)
                                                             gATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGcc                                                                                      >  2:1271239/1‑150 (MQ=255)
                                                              aTCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCg                                                                                     <  1:31681/150‑1 (MQ=255)
                                                              aTCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCATCGGGTTGTTGGGTGATGCGCCg                                                                                     >  1:1069029/1‑150 (MQ=255)
                                                                  aCCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTcc                                                                                 >  2:1610876/1‑150 (MQ=255)
                                                                   ccATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCct                                                                                >  1:588849/1‑150 (MQ=255)
                                                                   ccATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCct                                                                                >  2:401031/1‑150 (MQ=255)
                                                                            gCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGaa                                                                       >  2:193675/1‑150 (MQ=255)
                                                                              ttGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATg                                                                     >  2:1674529/1‑150 (MQ=255)
                                                                               tGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATgg                                                                    <  2:1178503/150‑1 (MQ=255)
                                                                                 gCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCg                                                                  >  1:1171616/1‑150 (MQ=255)
                                                                                     gcagATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGAtgt                                                              <  1:1674529/150‑1 (MQ=255)
                                                                                       agATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGAtgtgt                                                            <  1:253119/150‑1 (MQ=255)
                                                                                        gATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGTTGCGCCGTTCCTGCTGGTGAATGGCGAtgtgtg                                                           <  1:1610876/150‑1 (MQ=255)
                                                                                         aTCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACcggcggc                                                                                                                            >  1:702966/1‑84 (MQ=255)
                                                                                         aTCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACcggcggc                                                                                                                            <  2:702966/84‑1 (MQ=255)
                                                                                         aTCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGc                                                                   <  1:39820/141‑1 (MQ=255)
                                                                                         aTCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGc                                                                   >  2:39820/1‑141 (MQ=255)
                                                                                          tCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGa                                                         <  1:1757485/150‑1 (MQ=255)
                                                                                            gAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTACTGGTGAATGGCGATGTGTGGAcc                                                       <  2:498326/150‑1 (MQ=255)
                                                                                                  cATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTAc                                                 >  2:984061/1‑150 (MQ=255)
                                                                                                   aTCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACg                                                >  2:72645/1‑150 (MQ=255)
                                                                                                    tCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGa                                               <  1:984061/150‑1 (MQ=255)
                                                                                                        ggCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGAGTATTGCTCAAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTc                                           <  2:1506693/150‑1 (MQ=255)
                                                                                                            aCGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTcgcgc                                       <  1:1277327/150‑1 (MQ=255)
                                                                                                              ggTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTc                                                                                                                   >  2:1117027/1‑72 (MQ=255)
                                                                                                              ggTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTc                                                                                                                   <  1:1117027/72‑1 (MQ=255)
                                                                                                                   ccGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTggg                                                                                               >  1:1883718/1‑87 (MQ=255)
                                                                                                                   ccGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTggg                                                                                               <  2:1883718/87‑1 (MQ=255)
                                                                                                                   ccGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGTTTGCa                                <  1:1767121/150‑1 (MQ=255)
                                                                                                                            ccTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAg                                                                                                                >  2:28418/1‑61 (MQ=255)
                                                                                                                            ccTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAg                                                                                                                <  1:28418/61‑1 (MQ=255)
                                                                                                                            ccTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGTTTGCAGGCTCCCCt                       >  2:1107186/1‑150 (MQ=255)
                                                                                                                              tGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGtt                                    >  2:493765/1‑135 (MQ=255)
                                                                                                                              tGAGCATTCGCTATGCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGtt                                    <  1:493765/135‑1 (MQ=255)
                                                                                                                               gAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGTTTGCAGGCTCCCCTGCa                    >  2:42911/1‑150 (MQ=255)
                                                                                                                               gAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGTTTGCAGGCTCCCCTGCa                    >  1:1767128/1‑150 (MQ=255)
                                                                                                                                  cATTCGCTATTCACCTGAAGGTGAACCACTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGTTTGCAGGCTCCCCTGCAAgg                 >  1:1495780/1‑150 (MQ=255)
                                                                                                                                    ttCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGAACGACTACGACTTCGCGCGTGTGCAGGCTCCCCTGCAAGGcc               >  2:1287089/1‑150 (MQ=255)
                                                                                                                                      agcggttCACCTTCAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGttt                                   <  2:977239/123‑1 (MQ=255)
                                                                                                                                      agcggttCACCTTCAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGttt                                   >  1:977239/6‑128 (MQ=255)
                                                                                                                                       gCTATTCACCTGACGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGTTTGCAGGCTCCCCTGCAAGGCCTgg            <  1:401031/150‑1 (MQ=255)
                                                                                                                                       gCTATTCACCTGAAGGTGCACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGTTTGCAGGCTCCCCTGCAAGGCCCgg            <  2:673220/150‑1 (MQ=255)
                                                                                                                                       gCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGTTTGCAGGCTCCCCTGCAAGGCCTgg            <  1:193675/150‑1 (MQ=255)
                                                                                                                                                cTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGTTTGCAGGCTCCCCTGCAAGGCCTGGCTCACCTgg   <  1:72645/150‑1 (MQ=255)
                                                                                                                                                 tGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGTTTGCAGGCTCCCCTGCAAGGCCTGGCTCACCTGGt  >  2:895103/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCAGCCGCTGATCGAGTACCACTTGCGCGCCCTGGCGGCAGCGGGCGTGACCGAAGTGGTGATCAACCATGCCTGGCTTGGCCAGCAGATCGAAGACCATCTGGGCGACGGTAGCCGCTTCGGCCTGAGCATTCGCTATTCACCTGAAGGTGAACCGCTGGAAACCGGCGGCGGCATCTTCAAGGCCCTGCCTTTGTTGGGTGATGCGCCGTTCCTGCTGGTGAATGGCGATGTGTGGACCGACTACGACTTCGCGCGTTTGCAGGCTCCCCTGCAAGGCCTGGCTCACCTGGT  >  NC_002947/493658‑493952

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: ‑