Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,448,947 A→C T477P (ACG→CCG)  feaB → phenylacetaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,448,9470AC96.4% 373.0 / 1.3 111T477P (ACG→CCG) feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base C (57/49);  total (60/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 111 total aligned reads displayed.

GAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACG  >  NC_000913/1448851‑1449044
                                                                                                |                                                                                                 
gAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCgg                                                                                                >  2:2611620/1‑100 (MQ=255)
  atatAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGc                                                                                              >  1:2667388/1‑100 (MQ=255)
   tataGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGcc                                                                                             <  1:309361/100‑1 (MQ=255)
     naGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGt                                                                                           >  1:1281237/2‑100 (MQ=255)
       gCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGa                                                                                         >  2:2676057/1‑100 (MQ=255)
       gCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGCTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGa                                                                                         >  2:612054/1‑100 (MQ=255)
         gATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCGGGACCGGGCCGTGAtt                                                                                       <  2:1561002/100‑1 (MQ=255)
         gATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGAtt                                                                                       >  1:426076/1‑100 (MQ=255)
         gATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGAtt                                                                                       <  2:1306039/100‑1 (MQ=255)
          aTCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGAttt                                                                                      >  2:328629/1‑100 (MQ=255)
           tCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGAtttt                                                                                     <  2:458720/100‑1 (MQ=255)
            tgCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTg                                                                                    <  2:248388/99‑1 (MQ=255)
            cGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTg                                                                                    <  1:1212090/100‑1 (MQ=255)
            cGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTg                                                                                    <  1:328625/100‑1 (MQ=255)
             gCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTgg                                                                                   <  1:3073741/100‑1 (MQ=255)
             gCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTgg                                                                                   >  2:217886/1‑100 (MQ=255)
             gCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTga                                                                                   >  1:248385/1‑99 (MQ=255)
               ttACAGGCAGGGACGGTGTGGGTAAACAGCCATACGTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGcc                                                                                 <  1:3038264/100‑1 (MQ=255)
               ttACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGcc                                                                                 <  2:3089755/100‑1 (MQ=255)
                 aCAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGcccc                                                                               <  2:445072/100‑1 (MQ=255)
                 aCAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGcccc                                                                               <  1:272413/100‑1 (MQ=255)
                  caggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCg                                                                              <  2:1664871/100‑1 (MQ=255)
                  caggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCg                                                                              <  2:1328874/100‑1 (MQ=255)
                  caggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCg                                                                              <  1:1727425/100‑1 (MQ=255)
                  caggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCg                                                                              >  2:3092440/1‑100 (MQ=255)
                  caggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCg                                                                              >  1:1269997/1‑100 (MQ=255)
                  caggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGGTTGGTGGGATGAAGCAGTCAGGATCGCGCCGTGATTTTGGCCCCg                                                                              >  1:2945080/1‑100 (MQ=255)
                    ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGAc                                                                            <  2:296631/100‑1 (MQ=255)
                    ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGAc                                                                            <  1:3439277/100‑1 (MQ=255)
                    ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGAc                                                                            >  1:2747164/1‑100 (MQ=255)
                    ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGAc                                                                            >  2:1033544/1‑100 (MQ=255)
                    ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGAc                                                                            >  2:1683553/1‑100 (MQ=255)
                    ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGAc                                                                            >  2:2098649/1‑100 (MQ=255)
                    ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGAc                                                                            >  1:3117900/1‑100 (MQ=255)
                    ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTCACGCTAACTTACCGTTTGGTGGGATGAAGCATTCAGGACCGGGCCGTGATTTTGGCCCCGAc                                                                            >  2:1475762/1‑100 (MQ=255)
                      caggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGActg                                                                          >  1:617469/1‑100 (MQ=255)
                      caggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGActg                                                                          >  2:2237194/1‑100 (MQ=255)
                       aggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGActgg                                                                         <  1:2560110/100‑1 (MQ=255)
                         ggACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGActggct                                                                       <  2:1338194/100‑1 (MQ=255)
                         ggACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGActggct                                                                       <  2:3078208/100‑1 (MQ=255)
                          gACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActggctg                                                                      <  2:1529684/100‑1 (MQ=255)
                          gACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATCAAGCAGTCAGGACCGGGCCGTGATTTTGTCCCCGActggctg                                                                      <  2:2346029/100‑1 (MQ=255)
                            cGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGa                                                                    <  1:1520028/100‑1 (MQ=255)
                            cGGTGTGGGTAAACAGCCATAACTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGa                                                                    >  1:2817744/1‑100 (MQ=255)
                             gttgtgGGTAAACAGCCTTACGTTAATTGACGCTAACGTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGActggctn                                                                      <  1:1434916/95‑2 (MQ=255)
                             ggTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGan                                                                   <  2:2941581/100‑2 (MQ=255)
                               tgtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACgg                                                                 >  1:2292129/1‑100 (MQ=255)
                               tgtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACgg                                                                 >  2:459811/1‑100 (MQ=255)
                                gtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACgg                                                                 >  1:2916294/1‑99 (MQ=255)
                                gtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACgg                                                                 >  2:2715662/1‑99 (MQ=255)
                                gtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGt                                                                >  2:3470273/1‑100 (MQ=255)
                                gtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGt                                                                >  1:2916256/1‑100 (MQ=255)
                                  gggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTg                                                              >  1:2305993/1‑100 (MQ=255)
                                     tAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgtg                                                           <  2:400220/100‑1 (MQ=255)
                                      aaaCAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgtgt                                                          <  1:1182816/100‑1 (MQ=255)
                                       aaCAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCAGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGtgtg                                                         <  2:3142854/100‑1 (MQ=255)
                                         cAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGaa                                                       <  1:3413900/100‑1 (MQ=255)
                                         cAGCCATAAATTAATTGGCGCTACCTGACCGTTTGGTGGGCTGAAGCAGGCAGGGCCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgtgt                                                          <  1:1340630/97‑1 (MQ=255)
                                             cATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTa                                                   >  2:194942/1‑100 (MQ=255)
                                               tACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAg                                                 >  1:18191/1‑100 (MQ=255)
                                                 ccTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                               <  1:2069996/100‑1 (MQ=255)
                                                  cTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCg                                              >  2:114347/1‑100 (MQ=255)
                                                      aTTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtg                                           >  1:139442/1‑99 (MQ=255)
                                                       ttGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtg                                         <  1:612047/100‑1 (MQ=255)
                                                       ttGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtg                                         >  2:1281765/1‑100 (MQ=255)
                                                        tGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtgt                                        >  1:1300419/1‑100 (MQ=255)
                                                         gACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtgtg                                       <  2:1771202/100‑1 (MQ=255)
                                                           cGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTa                                     >  2:1654622/1‑100 (MQ=255)
                                                           cGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTa                                     <  2:2673539/100‑1 (MQ=255)
                                                            gCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTAc                                    <  2:878659/100‑1 (MQ=255)
                                                                 cTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTAt                               >  2:2586921/1‑100 (MQ=255)
                                                                 cTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTAt                               >  1:2967274/1‑100 (MQ=255)
                                                                 cTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTAt                               >  1:1721553/1‑100 (MQ=255)
                                                                   tACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTa                             >  2:282711/1‑100 (MQ=255)
                                                                       gTTTGGTGGGATGAAGCAGTCAGGACCGGGTCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCt                         <  2:1721576/100‑1 (MQ=255)
                                                                              gggATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGc                  >  1:1936607/1‑100 (MQ=255)
                                                                              gggATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGc                  >  1:1762600/1‑100 (MQ=255)
                                                                              gggATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGc                  >  2:1381203/1‑100 (MQ=255)
                                                                                gATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTc                >  2:2244758/1‑100 (MQ=255)
                                                                                gATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCGGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTc                <  1:3225061/100‑1 (MQ=255)
                                                                                 aTGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCa               >  1:566050/1‑100 (MQ=255)
                                                                                  tGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCAt              >  2:1051447/1‑100 (MQ=255)
                                                                                  tGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCAt              <  1:1688129/100‑1 (MQ=255)
                                                                                   gAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGCAACTAAGTCGGTGTCTGTACGGAATTCATCAGGTTCGCtaacc             >  2:2377008/1‑96 (MQ=255)
                                                                                   gAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATa             >  2:809295/1‑100 (MQ=255)
                                                                                   gAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATa             <  2:1682844/100‑1 (MQ=255)
                                                                                   gAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATa             <  2:1683091/100‑1 (MQ=255)
                                                                                   gAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATa             <  2:1801917/100‑1 (MQ=255)
                                                                                    aaGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATaa            >  1:1567123/1‑100 (MQ=255)
                                                                                    aaGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATaa            <  2:3289125/100‑1 (MQ=255)
                                                                                     aGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAg           <  1:2499523/100‑1 (MQ=255)
                                                                                         gTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaa       <  2:3119200/100‑1 (MQ=255)
                                                                                         gTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaa       >  2:3107072/1‑100 (MQ=255)
                                                                                         gTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaa       >  2:1515125/1‑100 (MQ=255)
                                                                                         gTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaa       >  2:2402618/1‑100 (MQ=255)
                                                                                          tCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaaa      <  1:3170260/100‑1 (MQ=255)
                                                                                           cAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaaaa     <  2:1483101/100‑1 (MQ=255)
                                                                                           cAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaaaa     >  2:2732977/1‑100 (MQ=255)
                                                                                             ggACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATGAGTAAAAAAc   >  2:3056985/1‑100 (MQ=255)
                                                                                              gAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACg  >  2:415411/1‑100 (MQ=255)
                                                                                                |                                                                                                 
GAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACG  >  NC_000913/1448851‑1449044

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