Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,188,043 (C)6→5 coding (372/1035 nt) yehA ← putative fimbrial adhesin YehA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,188,0380C.98.5% 269.1 / ‑2.5 65coding (377/1035 nt)yehAputative fimbrial adhesin YehA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (32/32);  minor base G (0/1);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAA  >  NC_000913/2187899‑2188166
                                                                                                                                           |                                                                                                                                
aTTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCCCTCCCGCGTGTCGATGATGAAGCCCTGGATAATGTCGCAACGCCATTAATCATCACTGAACcgccgc                                                                                                                                  <  2:245065/140‑1 (MQ=255)
aTTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACcgccgc                                                                                                                                  <  1:745299/140‑1 (MQ=255)
aTTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACcgccgc                                                                                                                                  <  1:255867/140‑1 (MQ=255)
aTTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACcgccgc                                                                                                                                  <  1:309654/140‑1 (MQ=255)
    cTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGccccc                                                                                                                              >  1:190310/1‑140 (MQ=255)
       tAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGcccccgca                                                                                                                           >  1:34474/1‑137 (MQ=255)
        aGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGcccccgcag                                                                                                                          <  2:777649/140‑5 (MQ=255)
         gCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGa                                                                                                                        >  1:706824/1‑140 (MQ=255)
          cTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGaa                                                                                                                       <  1:269764/140‑1 (MQ=255)
                    aTGCATCCCATGACTTGCCACCTAATAAACATAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACtt                                                                                                             >  2:301556/1‑140 (MQ=255)
                    aTGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACtt                                                                                                             >  2:685920/1‑140 (MQ=255)
                       cATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCAc                                                                                                               <  2:719608/135‑1 (MQ=255)
                       cATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCAc                                                                                                               >  1:719608/1‑135 (MQ=255)
                       cATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTg                                                                                                          >  1:7244/1‑140 (MQ=255)
                       cATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTg                                                                                                          >  2:446337/1‑140 (MQ=255)
                             aTGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTAc                                                                                                    <  2:7244/140‑1 (MQ=255)
                                aCTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTgg                                                                                                 <  2:133194/140‑1 (MQ=255)
                                    gCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGc                                                                                             >  2:27714/1‑140 (MQ=255)
                                    gCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCGGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGc                                                                                             <  1:425360/140‑1 (MQ=255)
                                        ccTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGgt                                                                                         <  1:125605/140‑1 (MQ=255)
                                        ccTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGgt                                                                                         <  1:286920/140‑1 (MQ=255)
                                          taataaACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGgtgt                                                                                       <  2:32270/140‑1 (MQ=255)
                                          taataaACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGgtgt                                                                                       <  1:455396/140‑1 (MQ=255)
                                           aataaACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTc                                                                                      >  1:764716/1‑140 (MQ=255)
                                                acacAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGt                                                                                 <  2:706824/140‑1 (MQ=255)
                                                 cacaGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTa                                                                                <  2:34474/140‑1 (MQ=255)
                                                 cacaGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTa                                                                                <  1:27714/140‑1 (MQ=255)
                                                 cacaGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTa                                                                                <  1:234378/140‑1 (MQ=255)
                                                   caGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAAc                                                                              <  2:190310/140‑1 (MQ=255)
                                                   caGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAAc                                                                              >  1:654649/1‑140 (MQ=255)
                                                   caGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAAc                                                                              >  2:118146/1‑140 (MQ=255)
                                                    aGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACt                                                                             >  2:361809/1‑140 (MQ=255)
                                                    aGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACt                                                                             >  2:797765/1‑140 (MQ=255)
                                                      tAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGGCGCAGGTGTCACCGTAACTGa                                                                           >  2:754718/1‑140 (MQ=255)
                                                      tAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGgtgt                                                                                       <  2:362839/128‑1 (MQ=255)
                                                      tAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGgtgt                                                                                       >  1:362839/1‑128 (MQ=255)
                                                      tAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGt                                                                                 <  1:352611/134‑1 (MQ=255)
                                                      tAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGt                                                                                 >  2:352611/1‑134 (MQ=255)
                                                      tAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGa                                                                           >  1:579580/1‑140 (MQ=255)
                                                        gAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAAt                                                                         <  1:711474/140‑1 (MQ=255)
                                                        gAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAAt                                                                         <  2:764716/140‑1 (MQ=255)
                                                        gAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACAGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAAt                                                                         <  2:740496/140‑1 (MQ=255)
                                                         aaaTTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATa                                                                        >  1:158687/1‑140 (MQ=255)
                                                            ttGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAAcc                                                                     <  1:446337/140‑1 (MQ=255)
                                                             tGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCa                                                                    <  2:654649/140‑1 (MQ=255)
                                                              gTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACAGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTGACCGTAACTGAATAACCAc                                                                   >  2:753425/1‑140 (MQ=255)
                                                                aCTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTg                                                                 >  2:780175/1‑140 (MQ=255)
                                                                         cTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTg                                                        >  1:562622/1‑140 (MQ=255)
                                                                            ccGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTc                                                     >  2:365796/1‑140 (MQ=255)
                                                                             cgcgTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTGACCGTAACTGAATAACCACTGGTGAAACTGGTCt                                                    <  1:753425/140‑1 (MQ=255)
                                                                                 tgtCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGatt                                                >  2:610041/1‑140 (MQ=255)
                                                                                 tgtCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGatt                                                >  2:62744/1‑140 (MQ=255)
                                                                                   tCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACCGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGCGAAACTGGTCTGattat                                              >  1:275220/1‑140 (MQ=255)
                                                                                      atgatgAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTg                                           <  1:623342/140‑1 (MQ=255)
                                                                                       tgatgaAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGa                                          <  1:301556/140‑1 (MQ=255)
                                                                                         atgaAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACa                                        <  1:685920/140‑1 (MQ=255)
                                                                                             aGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGt                                    >  1:425159/1‑140 (MQ=255)
                                                                                                  tGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGAtgt                               >  2:130779/1‑140 (MQ=255)
                                                                                                  tGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGAtgt                               <  1:95376/140‑1 (MQ=255)
                                                                                                         tGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGAtgtt                              >  2:638043/1‑134 (MQ=255)
                                                                                                         tGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGAtgtt                              <  1:638043/134‑1 (MQ=255)
                                                                                                         tGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGttttt                        >  2:603095/1‑140 (MQ=255)
                                                                                                          gTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTa                       <  2:158687/140‑1 (MQ=255)
                                                                                                          gTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTa                       <  1:201519/140‑1 (MQ=255)
                                                                                                          gTCGCAACACCATTAATCATCACTGAACCGCCGCCACCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTa                       <  2:579580/140‑1 (MQ=255)
                                                                                                              cAACACCATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTgg                   <  1:754718/140‑1 (MQ=255)
                                                                                                                    cATTAATCATCACTGATCCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCTCTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCAtt             <  1:118146/140‑1 (MQ=255)
                                                                                                                    cATTAATCATCACTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCAtt             >  1:600287/1‑140 (MQ=255)
                                                                                                                            atcaCTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGt     >  2:253639/1‑140 (MQ=255)
                                                                                                                            atcaCTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGt     <  1:62744/140‑1 (MQ=255)
                                                                                                                              caCTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTaa   >  2:79596/1‑140 (MQ=255)
                                                                                                                               aCTGAACCGCCGCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTaaa  >  2:678024/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                
ATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAA  >  NC_000913/2187899‑2188166

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