Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,462,003 G→A G233D (GGC→GAC)  fadL → long‑chain fatty acid outer membrane channel/bacteriophage T2 receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,462,0030GA99.7% 1461.1 / NA 342G233D (GGC→GAC) fadLlong‑chain fatty acid outer membrane channel/bacteriophage T2 receptor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (166/175);  minor base C (0/1);  total (166/176)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 342 total aligned reads displayed.

TGGCCAAACTCAGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGGCTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCAAAGGTAACTACAGC  >  NC_000913/2461905‑2462100
                                                                                                  |                                                                                                 
 ggCCAAACTCAGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTg                                                                                                >  1:4125891/1‑101 (MQ=255)
 ggCCAAACTCAGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTg                                                                                                <  2:4125642/101‑1 (MQ=255)
  gCCAAACTCAGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTgg                                                                                               >  2:4125892/1‑101 (MQ=255)
  gCCAAACTCAGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTgg                                                                                               <  2:4125788/101‑1 (MQ=255)
    cAAACTCAGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGaa                                                                                             <  1:4125323/101‑1 (MQ=255)
     aaaCTCAGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAAc                                                                                            >  2:4125896/1‑101 (MQ=255)
     aaaCTCAGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAAc                                                                                            <  1:4125595/101‑1 (MQ=255)
     aaaCTCAGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAAc                                                                                            <  2:4125479/101‑1 (MQ=255)
      aaCTCAGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACg                                                                                           >  1:4125900/1‑101 (MQ=255)
        cTCAGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGcc                                                                                         >  2:4125905/1‑101 (MQ=255)
          cAGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCgg                                                                                       <  2:4125455/101‑1 (MQ=255)
          cAGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCgg                                                                                       <  2:4125503/101‑1 (MQ=255)
          cAGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCgg                                                                                       >  2:4125908/1‑101 (MQ=255)
           aGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGa                                                                                      >  1:4125909/1‑101 (MQ=255)
             cAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAAt                                                                                    >  1:4125912/1‑101 (MQ=255)
              aaGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATc                                                                                   <  2:4125423/101‑1 (MQ=255)
               aGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATcc                                                                                  <  2:4125409/101‑1 (MQ=255)
                 ggCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTg                                                                                >  2:4125919/1‑101 (MQ=255)
                  gCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGt                                                                               <  2:4125387/101‑1 (MQ=255)
                  gCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGt                                                                               >  1:4125920/1‑101 (MQ=255)
                    aaGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTAt                                                                             <  1:4125129/101‑1 (MQ=255)
                    aaGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTAt                                                                             <  1:4125621/101‑1 (MQ=255)
                    aaGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTAt                                                                             >  1:4125925/1‑101 (MQ=255)
                     aGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATg                                                                            <  1:4125433/101‑1 (MQ=255)
                      gCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGa                                                                           <  2:4125317/101‑1 (MQ=255)
                      gCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGa                                                                           <  2:4125441/101‑1 (MQ=255)
                       cATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGaa                                                                          >  2:4125927/1‑101 (MQ=255)
                       cATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGaa                                                                          <  1:4125677/101‑1 (MQ=255)
                          tGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTg                                                                       <  1:4125451/101‑1 (MQ=255)
                          tGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTg                                                                       >  2:4125932/1‑101 (MQ=255)
                          tGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTg                                                                       <  1:4125360/101‑1 (MQ=255)
                            gCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATg                                                                            >  1:4125933/1‑94 (MQ=255)
                            gCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGa                                                                     <  1:4125326/101‑1 (MQ=255)
                              aGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATa                                                                   <  1:4125474/101‑1 (MQ=255)
                              aGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATa                                                                   <  2:4125153/101‑1 (MQ=255)
                               gCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATaa                                                                  >  2:4125939/1‑101 (MQ=255)
                                 tACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATaaaa                                                                >  1:4125942/1‑101 (MQ=255)
                                 tACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATaaaa                                                                <  2:4125567/101‑1 (MQ=255)
                                   ccGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAAt                                                              >  2:4125943/1‑101 (MQ=255)
                                      ccAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAAc                                                           >  1:4125947/1‑101 (MQ=255)
                                       cAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATa                                                             <  1:4125402/98‑1 (MQ=255)
                                       cAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAAcc                                                          >  2:4125952/1‑101 (MQ=255)
                                        aaCGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATaaaa                                                                <  1:4125186/94‑1 (MQ=255)
                                        aaCGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCg                                                         >  1:4125955/1‑101 (MQ=255)
                                          cGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCt                                                       >  2:4125959/1‑101 (MQ=255)
                                            gTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTAt                                                     <  1:4125637/101‑1 (MQ=255)
                                            gTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTAt                                                     >  1:4125963/1‑101 (MQ=255)
                                             tATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATg                                                    <  2:4125607/101‑1 (MQ=255)
                                             tATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATg                                                    >  2:4125964/1‑101 (MQ=255)
                                             tATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATg                                                    <  2:4125278/101‑1 (MQ=255)
                                               ttGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCa                                                  >  1:4125970/1‑101 (MQ=255)
                                                tGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCAc                                                 >  1:4125973/1‑101 (MQ=255)
                                                 gACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACt                                                <  1:4125512/101‑1 (MQ=255)
                                                 gACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACt                                                <  2:4125625/101‑1 (MQ=255)
                                                 gACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACt                                                >  2:4125975/1‑101 (MQ=255)
                                                  aCAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTg                                               >  2:4125979/1‑101 (MQ=255)
                                                      tAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCt                                           >  1:4125982/1‑101 (MQ=255)
                                                      tAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCt                                           <  2:4125532/101‑1 (MQ=255)
                                                       aaTACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTa                                          >  1:4125988/1‑101 (MQ=255)
                                                       aaTACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTa                                          <  1:4125609/101‑1 (MQ=255)
                                                       aaTACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTa                                          <  2:4125486/101‑1 (MQ=255)
                                                        aTACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTAc                                         <  2:4125657/101‑1 (MQ=255)
                                                        aTACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTAc                                         <  1:4125242/101‑1 (MQ=255)
                                                         tACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTAcc                                        >  2:4125990/1‑101 (MQ=255)
                                                          aCCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCg                                       >  1:4125992/1‑101 (MQ=255)
                                                           ccAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCg                                       <  2:4125525/100‑1 (MQ=255)
                                                            cAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGtt                                     <  1:4125497/101‑1 (MQ=255)
                                                               aaTCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTg                                  <  2:4125575/101‑1 (MQ=255)
                                                               aaTCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTg                                  >  2:4126001/1‑101 (MQ=255)
                                                               aaTCGCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTg                                  >  1:4125999/1‑101 (MQ=255)
                                                                   gCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGt                              <  2:4125557/101‑1 (MQ=255)
                                                                   gCTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGt                              >  1:4126005/1‑101 (MQ=255)
                                                                    cTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTg                             <  1:4125783/101‑1 (MQ=255)
                                                                    cTCATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTg                             >  2:4126007/1‑101 (MQ=255)
                                                                      cATCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGaa                           <  2:4125695/101‑1 (MQ=255)
                                                                       aTCTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGaaa                          <  1:4125550/101‑1 (MQ=255)
                                                                         cTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAAt                        <  1:4125377/101‑1 (MQ=255)
                                                                         cTGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAAt                        >  2:4126012/1‑101 (MQ=255)
                                                                          tGAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAAtt                       <  1:4125535/101‑1 (MQ=255)
                                                                           gAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTg                      >  2:4126017/1‑101 (MQ=255)
                                                                           gAACGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTg                      >  1:4126015/1‑101 (MQ=255)
                                                                             aCGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGAc                    >  2:4126022/1‑101 (MQ=255)
                                                                              cGGTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACt                   >  1:4126023/1‑101 (MQ=255)
                                                                               ggTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACtt                  >  1:4126027/1‑101 (MQ=255)
                                                                                gTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTc                 <  2:4125496/101‑1 (MQ=255)
                                                                                gTAACCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTc                 <  1:4125735/101‑1 (MQ=255)
                                                                                  aaCCAGTNGGGCTTTGACTGGAACGCCGGAATNCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCaa               >  1:4126033/1‑101 (MQ=255)
                                                                                  aaCCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCaa               >  2:4126032/1‑101 (MQ=255)
                                                                                   aCCAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCaaa              <  1:4125343/101‑1 (MQ=255)
                                                                                    ccAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCAAAg             <  1:4125585/101‑1 (MQ=255)
                                                                                    ccAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCAAAg             <  1:4125460/101‑1 (MQ=255)
                                                                                     cAGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCAAAgg            >  1:4126039/1‑101 (MQ=255)
                                                                                      aGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCAAAGGt           >  2:4126041/1‑101 (MQ=255)
                                                                                      aGTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCAAAGGt           >  1:4126044/1‑101 (MQ=255)
                                                                                       gTGGGGCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCAAAGGTa          <  2:4125742/101‑1 (MQ=255)
                                                                                          gggCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCAAAGGTAACt       >  2:4126048/1‑101 (MQ=255)
                                                                                           ggCTTTGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCAAAGGTAACTa      >  1:4126050/1‑101 (MQ=255)
                                                                                              tttGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCAAAGGTAACTAcag   >  1:4126056/1‑101 (MQ=255)
                                                                                               ttGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCAAAGGTAACTAcagc  <  2:4125469/101‑1 (MQ=255)
                                                                                               ttGACTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCAAAGGTAACTAcagc  <  2:4125356/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
TGGCCAAACTCAGCAAGGGCAAGCATTGGCAGCTACCGCCAACGGTATTGACAGTAATACCAAAATCGCTCATCTGAACGGTAACCAGTGGGGCTTTGGCTGGAACGCCGGAATCCTGTATGAACTGGATAAAAATAACCGCTATGCACTGACCTACCGTTCTGAAGTGAAAATTGACTTCAAAGGTAACTACAGC  >  NC_000913/2461905‑2462100

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