Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,082,709 C→A 5.0% intergenic (+38/+168) loiP → / ← speB Phe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; Era‑binding heat shock protein/agmatinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,082,7090CA 5.0% 1107.8 / 31.2 298intergenic (+38/+168)loiP/speBPhe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; Era‑binding heat shock protein/agmatinase
Reads supporting (aligned to +/- strand):  ref base C (150/133);  new base A (6/9);  total (156/142)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 308 total aligned reads displayed.

CTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGT  >  NC_000913/3082561‑3082858
                                                                                                                                                    |                                                                                                                                                     
 tGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAc                                                                                                                                                     >  2:2655427/1‑150 (MQ=255)
 tGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAc                                                                                                                                                     >  2:2468370/1‑150 (MQ=255)
 tGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAc                                                                                                                                                     >  1:1977368/1‑150 (MQ=255)
   gCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAcgc                                                                                                                                                   >  1:3819668/1‑150 (MQ=255)
   gCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAcgc                                                                                                                                                   <  1:2128714/150‑1 (MQ=255)
   gCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAcgc                                                                                                                                                   >  2:894106/1‑150 (MQ=255)
     aaaaCTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAcgccg                                                                                                                                                 >  1:1481548/1‑150 (MQ=255)
     aaaaCTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAcgccg                                                                                                                                                 >  2:3622255/1‑150 (MQ=255)
     aaaaCTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGATGTCGAACGTAGGTCAGATAAGGCGTTCAcgccg                                                                                                                                                 >  2:253101/1‑150 (MQ=255)
       aaCTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCa                                                                                                                                               >  1:3863064/1‑150 (MQ=255)
        aCTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATGCGGCGTGAAcgccttat                                                                                                                                              <  2:1928892/150‑5 (MQ=255)
                    tCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATg                                                                                                                                  <  1:623153/150‑1 (MQ=255)
                     cGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATgg                                                                                                                                 <  2:62030/150‑1 (MQ=255)
                     cGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATgg                                                                                                                                 <  2:350554/150‑1 (MQ=255)
                      gCCAAAGCTCAATGTTTGACGACCACCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCTTGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGc                                                                                                                                <  2:3433018/150‑1 (MQ=255)
                      gCCAAAGCTCAATGTTTGACGACCACCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGc                                                                                                                                <  2:2198961/150‑1 (MQ=255)
                        cAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGcac                                                                                                                              >  1:3395898/1‑150 (MQ=255)
                         aaaGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGCTCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGcaca                                                                                                                             >  1:2959521/1‑150 (MQ=255)
                         aaaGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGcaca                                                                                                                             <  2:2381927/150‑1 (MQ=255)
                         aaaGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGcaca                                                                                                                             <  1:3680346/150‑1 (MQ=255)
                            gCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACAttt                                                                                                                          <  2:3882605/150‑1 (MQ=255)
                              tCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGt                                                                                                                        <  1:1623488/150‑1 (MQ=255)
                               cAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGtt                                                                                                                       >  1:722648/1‑150 (MQ=255)
                                  tGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCAAATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTAcc                                                                                                                    <  2:2110604/150‑1 (MQ=255)
                                   gTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCt                                                                                                                   >  2:1573616/1‑150 (MQ=255)
                                   gTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCt                                                                                                                   >  2:1970751/1‑150 (MQ=255)
                                    tttGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCtt                                                                                                                  >  1:4103844/1‑150 (MQ=255)
                                    tttGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGCGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGACAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCtt                                                                                                                  <  2:2287999/150‑1 (MQ=255)
                                         cgacCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTgcg                                                                                                             <  2:2882918/150‑1 (MQ=255)
                                            ccATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGCTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCAt                                                                                                          >  2:1394875/1‑150 (MQ=255)
                                                 cTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCg                                                                                                     <  1:1515868/150‑1 (MQ=255)
                                                 cTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCg                                                                                                     <  1:1893690/150‑1 (MQ=255)
                                                   gCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGga                                                                                                   <  1:2577220/150‑1 (MQ=255)
                                                    cATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgat                                                                                                  <  1:3485326/150‑1 (MQ=255)
                                                      tCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgc                                                                                                >  1:3794117/1‑150 (MQ=255)
                                                          ccgAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATTTGATTTAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgcgatg                                                                                            <  2:3346042/150‑1 (MQ=255)
                                                           cgAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgcgatgc                                                                                           >  1:3754196/1‑150 (MQ=255)
                                                               cgcgCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGgc                                                                                       <  1:354322/150‑1 (MQ=255)
                                                                gcgcCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGgcg                                                                                      <  1:3118419/150‑1 (MQ=255)
                                                                  gcCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCa                                                                                    <  1:1181953/150‑1 (MQ=255)
                                                                            ttCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCgg                                                                          >  1:2917885/1‑150 (MQ=255)
                                                                             tCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGc                                                                         <  1:3324458/150‑1 (MQ=255)
                                                                             tCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGc                                                                         <  2:1000586/150‑1 (MQ=255)
                                                                              cgcgATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGcc                                                                        <  2:2611806/150‑1 (MQ=255)
                                                                                 gATCGCATGAGCGCGGATGGTATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTAc                                                                     <  1:2286654/150‑1 (MQ=255)
                                                                                   tCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACgg                                                                   <  2:513982/150‑1 (MQ=255)
                                                                                     gCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGt                                                                 >  1:1368109/1‑150 (MQ=255)
                                                                                       aTGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGc                                                               >  2:3149709/1‑150 (MQ=255)
                                                                                          aGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGa                                                            <  1:481007/150‑1 (MQ=255)
                                                                                            cgcgGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAAc                                                          <  1:907344/150‑1 (MQ=255)
                                                                                              cgGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGACTGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGt                                                        <  1:1657177/150‑1 (MQ=255)
                                                                                              cgGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGt                                                        <  1:816061/150‑1 (MQ=255)
                                                                                               ggATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTa                                                       >  1:3639345/1‑150 (MQ=255)
                                                                                                 aTGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAgg                                                     >  2:2827803/1‑150 (MQ=255)
                                                                                                   ggggTTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTc                                                   <  2:181107/150‑1 (MQ=255)
                                                                                                   gggATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTc                                                   <  2:1132682/150‑1 (MQ=255)
                                                                                                   gggATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTc                                                   <  2:749587/150‑1 (MQ=255)
                                                                                                     gATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCgg                                                 >  2:3043926/1‑150 (MQ=255)
                                                                                                     gATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCgg                                                 >  2:2972789/1‑150 (MQ=255)
                                                                                                      aTTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGa                                                >  2:3766012/1‑150 (MQ=255)
                                                                                                       ttAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGAt                                               >  2:3694943/1‑150 (MQ=255)
                                                                                                          aGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAg                                            <  1:3946994/150‑1 (MQ=255)
                                                                                                          aGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAg                                            >  1:2250536/1‑150 (MQ=255)
                                                                                                           gTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAgg                                           >  2:1752256/1‑150 (MQ=255)
                                                                                                           gTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAgg                                           <  2:3207177/150‑1 (MQ=255)
                                                                                                           gTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAgg                                           >  1:2449399/1‑150 (MQ=255)
                                                                                                            tAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGc                                          <  1:3207312/150‑1 (MQ=255)
                                                                                                              aaGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGt                                        >  1:2657511/1‑150 (MQ=255)
                                                                                                              aaGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGt                                        <  1:1443624/150‑1 (MQ=255)
                                                                                                              aaGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGt                                        <  1:1092653/150‑1 (MQ=255)
                                                                                                                       ggtgTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgccgc                               <  1:1154885/150‑1 (MQ=255)
                                                                                                                          gtCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATc                            <  2:4169826/150‑1 (MQ=255)
                                                                                                                             gAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                         >  2:1873862/1‑150 (MQ=255)
                                                                                                                             gAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                         <  2:1496553/150‑1 (MQ=255)
                                                                                                                             gAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                         <  2:406628/150‑1 (MQ=255)
                                                                                                                             gAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                         <  2:3794138/150‑1 (MQ=255)
                                                                                                                              aaCGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGAc                        >  1:767560/1‑150 (MQ=255)
                                                                                                                               aCGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCAGTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGATCACGCCGCATCCGACa                       >  2:3643230/1‑150 (MQ=255)
                                                                                                                                 gTAGGTCAGATAAGGCGTTCACGCCGCATCTGTCATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATCAGGCGTTCACGCCGCATCCGACagt                     >  2:2766712/1‑148 (MQ=255)
                                                                                                                                    tgTCAGATGCGGCGTGAACGCCTTATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAAt                  >  2:913131/2‑150 (MQ=38)
                                                                                                                                    tgTCAGATGCGGCGTGAACGCCTTATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAAt                  >  2:653646/2‑150 (MQ=38)
                                                                                                                                     gTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATg                 >  2:814140/1‑150 (MQ=255)
                                                                                                                                     gTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATg                 >  2:108624/1‑150 (MQ=255)
                                                                                                                                     gTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATg                 >  2:3240535/1‑150 (MQ=255)
                                                                                                                                      tCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATgg                >  1:1768055/1‑150 (MQ=255)
                                                                                                                                      tCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATgg                <  2:1164148/150‑1 (MQ=255)
                                                                                                                                      tCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATgg                <  2:1170190/150‑1 (MQ=255)
                                                                                                                                      tCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATgg                >  2:1661427/1‑150 (MQ=255)
                                                                                                                                      tCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATgg                >  2:479737/1‑150 (MQ=255)
                                                                                                                                       cAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGc               <  1:1573616/150‑1 (MQ=255)
                                                                                                                                       cAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGc               <  1:1970650/150‑1 (MQ=255)
                                                                                                                                          aTAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACg            >  1:299358/1‑150 (MQ=255)
                                                                                                                                            aaGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtt          >  1:4046161/1‑150 (MQ=255)
                                                                                                                                            aaGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtt          <  1:2843342/150‑1 (MQ=255)
                                                                                                                                            aaGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtt          <  1:496832/150‑1 (MQ=255)
                                                                                                                                            aaGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATGAATGGCACGtt          >  1:2673733/1‑150 (MQ=255)
                                                                                                                                               gCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTa       <  2:2013875/150‑1 (MQ=255)
                                                                                                                                                  ttCACGCCGCATCTGACATGAATGGCACAGTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGAGAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTAccc    >  2:3570506/1‑150 (MQ=255)
                                                                                                                                                   tCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCg   >  2:3999860/1‑150 (MQ=255)
                                                                                                                                                    cACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGt  <  1:95713/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGT  >  NC_000913/3082561‑3082858

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