Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,052,980 C→T 5.2% L58L (CTG→CTA)  PP_3572 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,052,9800CT 5.2% 225.3 / 2.1 76L58L (CTG→CTA) PP_3572hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/41);  new base T (2/2);  total (33/43)
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.57e-01

TTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATACCGACCTGGCCGGTTTCCTGGCA  >  NC_002947/4052832‑4053129
                                                                                                                                                    |                                                                                                                                                     
ttGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCa                                                                                                                                                      >  1:1916320/1‑150 (MQ=255)
       gtgtgtGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATc                                                                                                                                               <  2:634972/150‑1 (MQ=255)
            tgCCCAGTGCTTCACGCCCGGTGAACAGGGCCACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGaa                                                                                                                                          <  1:1832570/150‑1 (MQ=255)
             gCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGCGAAt                                                                                                                                         >  2:328951/1‑150 (MQ=255)
               ccAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCt                                                                                                                                       >  2:57296/1‑150 (MQ=255)
                   tGCTTCCCGCCCGGTGAACAGGGCCACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCgg                                                                                                                                   <  2:1915163/150‑1 (MQ=255)
                   tGCTTCCCGCCCGGTGAACAGGGCCACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCgg                                                                                                                                   <  1:1698980/150‑1 (MQ=255)
                      ttCCCGCCCGGTGAACAGGGCCACGTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGcc                                                                                                                                >  1:58324/1‑150 (MQ=255)
                      ttCCCGCCCGGTGAACAGGGCCACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGcc                                                                                                                                >  2:1402714/1‑150 (MQ=255)
                         ccgcccgGTGAACAGGGCCACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAg                                                                                                                             >  2:1714266/1‑150 (MQ=255)
                         ccgcccgGTGAACAGGGCCACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAg                                                                                                                             <  2:793984/150‑1 (MQ=255)
                                 tGAACAGGGCCACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCa                                                                                                                     >  2:315722/1‑150 (MQ=255)
                                   aaCAGGGCCACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAgg                                                                                                                   <  1:600223/150‑1 (MQ=255)
                                    agaGGGCCACCGGGCCGGGGCTGTCGATCACCGTGACGTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGTGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAggg                                                                                                                  <  2:166982/148‑1 (MQ=255)
                                       gggCCACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGa                                                                                                               <  1:1329540/150‑1 (MQ=255)
                                           cACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGtt                                                                                                           <  1:1692511/150‑1 (MQ=255)
                                           cACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGtt                                                                                                           <  1:146834/150‑1 (MQ=255)
                                            aCCTGGCCGTGGCTGTAGATCACTGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTc                                                                                                          <  1:1038579/150‑1 (MQ=255)
                                              cTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTg                                                                                                        <  1:828524/150‑1 (MQ=255)
                                              cTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTATCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTg                                                                                                        <  1:1732464/150‑1 (MQ=255)
                                               tGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTgg                                                                                                       <  2:1788545/150‑1 (MQ=255)
                                                  ccGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTc                                                                                                    <  1:315722/150‑1 (MQ=255)
                                                    gTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCga                                                                                                  <  1:836513/150‑1 (MQ=255)
                                                        cTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCCGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGGATCTGCGGGCCCATCGCCGGCCCGGTTCGCTTCTGGGTCGAGCCt                                                                                              >  2:142558/1‑150 (MQ=255)
                                                         tGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTc                                                                                                    <  2:559774/143‑1 (MQ=255)
                                                         tGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTc                                                                                                    >  1:559774/1‑143 (MQ=255)
                                                             gATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTg                                                                                                        <  1:1321991/135‑1 (MQ=255)
                                                             gATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTg                                                                                                        >  2:1321991/1‑135 (MQ=255)
                                                              aTCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGccc                                                                                        >  1:1007026/1‑150 (MQ=255)
                                                                 aCCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCcaca                                                                                     >  2:531107/1‑150 (MQ=255)
                                                                  ccGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTg                                                                                                                <  1:1481977/122‑1 (MQ=255)
                                                                  ccGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTg                                                                                                                >  2:1481977/1‑122 (MQ=255)
                                                                           ggCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGgcacgca                                                                           <  2:1432865/150‑1 (MQ=255)
                                                                                tACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAg                                                                      <  1:1598576/150‑1 (MQ=255)
                                                                                 aCTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCTAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCgaga                                                                                                <  1:1703743/123‑1 (MQ=255)
                                                                                 aCTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCTAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCgaga                                                                                                <  1:1703758/123‑1 (MQ=255)
                                                                                 aCTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCTAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCgaga                                                                                                >  2:1703758/1‑123 (MQ=255)
                                                                                 aCTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCTAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCgaga                                                                                                >  2:1703743/1‑123 (MQ=255)
                                                                                 aCTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAgg                                                                     <  1:1404299/150‑1 (MQ=255)
                                                                                  cTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAggg                                                                    <  2:706418/150‑1 (MQ=255)
                                                                                  cTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAggg                                                                    >  1:1820450/1‑150 (MQ=255)
                                                                                  cTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAggg                                                                    <  1:1714266/150‑1 (MQ=255)
                                                                                   tCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATAGATAATCTGCGGGCTCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAgggg                                                                   >  1:1176189/1‑150 (MQ=255)
                                                                                   tCGCTCCAGCCATTGGCACTGAGCACCCTGCCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAgggg                                                                   >  2:1116201/1‑150 (MQ=255)
                                                                                     gCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCa                                                                 >  1:1150815/1‑150 (MQ=255)
                                                                                         cAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGa                                                             <  1:57296/150‑1 (MQ=255)
                                                                                         cAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGa                                                             <  2:579706/150‑1 (MQ=255)
                                                                                          aGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGAGTCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAg                                                            >  1:320063/1‑150 (MQ=255)
                                                                                           gCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGc                                                           <  1:1779470/150‑1 (MQ=255)
                                                                                             cATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGACGGCACCTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCAc                                                         <  2:724226/150‑1 (MQ=255)
                                                                                              aTTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCAcc                                                        <  2:1916320/150‑1 (MQ=255)
                                                                                                  gCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGCTGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACg                                                    >  2:518262/1‑150 (MQ=255)
                                                                                                    aCTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACggt                                                  >  2:798091/1‑149 (MQ=255)
                                                                                                      tGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGcc                                                                                         >  2:971047/1‑109 (MQ=255)
                                                                                                      tGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGcc                                                                                         <  1:971047/109‑1 (MQ=255)
                                                                                                             ccTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATTCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTgcgc                                                                                          <  2:1516853/101‑1 (MQ=255)
                                                                                                             ccTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATTCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTgcgc                                                                                          >  1:1516853/1‑101 (MQ=255)
                                                                                                                 tCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAg                                     <  1:1854667/150‑1 (MQ=255)
                                                                                                                   cagcgcagcgGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCt                                   >  2:1719348/1‑150 (MQ=255)
                                                                                                                         agcgGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGc                             >  1:1072668/1‑150 (MQ=255)
                                                                                                                               gggcggCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATAcc                       >  1:91219/1‑150 (MQ=255)
                                                                                                                                 gcggcAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGTGCTGATGTCGATACCGa                     <  2:117952/150‑1 (MQ=255)
                                                                                                                                  cggcAACTGGCCATGCTCCGGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATACCGAc                    <  1:765396/150‑1 (MQ=255)
                                                                                                                                        cTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACTAGCTGATGGCGATACCGACCTGGcc              <  1:1586175/150‑1 (MQ=255)
                                                                                                                                           gCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATACCGACCTGGCCGGt           >  2:1091781/1‑150 (MQ=255)
                                                                                                                                             cATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAgg                                                                     >  2:1555542/1‑90 (MQ=255)
                                                                                                                                             cATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAgg                                                                     <  1:1555542/90‑1 (MQ=255)
                                                                                                                                             cATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATACCGAc                    >  2:124338/1‑139 (MQ=255)
                                                                                                                                             cATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATACCGAc                    <  1:124338/139‑1 (MQ=255)
                                                                                                                                              aTGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATACCGAc                    >  1:871145/1‑138 (MQ=255)
                                                                                                                                              aTGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATACCGAc                    <  2:871145/138‑1 (MQ=255)
                                                                                                                                                gCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATACCGACCTGGCCGGTTTCCt      <  1:767553/150‑1 (MQ=255)
                                                                                                                                                gCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATACCGACCTGGCCGGTTTCCt      <  2:1176189/150‑1 (MQ=255)
                                                                                                                                                 ctacAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATACCGACCTGGCCGGTTTCCTg     <  2:1322956/147‑1 (MQ=255)
                                                                                                                                                 cTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATCGCTATACCGACCTGGCCGGTTTCCTg     <  1:1402714/150‑1 (MQ=255)
                                                                                                                                                  tCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATCCCGACCTGGACGGTTCCCTgc    >  2:1541743/1‑149 (MQ=255)
                                                                                                                                                  tCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATACCGACCTGGCCGGTTTCCTgg    >  1:658935/1‑150 (MQ=255)
                                                                                                                                                   ccAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATACCGACCTGGCCGGTTTCCTGGc   >  1:1770665/1‑150 (MQ=255)
                                                                                                                                                    cAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATACCGACCTGGCCGGTTTCCTGGCa  >  2:947576/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGCCGTGGCTGTCGATCACCGTGACCTGGCGGTACTCGCTCCAGCCATTGGCACTGAGCACCCTGTCCAGCGCAGCGGCGGGCGGCAACTGGCCATGCTCCAGGGCATCGAGAATCTGCGGGCCCAGCGCCGGCAGGGTGATGTTCTGGGTCGAGACTGCGCCCACACCGGCACGCACCCAGGGGCAACGAGCACCCACGGCGATGCTCGACGAGCTGATGGCGATACCGACCTGGCCGGTTTCCTGGCA  >  NC_002947/4052832‑4053129

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