Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,958,2720TG22.7% 503.2 / 141.3 238intergenic (‑140/+248)torY/cutCTMAO reductase III (TorYZ), cytochrome c‑type subunit/copper homeostasis protein
Reads supporting (aligned to +/- strand):  ref base T (97/87);  new base G (26/28);  total (123/115)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 239 total aligned reads displayed.

TAATTAAATATACAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAATATATACAATATTTTTAGCA  >  NC_000913/1958174‑1958370
                                                                                                  |                                                                                                  
taATTAAATATACAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAtt                                                                                                   <  2:4433540/100‑1 (MQ=255)
  aTTAAATATACAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGc                                                                                                 <  2:1650897/100‑1 (MQ=255)
  aTTAAATATACAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGc                                                                                                 <  1:1404949/100‑1 (MQ=255)
  aTTAAATATACAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGc                                                                                                 >  2:310038/1‑100 (MQ=255)
   ttAAATATACAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCt                                                                                                >  1:3719824/1‑100 (MQ=255)
     aaaTATACAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTaa                                                                                              >  1:5133058/1‑100 (MQ=255)
      aaTATACAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAg                                                                                             <  2:3788872/100‑1 (MQ=255)
        tataCAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACATTCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGa                                                                                           >  2:4972782/1‑100 (MQ=255)
         ataCAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGaa                                                                                          <  1:2906986/100‑1 (MQ=255)
          taCAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAAt                                                                                         >  1:4350166/1‑100 (MQ=255)
          taCAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAAt                                                                                         >  2:5264911/1‑100 (MQ=255)
           aCAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATa                                                                                        >  1:3752367/1‑100 (MQ=255)
           aCAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATa                                                                                        <  2:4114969/100‑1 (MQ=255)
            cAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATaa                                                                                       >  1:1726654/1‑100 (MQ=255)
            cAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATaa                                                                                       >  2:969126/1‑100 (MQ=255)
             aaCGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAAc                                                                                      >  1:2124220/1‑100 (MQ=255)
               ctaaTTAACCGCCTGCATAATTAGGCACAACACTGCCTGATACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAAccc                                                                                    <  2:1510006/98‑1 (MQ=255)
               cGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAAccc                                                                                    <  2:3539675/100‑1 (MQ=255)
                 aaTTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAg                                                                                  <  1:3967949/100‑1 (MQ=255)
                 aaTTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAg                                                                                  <  1:1320371/100‑1 (MQ=255)
                   ttAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGtt                                                                                <  1:5220096/100‑1 (MQ=255)
                    tAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTg                                                                               <  2:1678603/100‑1 (MQ=255)
                    tAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTg                                                                               <  1:3241766/100‑1 (MQ=255)
                     aaCCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGCTAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGc                                                                              >  1:4740387/1‑100 (MQ=255)
                     aaCCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGc                                                                              <  1:16649/100‑1 (MQ=255)
                      aCCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCt                                                                             >  2:1884955/1‑100 (MQ=255)
                      aCCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCt                                                                             >  2:5099967/1‑100 (MQ=255)
                           ttGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTa                                                                        <  1:1260819/100‑1 (MQ=255)
                           ttGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTa                                                                        <  1:1900884/100‑1 (MQ=255)
                            tGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAg                                                                       <  1:1175572/100‑1 (MQ=255)
                            tGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAg                                                                       >  2:858271/1‑100 (MQ=255)
                            tGCATAATTAGGCACAACAATCCCTGAAACAATCGATAAAGAATATGATTTATTACAAGGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAg                                                                       >  1:4754947/1‑100 (MQ=255)
                              cATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGaa                                                                     >  2:2099685/1‑100 (MQ=255)
                               aTAATTAGGCACGACACTGCCTGAAACAATGGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAat                                                                    <  2:4518954/100‑1 (MQ=255)
                               aTAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAat                                                                    >  2:1308775/1‑100 (MQ=255)
                               aTAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAat                                                                    <  2:4013666/100‑1 (MQ=255)
                               aTAATTAGGCACAACACTGCCTGAAACAAGCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAat                                                                    <  2:1799879/100‑1 (MQ=255)
                                 aaTTAGGCACAACACCGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGTatat                                                                  <  2:4655289/100‑1 (MQ=255)
                                    tAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAatatata                                                               <  1:5042743/100‑1 (MQ=255)
                                    tAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAatatata                                                               <  1:1821614/100‑1 (MQ=255)
                                    tAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAatatata                                                               <  1:3771892/100‑1 (MQ=255)
                                      ggCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATAtt                                                             >  2:2289803/1‑100 (MQ=255)
                                        cacaACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTa                                                           <  2:2406178/100‑1 (MQ=255)
                                         acaacgCTGGCTGGAGCAATCGGTAGAGAATGTTGGTTATGACAATGTATTCATTAATTGCGAAGGAATAACCCAGTTGCTATTTAGAttatttatttat                                                          >  2:1099807/1‑100 (MQ=255)
                                         acaacaCTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTAt                                                          >  2:1398337/1‑100 (MQ=255)
                                          caacaCTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATc                                                         >  2:2018230/1‑100 (MQ=255)
                                          caacaCTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATc                                                         <  2:3881960/100‑1 (MQ=255)
                                              acTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGtt                                                     <  1:1479002/100‑1 (MQ=255)
                                               cTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAACGAATAGCCCGGTTGCTATTTAGAATATATATTTATCAGttt                                                    <  2:3484644/100‑1 (MQ=255)
                                                tGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGtttt                                                   <  1:56502/100‑1 (MQ=255)
                                                 gCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTa                                                  >  2:4574652/1‑100 (MQ=255)
                                                 gCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTa                                                  >  2:2375930/1‑100 (MQ=255)
                                                 gCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTa                                                  >  1:1333750/1‑100 (MQ=255)
                                                 gCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTa                                                  >  2:3797862/1‑100 (MQ=255)
                                                 gCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTa                                                  >  1:1014307/1‑100 (MQ=255)
                                                 gCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTa                                                  >  2:1792214/1‑100 (MQ=255)
                                                  ccTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAg                                                 >  1:4975774/1‑100 (MQ=255)
                                                  ccTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAg                                                 <  2:145838/100‑1 (MQ=255)
                                                  ccTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAg                                                 <  1:3003768/100‑1 (MQ=255)
                                                    tGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTa                                               >  1:3402517/1‑100 (MQ=255)
                                                    tGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTa                                               >  2:2943776/1‑100 (MQ=255)
                                                      aaaCAATCGATAAAGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAAt                                             <  1:4101583/100‑1 (MQ=255)
                                                         cAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAAtttn                                          <  1:4436222/100‑2 (MQ=255)
                                                          aaTCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTaa                                         >  2:3604386/1‑100 (MQ=255)
                                                            tCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAAt                                       <  2:5317535/100‑1 (MQ=255)
                                                              gATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATcc                                     >  2:4344795/1‑100 (MQ=255)
                                                              gATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATcc                                     >  1:2650864/1‑100 (MQ=255)
                                                              gATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATcc                                     <  2:589451/100‑1 (MQ=255)
                                                              gATAAAGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATcc                                     <  2:1000872/100‑1 (MQ=255)
                                                                 aaaGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCAt                                  >  1:4845603/1‑100 (MQ=255)
                                                                  aaGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGCATAACCCAGTTGCTATTTAGAAAATATATTTATCAGTTTTAGTAATTTAAATCCCAta                                 >  1:3337877/1‑100 (MQ=255)
                                                                  aaGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCAta                                 >  1:1526455/1‑100 (MQ=255)
                                                                  aaGAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCAta                                 <  1:2811233/100‑1 (MQ=255)
                                                                    gAATATGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCAtaat                               <  2:82312/100‑1 (MQ=255)
                                                                       tatGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCAtaattaa                            <  2:1474359/100‑1 (MQ=255)
                                                                       tatGATTTATTACAATGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCAtaattaa                            <  1:5498527/100‑1 (MQ=255)
                                                                        atGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCAtaattaat                           >  2:3394322/1‑100 (MQ=255)
                                                                         tGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATNTATANTTATCAGNTTTAGTAATTTAAATCCCATAATTAAtg                          <  1:3527500/100‑1 (MQ=255)
                                                                          gATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAAtgt                         >  2:4842018/1‑100 (MQ=255)
                                                                          gATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAAtgt                         <  1:2313349/100‑1 (MQ=255)
                                                                          gATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAAtgt                         >  2:3240337/1‑100 (MQ=255)
                                                                          gATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAAtgt                         <  1:2563694/100‑1 (MQ=255)
                                                                             ttattaCAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAGATCCCATAATTAATGTGaa                      <  2:5033646/100‑1 (MQ=255)
                                                                               attaCAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTACTTTAAATCCCATAATTAATGTGAata                    >  1:2068396/1‑100 (MQ=255)
                                                                                 taCAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAatata                  >  2:4078742/1‑100 (MQ=255)
                                                                                 taCAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAatata                  <  2:5465449/100‑1 (MQ=255)
                                                                                  aCAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAatatat                 >  2:1763950/1‑100 (MQ=255)
                                                                                   cAATGTAATTATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAAtatata                <  2:2113868/100‑1 (MQ=255)
                                                                                   cAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAAtatata                <  1:2931831/100‑1 (MQ=255)
                                                                                      tGTAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAATATATACaa             >  2:427505/1‑100 (MQ=255)
                                                                                        tAATCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAATATATACAata           <  1:669772/100‑1 (MQ=255)
                                                                                         aaTCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAATATATACAatat          <  2:713579/100‑1 (MQ=255)
                                                                                           tCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAATATATACAATAttt        >  2:271075/1‑100 (MQ=255)
                                                                                           tCATTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAATATATACAATAttt        <  2:2602394/100‑1 (MQ=255)
                                                                                             aTTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAATATATACAATAttttt      >  1:5392596/1‑100 (MQ=255)
                                                                                             aTTAAGTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAATATATACAATAttttt      <  1:2172276/100‑1 (MQ=255)
                                                                                              ttAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAATATATACAATATTTTTa     <  2:1066287/100‑1 (MQ=255)
                                                                                              ttAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAATATATACAATATTTTTa     <  2:2875808/100‑1 (MQ=255)
                                                                                                aaTTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAATATATACAATATTTTTAGc   >  2:2158129/1‑100 (MQ=255)
                                                                                                 aTTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAATATATACAATATTTTTAGca  <  2:5401286/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
TAATTAAATATACAACGAATTAACCGCTTGCATAATTAGGCACAACACTGCCTGAAACAATCGATAAAGAATATGATTTATTACAATGTAATCATTAATTGCTAAGGAATAACCCAGTTGCTATTTAGAATATATATTTATCAGTTTTAGTAATTTAAATCCCATAATTAATGTGAATATATACAATATTTTTAGCA  >  NC_000913/1958174‑1958370

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