Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,901,422 T→C K267E (AAA→GAA)  bglH ← carbohydrate‑specific outer membrane porin, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,901,4220TC99.7% 1353.9 / ‑4.5 339K267E (AAA→GAA) bglHcarbohydrate‑specific outer membrane porin, cryptic
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (174/164);  total (174/165)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

Only 100 of 346 total aligned reads displayed.

TTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCA  >  NC_000913/3901322‑3901519
                                                                                                    |                                                                                                 
ttATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAAtttt                                                                                                   >  2:6649114/1‑101 (MQ=255)
  ataGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCAc                                                                                                 <  1:6648412/101‑1 (MQ=255)
    aGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCAccc                                                                                               <  1:6648464/101‑1 (MQ=255)
      cGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCgg                                                                                             <  1:6648584/101‑1 (MQ=255)
      cGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCgg                                                                                             >  2:6649124/1‑101 (MQ=255)
        tAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGAc                                                                                           <  1:6648658/101‑1 (MQ=255)
         aaaTCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGAcc                                                                                          >  1:6649128/1‑101 (MQ=255)
         aaaTCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGAcc                                                                                          >  2:6649127/1‑101 (MQ=255)
           aTCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCg                                                                                         >  2:6649133/1‑100 (MQ=255)
           aTCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGa                                                                                        <  2:6648591/101‑1 (MQ=255)
              aaTGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACtt                                                                                     >  2:6649137/1‑101 (MQ=255)
               aTGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACttt                                                                                    >  1:6649139/1‑101 (MQ=255)
                tggtattggtatTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTc                                                                                   <  2:6648575/101‑1 (MQ=255)
                 ggtattggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTcc                                                                                  >  2:6649144/1‑101 (MQ=255)
                   tattggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCATCCGGACCGACTTTCCAg                                                                                <  1:6648532/101‑1 (MQ=255)
                   tattggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAg                                                                                <  2:6648656/101‑1 (MQ=255)
                    attggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGt                                                                               <  1:6648695/101‑1 (MQ=255)
                      tggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGttt                                                                             <  2:6648702/101‑1 (MQ=255)
                      tggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGttt                                                                             <  2:6648629/101‑1 (MQ=255)
                       ggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGtttt                                                                            >  1:6649149/1‑101 (MQ=255)
                         tattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTcc                                                                          >  1:6649153/1‑101 (MQ=255)
                         tattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTcc                                                                          >  2:6649152/1‑101 (MQ=255)
                           ttAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAg                                                                        >  2:6649156/1‑101 (MQ=255)
                            tAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGa                                                                       >  2:6649162/1‑101 (MQ=255)
                            tAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGa                                                                       >  2:6649165/1‑101 (MQ=255)
                              aTCTGCTGTTTGTTTTGCAGGCTGCGATCGTAGTCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGAcc                                                                     >  2:6649170/1‑101 (MQ=255)
                               tCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCt                                                                    <  1:6648634/101‑1 (MQ=255)
                               tCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCt                                                                    >  1:6649172/1‑101 (MQ=255)
                               tCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCt                                                                    <  2:6648426/101‑1 (MQ=255)
                                ctgctgTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTa                                                                   >  2:6649174/1‑101 (MQ=255)
                                 tgctgTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTac                                                                  <  2:6648451/101‑1 (MQ=255)
                                   ctgTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTacac                                                                >  1:6649179/1‑101 (MQ=255)
                                    tgtttgttTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACAcc                                                               >  2:6649181/1‑101 (MQ=255)
                                     gtttgtttTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACAccc                                                              <  1:6648863/101‑1 (MQ=255)
                                      tttgtttTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCg                                                             >  2:6649185/1‑101 (MQ=255)
                                        tgtttTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCt                                                           <  1:6648797/101‑1 (MQ=255)
                                         gtttTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTg                                                          <  2:6648384/101‑1 (MQ=255)
                                          ttttGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGc                                                         >  2:6649191/1‑101 (MQ=255)
                                          ttttGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGc                                                         <  1:6648757/101‑1 (MQ=255)
                                           tttGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCg                                                        >  1:6649192/1‑101 (MQ=255)
                                           tttGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCg                                                        <  2:6648998/101‑1 (MQ=255)
                                            ttGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCgg                                                       >  1:6649197/1‑101 (MQ=255)
                                             tGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGc                                                      >  2:6649198/1‑101 (MQ=255)
                                              gCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCa                                                     <  2:6648665/101‑1 (MQ=255)
                                               cAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCg                                                                                              >  1:6649204/1‑59 (MQ=255)
                                               cAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCAt                                                    <  1:6648827/101‑1 (MQ=255)
                                                aGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATc                                                   <  2:6648510/101‑1 (MQ=255)
                                                aGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATc                                                   >  2:6649205/1‑101 (MQ=255)
                                                aGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATc                                                   >  2:6649209/1‑101 (MQ=255)
                                                 ggCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCa                                                  <  1:6648602/101‑1 (MQ=255)
                                                  gCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAg                                                 <  2:6648719/101‑1 (MQ=255)
                                                   cTGCGATCATAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGt                                                >  1:6649214/1‑101 (MQ=255)
                                                     gCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTAc                                              <  2:6648762/101‑1 (MQ=255)
                                                     gCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTAc                                              >  2:6649216/1‑101 (MQ=255)
                                                       gATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGc                                            >  1:6649220/1‑101 (MQ=255)
                                                        aTCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCt                                           >  2:6649222/1‑101 (MQ=255)
                                                        aTCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCt                                           <  1:6648673/101‑1 (MQ=255)
                                                         tCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTg                                          >  2:6649225/1‑101 (MQ=255)
                                                          cGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGc                                         <  2:6648635/101‑1 (MQ=255)
                                                          cGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGc                                         >  1:6649228/1‑101 (MQ=255)
                                                           gTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGc                                                         <  1:6648692/84‑1 (MQ=255)
                                                           gTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCg                                        <  2:6648465/101‑1 (MQ=255)
                                                              atcatcGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGttt                                     >  2:6649235/1‑101 (MQ=255)
                                                              atcatcGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGttt                                     >  1:6649234/1‑101 (MQ=255)
                                                                catcGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTc                                   >  2:6649239/1‑101 (MQ=255)
                                                                  tcGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCa                                 <  2:6648184/101‑1 (MQ=255)
                                                                  tcGACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCa                                 >  1:6649242/1‑101 (MQ=255)
                                                                    gACATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGt                               <  1:6648352/101‑1 (MQ=255)
                                                                      cATCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCa                             >  1:6649245/1‑101 (MQ=255)
                                                                        tCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAg                           >  2:6649249/1‑101 (MQ=255)
                                                                        tCTTCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAg                           <  2:6648497/101‑1 (MQ=255)
                                                                           tCGCGAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCAt                        >  2:6649255/1‑101 (MQ=255)
                                                                            cgcgAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATc                       <  2:6648778/101‑1 (MQ=255)
                                                                             gcgAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCt                      <  2:6648724/101‑1 (MQ=255)
                                                                             gcgAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCt                      >  1:6649259/1‑101 (MQ=255)
                                                                              cgAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTg                     <  2:6648639/101‑1 (MQ=255)
                                                                              cgAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTg                     >  1:6649262/1‑101 (MQ=255)
                                                                              cgAACCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTg                     <  2:6648793/101‑1 (MQ=255)
                                                                                 aCCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGAt                  >  2:6649266/1‑101 (MQ=255)
                                                                                 aCCAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGAt                  <  2:6648559/101‑1 (MQ=255)
                                                                                  ccAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGAtt                 >  1:6649268/1‑101 (MQ=255)
                                                                                  ccAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGAtt                 <  1:6649020/101‑1 (MQ=255)
                                                                                   cAGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGAttt                >  2:6649271/1‑101 (MQ=255)
                                                                                    aGCGCGATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTc               <  2:6648533/101‑1 (MQ=255)
                                                                                      cgcgATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGc                                                         >  1:6649276/1‑57 (MQ=255)
                                                                                       gcgATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAAt            >  2:6649277/1‑101 (MQ=255)
                                                                                         gATATCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAAttt          <  2:6648803/101‑1 (MQ=255)
                                                                                          atatCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAAtttt         <  1:6648838/101‑1 (MQ=255)
                                                                                          atatCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAAtttt         <  2:6648944/101‑1 (MQ=255)
                                                                                           tatCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTc        <  2:6648828/101‑1 (MQ=255)
                                                                                           tatCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTc        >  2:6649283/1‑101 (MQ=255)
                                                                                            atCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAAtttt         >  1:6649287/1‑99 (MQ=255)
                                                                                            atCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCg       <  2:6648788/101‑1 (MQ=255)
                                                                                             tCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCTAGTCAAGCATCTGGATTTCAATTTTCgg      >  2:6649288/1‑101 (MQ=255)
                                                                                             tCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCgg      <  2:6648923/101‑1 (MQ=255)
                                                                                             tCAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCgg      <  1:6648744/101‑1 (MQ=255)
                                                                                              cAATTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGgc     <  2:6648708/101‑1 (MQ=255)
                                                                                               aaTTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAAtttt         >  2:6649293/1‑96 (MQ=255)
                                                                                                aTTTCACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGgcgc   >  1:6649294/1‑101 (MQ=255)
                                                                                                 tttcACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCa  <  2:6648910/101‑1 (MQ=255)
                                                                                                    |                                                                                                 
TTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCA  >  NC_000913/3901322‑3901519

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