Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,821,936 A→T D107E (GAT→GAA)  osmE ← osmotically‑inducible lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,821,9360AT96.9% 240.6 / ‑0.6 65D107E (GAT→GAA) osmEosmotically‑inducible lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base T (30/33);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

CCTTGACCTGTGGTTATGACCCTCTTCACATTTCGGGCAAATAATTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCAC  >  NC_000913/1821788‑1822082
                                                                                                                                                    |                                                                                                                                                  
ccTTGACCTGTGGTTATGACCCTCTTCACATTTCGGGCAAATAATTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTc                                                                                                                                                  <  1:710993/151‑1 (MQ=255)
    gACCTGTGGTTATGACCCTCTTCACATTTCGGGCAAATAATTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAgtg                                                                                                                                              >  2:510293/1‑151 (MQ=255)
     aCCTGTGGTTATGACCCTCTTCACATTTCGGGCAAATAATTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAgtgt                                                                                                                                             >  2:333778/1‑151 (MQ=255)
       cTGTGGTTATGACCCTCTTCACATTTCGGGCAAATAATTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCg                                                                                                                                           >  1:112453/1‑151 (MQ=255)
           ggTTATGACCCTCTTCACATTTCGGGCAAATAATTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGATCAGTGTCGTAtt                                                                                                                                       <  2:628789/151‑1 (MQ=255)
            gTTATGACCCTCTTCACATTTCGGGCAAATAATTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTc                                                                                                                                      <  1:529126/151‑1 (MQ=255)
                  aCCCTCTTCACATTTCGGGCAAATAATTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGcac                                                                                                                                 >  1:359517/1‑150 (MQ=255)
                         tCACATTTCGGGCAAATAATTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCt                                                                                                                           <  2:359563/149‑1 (MQ=255)
                             aTTTCGGGCAAATAATTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAAc                                                                                                                     >  2:97063/1‑151 (MQ=255)
                             aTTTCGGGCAAATAATTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAAc                                                                                                                     <  1:3854/151‑1 (MQ=255)
                                   ggCAAATAATTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGt                                                                                                               >  2:123659/1‑151 (MQ=255)
                                        ataatTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTtgatg                                                                                                          >  1:435199/1‑151 (MQ=255)
                                        ataatTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTtgatg                                                                                                          <  2:505814/151‑1 (MQ=255)
                                         taatTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTtgatg                                                                                                          >  2:679357/1‑150 (MQ=255)
                                              cgcgAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGt                                                                                                    <  1:116026/151‑1 (MQ=255)
                                               gcgAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGt                                                                                                    >  1:325678/1‑150 (MQ=255)
                                                   aGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGt                                                                                               <  2:495541/151‑1 (MQ=255)
                                                     aTAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTat                                                                                             >  2:170416/1‑151 (MQ=255)
                                                     aTAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTat                                                                                             <  1:157871/151‑1 (MQ=255)
                                                      tAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTatc                                                                                            >  1:40054/1‑151 (MQ=255)
                                                       taTTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTatc                                                                                            >  1:506385/2‑150 (MQ=255)
                                                                  cAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGt                                                                                                               <  1:470103/120‑1 (MQ=255)
                                                                  cAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGt                                                                                                               >  2:470154/1‑120 (MQ=255)
                                                                  cAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGAc                                                                                <  2:321972/151‑1 (MQ=255)
                                                                       cAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGt                                                                           <  1:91687/151‑1 (MQ=255)
                                                                       cAAAAAAAACCGGCGAAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGt                                                                           <  1:321610/151‑1 (MQ=255)
                                                                            aaaaCCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGtt                                                                      >  2:746398/1‑151 (MQ=255)
                                                                             aaaCCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGttt                                                                     <  2:599610/151‑1 (MQ=255)
                                                                              aaCCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTc                                                                    <  1:402268/151‑1 (MQ=255)
                                                                                 cGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGGCGAATTGATTTCTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAATTGATGACATTACCGGTATCATCTTACGCGACAAAGTAGGTTTCttc                                                                 <  1:395652/151‑3 (MQ=255)
                                                                                    cgcAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCttt                                                              >  2:111891/1‑151 (MQ=255)
                                                                                    cgcAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCttt                                                              >  1:360753/1‑151 (MQ=255)
                                                                                    cgcAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCttt                                                              >  2:694850/1‑151 (MQ=255)
                                                                                      cAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCtt                                                               <  2:514871/148‑1 (MQ=255)
                                                                                       aaTGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTAcc                                                           >  2:586974/1‑151 (MQ=255)
                                                                                          ggCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATc                                                        >  2:446419/1‑151 (MQ=255)
                                                                                             cGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACg                                                     <  2:664598/151‑1 (MQ=255)
                                                                                               gTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGtt                                                   <  1:544079/151‑1 (MQ=255)
                                                                                               gTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGtt                                                   <  2:183453/151‑1 (MQ=255)
                                                                                                 ttCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTg                                                  <  2:290274/150‑1 (MQ=255)
                                                                                                 ttCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTg                                                  <  2:625471/150‑1 (MQ=255)
                                                                                                        gttACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGa                                          <  2:173966/151‑1 (MQ=255)
                                                                                                         ttACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGa                                          <  1:467781/150‑1 (MQ=255)
                                                                                                          tACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATg                                        >  2:319899/1‑151 (MQ=255)
                                                                                                           aCTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATg                                        <  1:562284/150‑1 (MQ=255)
                                                                                                            cTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGt                                       >  2:567825/1‑150 (MQ=255)
                                                                                                              cAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTa                                      <  2:442020/149‑1 (MQ=255)
                                                                                                              cAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAg                                     <  2:353869/150‑1 (MQ=255)
                                                                                                              cAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAg                                     <  2:754022/150‑1 (MQ=255)
                                                                                                                 gCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCt                                 <  1:510241/151‑1 (MQ=255)
                                                                                                                   tttCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTgg                               >  1:638869/1‑151 (MQ=255)
                                                                                                                    ttCAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGc                              >  1:521567/1‑151 (MQ=255)
                                                                                                                      cAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTTTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGc                              <  2:577270/149‑1 (MQ=255)
                                                                                                                      cAGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCa                             >  1:127727/1‑150 (MQ=255)
                                                                                                                      cAGACGAAATGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCTGAGTTGATGACATGTCCGGTATCATCTATCGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCATGATGTAGGTCTGGCa                             <  1:187500/150‑1 (MQ=255)
                                                                                                                       aGACGAATTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCa                             >  1:567978/1‑149 (MQ=255)
                                                                                                                            aaTTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACcg                      <  1:494127/151‑1 (MQ=255)
                                                                                                                            aaTTGATTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACcg                      >  2:146646/1‑151 (MQ=255)
                                                                                                                                 aTTACTTCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGc                 <  1:170383/151‑1 (MQ=255)
                                                                                                                                      ttCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATgg             >  2:650720/1‑150 (MQ=255)
                                                                                                                                       tCGCAGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGa            <  1:746323/150‑1 (MQ=255)
                                                                                                                                           aGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATg       >  1:127954/1‑151 (MQ=255)
                                                                                                                                           aGCCTGTGGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATg       >  2:451248/1‑151 (MQ=255)
                                                                                                                                                gtgGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTc    <  1:659637/149‑1 (MQ=255)
                                                                                                                                                gtgGTTCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCAc  >  2:415817/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCTTGACCTGTGGTTATGACCCTCTTCACATTTCGGGCAAATAATTCGCGAAGATAATTAAGAACTCAGGCCAAAAAAAACCGGCGCAATGGCCGGTTTCCGTTGTTACTCAAGCTTTCAGACGAATTGATTACTTCGCAGCCTGTGGATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCAC  >  NC_000913/1821788‑1822082

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑