Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,373,180 T→C 9.2% W97R (TGG→CGG)  arnF → undecaprenyl phosphate‑alpha‑L‑ara4N exporter; flippase ArnEF subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,373,1800TC 9.2% 390.9 / 20.5 129W97R (TGG→CGG) arnFundecaprenyl phosphate‑alpha‑L‑ara4N exporter; flippase ArnEF subunit
Reads supporting (aligned to +/- strand):  ref base T (61/55);  new base C (8/4);  total (70/59)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 129 total aligned reads displayed.

GATACTGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCGGGGCTATA  >  NC_000913/2373035‑2373325
                                                                                                                                                   |                                                                                                                                                 
gATACTGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTggg                                                                                                                                                >  4:234423/1‑149 (MQ=255)
  tACTGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGa                                                                                                                                               <  2:700132/148‑1 (MQ=255)
  tACTGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGa                                                                                                                                               <  1:15909/148‑1 (MQ=255)
  gattGTTACTTCGGGCTCCTGGGTTTTCTCTTTTCCGTTTTTGGTTGGTATAAAACGCTTATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGa                                                                                                                                               <  2:264634/145‑1 (MQ=255)
   aCTGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGaa                                                                                                                                              <  1:723240/148‑1 (MQ=255)
     tGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAgg                                                                                                                                            <  1:567901/148‑1 (MQ=255)
     tGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTAGGAAggg                                                                                                                                           >  3:778717/1‑149 (MQ=255)
     tGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCAGCCGGGAAgta                                                                                                                                           <  3:558087/149‑3 (MQ=255)
      gTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGa                                                                                                                                          >  4:598526/1‑149 (MQ=255)
      gTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGa                                                                                                                                          <  3:243990/149‑1 (MQ=255)
       ttACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAggg                                                                                                                                           >  2:651953/1‑147 (MQ=255)
       ttACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAggg                                                                                                                                           <  1:651953/147‑1 (MQ=255)
       atACTCGGGCTGCTGGGTTATCTGCTTTTCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGAc                                                                                                                                         <  4:778717/148‑1 (MQ=255)
        tACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCTTAAATTGCGGCAAGCAAAAGCCTTTCCTTTCTCAGTATGGGGTTTGTGGTGGTATGGAATTCTTCCATTGTACTTCCCGGCATGGAAggaatc                                                                                                                                        >  1:264634/1‑144 (MQ=255)
                   gctgGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTg                                                                                                                               >  1:213528/1‑148 (MQ=255)
                     tgGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGaaa                                                                                                                            <  3:51415/149‑1 (MQ=255)
                     tgGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGaaa                                                                                                                            <  4:179258/149‑1 (MQ=255)
                         ttATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCac                                                                                                                        <  3:598526/149‑1 (MQ=255)
                           aTCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCacta                                                                                                                      >  4:600472/1‑149 (MQ=255)
                            tCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCactac                                                                                                                     >  3:520561/1‑149 (MQ=255)
                               gCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTAGGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTg                                                                                                                   >  2:17992/1‑148 (MQ=255)
                                 tttCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGa                                                                                                                >  3:224987/1‑149 (MQ=255)
                                  ttCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGaa                                                                                                                             <  2:564116/135‑1 (MQ=255)
                                   tCCGTTTTTTGTTGGTATAAAACGCTTCATAATCTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCAATACTGGGAGt                                                                                                              >  4:775254/1‑149 (MQ=255)
                                   tCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAg                                                                                                               >  1:445560/1‑148 (MQ=255)
                                        tttttGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCt                                                                                                         <  3:837319/149‑1 (MQ=255)
                                          tttGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGt                                                                                                       <  3:600472/149‑1 (MQ=255)
                                           ttgttgGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGta                                                                                                      <  3:386103/149‑1 (MQ=255)
                                              ttgGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGtatta                                                                                                   >  4:330879/1‑149 (MQ=255)
                                              ttgGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGtatta                                                                                                   <  3:330879/149‑1 (MQ=255)
                                               tgGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGtatta                                                                                                   <  2:274122/148‑1 (MQ=255)
                                                  tataAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTCTGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAg                                                                                               >  3:233264/1‑149 (MQ=255)
                                                  tataAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTg                                                                                                                   >  2:363503/1‑129 (MQ=255)
                                                  tataAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTg                                                                                                                   <  1:363503/129‑1 (MQ=255)
                                                  tataAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGa                                                                                                >  2:860403/1‑148 (MQ=255)
                                                    taAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTTGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCg                                                                                             >  3:811631/1‑149 (MQ=255)
                                                    taAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCg                                                                                             >  4:188278/1‑149 (MQ=255)
                                                     aaaaCGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCgg                                                                                            >  4:315912/1‑149 (MQ=255)
                                                      aaaCGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTACTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCggg                                                                                           >  4:825069/1‑149 (MQ=255)
                                                        aCGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGtt                                                                                         <  3:188278/149‑1 (MQ=255)
                                                        aCGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGt                                                                                          >  2:890696/1‑148 (MQ=255)
                                                          gCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTc                                                                                                                                   >  3:750354/1‑105 (MQ=255)
                                                          gCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTc                                                                                                                                   <  4:750354/105‑1 (MQ=255)
                                                           cTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGAt                                                                                      <  3:234423/149‑1 (MQ=255)
                                                           cTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGAt                                                                                      <  4:84335/149‑1 (MQ=255)
                                                           cTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAAGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGa                                                                                       <  1:427587/148‑1 (MQ=255)
                                                               aTAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACtttt                                                                                                                                     >  4:500406/1‑98 (MQ=255)
                                                               aTAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTg                                                                                  <  4:346851/149‑1 (MQ=255)
                                                                tAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCt                                                                                   >  1:633571/1‑147 (MQ=255)
                                                                    cTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGAtttt                                                                             <  4:224987/149‑1 (MQ=255)
                                                                     ttGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTc                                                                            >  3:80832/1‑149 (MQ=255)
                                                                       gcgcTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCt                                                                          <  3:315912/149‑1 (MQ=255)
                                                                        tactgaGCAAGGCATAGGCTTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTg                                                                         >  3:224685/6‑149 (MQ=255)
                                                                         gcTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGc                                                                        >  3:22562/1‑149 (MQ=255)
                                                                         gcTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGc                                                                        >  4:354240/1‑149 (MQ=255)
                                                                          cTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGc                                                                        <  2:445560/148‑1 (MQ=255)
                                                                               cAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACg                                                                   >  1:507667/1‑148 (MQ=255)
                                                                                  aaGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACa                                                                >  4:218383/1‑149 (MQ=255)
                                                                                     ccTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACaaaa                                                             <  4:80832/149‑1 (MQ=255)
                                                                                        aTGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAaaca                                                           <  1:737344/148‑1 (MQ=255)
                                                                                         tGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAaacaa                                                          <  2:507667/148‑1 (MQ=255)
                                                                                            cTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGaa                                                                                                                             <  2:105603/78‑1 (MQ=255)
                                                                                            cTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGaa                                                                                                                             >  1:105603/1‑78 (MQ=255)
                                                                                             ttGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGtt                                                      <  1:559083/148‑1 (MQ=255)
                                                                                               gCTCAGTATGAGTTATATGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACt                                                   <  3:771520/149‑1 (MQ=255)
                                                                                                 tCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGaaa                                                                                                                            <  4:21051/74‑1 (MQ=255)
                                                                                                 tCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGaaa                                                                                                                            >  3:21051/1‑74 (MQ=255)
                                                                                                 tCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTg                                                  <  1:905230/148‑1 (MQ=255)
                                                                                                  cAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGa                                                 >  1:741090/1‑148 (MQ=255)
                                                                                                     tATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGttt                                             >  3:866966/1‑149 (MQ=255)
                                                                                                      aTGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGttt                                             >  1:515169/1‑148 (MQ=255)
                                                                                                                     gTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTaca                             >  3:61721/1‑149 (MQ=255)
                                                                                                                     gTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTaca                             >  4:702788/1‑149 (MQ=255)
                                                                                                                          gATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGcc                        >  3:9555/1‑149 (MQ=255)
                                                                                                                           aTTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGcc                        >  1:261431/1‑148 (MQ=255)
                                                                                                                           aTTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCg                       >  3:137607/1‑149 (MQ=255)
                                                                                                                           aTTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCg                       >  4:184845/1‑149 (MQ=255)
                                                                                                                             tGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGt                      >  2:478501/1‑148 (MQ=255)
                                                                                                                               cTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTcc                   >  4:417673/1‑149 (MQ=255)
                                                                                                                               cTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTcc                   <  4:22562/149‑1 (MQ=255)
                                                                                                                                ttCAATGGTACTTCCCGGCTGGGAAGGGACTTTCTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACGTTACACTGCCGTTccc                  >  3:572473/1‑149 (MQ=255)
                                                                                                                                  cAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCa                 >  1:112596/1‑148 (MQ=255)
                                                                                                                                  cAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCAc                <  3:184845/149‑1 (MQ=255)
                                                                                                                                  cAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCAc                <  3:262829/149‑1 (MQ=255)
                                                                                                                                    aTGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACt               >  1:271447/1‑148 (MQ=255)
                                                                                                                                      ggTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCa             <  2:93789/148‑1 (MQ=255)
                                                                                                                                       tcccTTCCCAGCCGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATc           >  3:136122/4‑149 (MQ=255)
                                                                                                                                       gTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCAt            <  1:890696/148‑1 (MQ=255)
                                                                                                                                        cccTTCCCAGCCGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATc           >  2:719086/3‑148 (MQ=255)
                                                                                                                                         ccTTCCCAGCCGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCg          >  1:607435/2‑148 (MQ=255)
                                                                                                                                         ccTTCCCAGCCGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCg          >  2:906854/2‑148 (MQ=255)
                                                                                                                                          cTTCCCAGCCGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGtt                     <  2:335315/136‑1 (MQ=255)
                                                                                                                                          cTTCCCAGCCGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGtt                     >  1:335315/1‑136 (MQ=255)
                                                                                                                                          cTTCCCAGCCGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTc                    <  3:183290/137‑1 (MQ=255)
                                                                                                                                          cTTCCCAGCCGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCgg         >  1:844130/1‑148 (MQ=255)
                                                                                                                                           ttCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTGGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAATGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCgggg       >  4:555820/1‑149 (MQ=255)
                                                                                                                                           ttCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGa                                                                  >  4:835269/1‑90 (MQ=255)
                                                                                                                                           ttCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCGGCCCACGa                                                                  <  3:835269/90‑1 (MQ=255)
                                                                                                                                            tCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCgggg       >  2:398488/1‑148 (MQ=255)
                                                                                                                                             cccGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCGGGGc      >  2:891229/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GATACTGTTACTCGGGCTGCTGGGTTATCTGCTTTCCGTTTTTTGTTGGTATAAAACGCTTCATAAACTTGCGCTAAGCAAAGCCTATGCCTTGCTCAGTATGAGTTATGTGCTGGTATGGATTGCTTCAATGGTACTTCCCGGCTGGGAAGGGACTTTTTCGCTGAAAGCACTACTGGGAGTAGCCTGTATTATGAGCGGGTTGATGCTGATTTTCCTGCCCACGACAAAACAACGTTACTGAGTTTTCCCTGCCACTTTACACTGCCGTTCCCACTCATCGGGGCTATA  >  NC_000913/2373035‑2373325

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