Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,435,737 G→A 23.7% intergenic (‑118/+23) uspF ← / ← ompN stress‑induced protein, ATP‑binding protein/outer membrane pore protein N, non‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,435,7370GA23.7% 184.3 / 55.3 106intergenic (‑118/+23)uspF/ompNstress‑induced protein, ATP‑binding protein/outer membrane pore protein N, non‑specific
Reads supporting (aligned to +/- strand):  ref base G (32/48);  new base A (13/12);  total (46/60)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-01

Only 100 of 106 total aligned reads displayed.

TCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATC  >  NC_000913/1435639‑1435834
                                                                                                  |                                                                                                 
tCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGc                                                                                                  <  2:56962/100‑1 (MQ=255)
   gCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAACCCGCCAt                                                                                               >  1:748063/1‑100 (MQ=255)
    cATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATa                                                                                              <  1:1154092/100‑1 (MQ=255)
     aTAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGTTTTAAGGAGGAAAAATAAGAAGAAAAAGTCCTCTATaa                                                                                             >  2:1424424/1‑100 (MQ=255)
        aaGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAg                                                                                          >  1:2745434/1‑100 (MQ=255)
            gtAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCggg                                                                                      <  1:2227389/100‑1 (MQ=255)
             tAACGCCAGATGACACTTTTTGTGTAATGACGGAGATCAAATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGc                                                                                     >  2:2240560/1‑100 (MQ=255)
              aaCGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCa                                                                                    >  1:1172110/1‑100 (MQ=255)
              aaCGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCa                                                                                    >  2:908099/1‑100 (MQ=255)
               aCGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGCAGAATAAGCAGAAAAAGCCCGCCATAACAGCggg                                                                                      >  2:2574186/1‑97 (MQ=255)
                cGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCagg                                                                                  >  2:940685/1‑100 (MQ=255)
                  ccAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCaggag                                                                                <  1:1895468/100‑1 (MQ=255)
                    aGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAggagga                                                                              <  1:1852459/100‑1 (MQ=255)
                       tGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAttt                                                                           <  2:434278/100‑1 (MQ=255)
                       tGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAttt                                                                           >  2:6396/1‑100 (MQ=255)
                        gACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTa                                                                          <  2:355910/100‑1 (MQ=255)
                         acacTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTa                                                                          <  2:788657/99‑1 (MQ=255)
                         acacTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAg                                                                         <  2:493352/100‑1 (MQ=255)
                          cacTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGa                                                                        <  1:743958/100‑1 (MQ=255)
                          cacTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGa                                                                        >  2:1274062/1‑100 (MQ=255)
                           acTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGaa                                                                       <  2:416356/100‑1 (MQ=255)
                           acTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGa                                                                        <  1:2562769/99‑1 (MQ=255)
                            cTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAAc                                                                      <  2:733939/100‑1 (MQ=255)
                             tttttGTGTAATGATGGAGTTCACATTTTTAATTTAAATCAAAGGAGAAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACt                                                                     <  1:1386258/100‑1 (MQ=255)
                             tttttGTGTAATGACGGAGTTCACATTTTTACTTTAGAGCATAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAAc                                                                      <  2:2134428/99‑1 (MQ=255)
                              ttttGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                    <  2:447756/100‑1 (MQ=255)
                              ttttGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                    <  2:44447/100‑1 (MQ=255)
                              ttttGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCAGGCAGGAGGATTTAGAACTg                                                                    >  1:44447/1‑100 (MQ=255)
                                ttGTGTAATGTCGTAGTTCATATTTTTAATTTAGATTAAAGAATGAAGAATAAGCAGAAAAGGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGAt                                                                  <  1:487508/100‑1 (MQ=255)
                                 tgtgTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATa                                                                 >  2:2148823/1‑100 (MQ=255)
                                   tgtAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATaaa                                                               <  2:2002553/100‑1 (MQ=255)
                                   tgtAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATaaa                                                               >  1:1718822/1‑100 (MQ=255)
                                   tgtAATGACGGAGTTCACATTTTGAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAATGCCCGCCATGACAGCGGGCAGGAGGATTTAGAACTGaca                                                                 <  2:2344045/98‑3 (MQ=255)
                                     tAATGACGGAGTTCACATTTTTAATTTATATCAAAGGAGTAAGAATAAGCAGAAAAATCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAcc                                                             <  1:758838/100‑1 (MQ=255)
                                     tAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCACGCCATAACAGCGGGCAGGAGGATTTAGAACTGATaaa                                                               >  2:2123372/1‑98 (MQ=255)
                                         gACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAc                                                         >  2:596776/1‑100 (MQ=255)
                                         gACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAc                                                         <  2:1512522/100‑1 (MQ=255)
                                         gACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAc                                                         <  2:35051/100‑1 (MQ=255)
                                          acttaGTTTACATTTCTAATTTAGATCAAAGGAAGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCTGGCAGGAGGATTTAGTACTGATAAACCAGAcc                                                        <  2:740360/96‑1 (MQ=255)
                                          aCGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAcc                                                        <  1:30110/100‑1 (MQ=255)
                                          aCGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGTTTTAGAACTGATAAACCAGAcc                                                        <  2:1760107/100‑1 (MQ=255)
                                            gaaGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGAAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAGTTAGAACTGATAAACCAGa                                                          <  2:2012359/94‑1 (MQ=255)
                                               gTTCACATTTTTAATTTAGATTAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAg                                                   <  1:1191229/100‑1 (MQ=255)
                                                 tCACATTTTTAATTTATATCAAAGGAGAAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCg                                                 <  1:737662/100‑1 (MQ=255)
                                                 tCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCg                                                 <  2:1327650/100‑1 (MQ=255)
                                                  cacaTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGa                                                >  2:202182/1‑100 (MQ=255)
                                                  cacaTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAAACCTAAAGCGa                                                >  1:467025/1‑100 (MQ=255)
                                                   acaTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCACCATAACAGCGTGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGAc                                               <  2:2987376/100‑1 (MQ=255)
                                                   acaTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCACCATAACAGCGGGCAGAAGGATTTAGAACTGATAAACCAGACCTAAAGCGAc                                               >  1:546891/1‑100 (MQ=255)
                                                    caTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACa                                              >  1:1260906/1‑100 (MQ=255)
                                                    caTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACa                                              <  1:1821488/100‑1 (MQ=255)
                                                    caTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACa                                              <  2:2982536/100‑1 (MQ=255)
                                                    caTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACa                                              <  2:2939773/100‑1 (MQ=255)
                                                    caTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACa                                              <  2:1268420/100‑1 (MQ=255)
                                                    caTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAGCGACaa                                             >  2:2622180/1‑100 (MQ=255)
                                                     aTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACaa                                             <  2:2476320/100‑1 (MQ=255)
                                                     aTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACaa                                             >  2:1094408/1‑100 (MQ=255)
                                                      tttttAATTTAGATAAAAGTAGTAAGAATAAGCTTAAAAATCCCGCCATAACAGCTTGCAGTTGAATTTAAAACTTATAAACCAGACCTAAAGCTACAat                                            <  2:2409032/100‑1 (MQ=255)
                                                        tttAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAatat                                          <  2:1000120/100‑1 (MQ=255)
                                                          tAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatca                                        >  1:1256174/1‑100 (MQ=255)
                                                          tAATTTAGATAAAAGGAGGAAGAATAAGCAGAAAAAGACCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatca                                        <  1:2660051/100‑1 (MQ=255)
                                                            aTTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatcatc                                      <  1:2365463/100‑1 (MQ=255)
                                                              ttAGATCAAAGGAGGAAGAATAAGAAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCgg                                    <  1:101412/100‑1 (MQ=255)
                                                                aGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTa                                  <  1:404729/100‑1 (MQ=255)
                                                                aGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTa                                  <  1:2296654/100‑1 (MQ=255)
                                                                 gATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGt                                   >  2:2449441/1‑98 (MQ=255)
                                                                 gATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCACCATCACAGCGGGCTGGAGGATTTATAACTGATAAAACAAAACAAAAACGACAAAATCATCGGTa                                  >  1:205057/1‑99 (MQ=255)
                                                                 gATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTag                                 >  2:2663943/1‑100 (MQ=255)
                                                                  aTCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTaga                                <  1:234709/100‑1 (MQ=255)
                                                                  aTCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGAATTTAGAACTGATAAACCAGACCTAAAACGACAATATCATCGGTaga                                >  1:338063/1‑100 (MQ=255)
                                                                   tccacggaggaAGAATAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGCGGATTTAGAACTCATAAACCAGACCCAAAGCGAccatcccatc                                      >  2:2583396/6‑93 (MQ=255)
                                                                   tCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTagag                               >  2:2765592/1‑100 (MQ=255)
                                                                   tCAAAGGAGGAAGAATAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTagag                               <  1:1728032/100‑1 (MQ=255)
                                                                     aaaGGAGGAAGAATAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGAt                             >  1:2641424/1‑100 (MQ=255)
                                                                     aaaGGAGGAAGAATAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGAt                             >  1:2641431/1‑100 (MQ=255)
                                                                         ggggaAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAACTAAAGCGACAATATCATCGGTAGAGATGCCa                         >  1:2044597/3‑100 (MQ=255)
                                                                            gaagaaTAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTg                      <  2:528738/100‑1 (MQ=255)
                                                                            gaagaaTAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTg                      <  2:215548/100‑1 (MQ=255)
                                                                             aagaaTAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTgg                     <  1:3004619/100‑1 (MQ=255)
                                                                             aagaaTAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCAtt                       >  1:1519040/1‑98 (MQ=255)
                                                                              agaaTAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGc                    <  2:2586114/100‑1 (MQ=255)
                                                                               gaaTAAGTAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGc                    <  1:1984294/99‑1 (MQ=255)
                                                                                aaTAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAg                  <  2:938570/100‑1 (MQ=255)
                                                                                  tAAGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCg                >  1:1345307/1‑100 (MQ=255)
                                                                                   aaGCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGt               <  1:269376/100‑1 (MQ=255)
                                                                                    aGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTa              >  1:347638/1‑100 (MQ=255)
                                                                                    aGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTa              >  1:2413058/1‑100 (MQ=255)
                                                                                     gCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAg             <  2:2354294/100‑1 (MQ=255)
                                                                                     gCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAg             >  1:1357369/1‑100 (MQ=255)
                                                                                     gCAGAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAg             >  1:2019702/1‑100 (MQ=255)
                                                                                       aGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGa            <  1:1359600/99‑1 (MQ=255)
                                                                                        gAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAg          >  2:1961036/1‑100 (MQ=255)
                                                                                        gAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATATACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAg          <  1:848212/100‑1 (MQ=255)
                                                                                        gAAAAAGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCAATTGGCAGCGTAGAAg          >  1:1203516/1‑100 (MQ=255)
                                                                                           aaaGCCCGTGATATCAGCGGGCAGGAGGATTTAGAACTGATAAACGAGAACTAAAGCGACAATATTGTCGGTAGAAATGCCATTGGCAGCGTAGaa           >  2:1340902/1‑96 (MQ=255)
                                                                                            aaGCCCACCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGt      <  1:2232546/100‑1 (MQ=255)
                                                                                             aGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTc     >  2:2435083/1‑100 (MQ=255)
                                                                                              gCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCa    >  2:158058/1‑100 (MQ=255)
                                                                                              gCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCa    <  1:286945/100‑1 (MQ=255)
                                                                                                ccGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATc  >  1:3015897/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
TCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATC  >  NC_000913/1435639‑1435834

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