Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 588,485 G→A R824C (CGC→TGC)  nfrA ← bacteriophage N4 receptor, outer membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913588,4850GA100.0% 534.4 / NA 137R824C (CGC→TGC) nfrAbacteriophage N4 receptor, outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (65/72);  total (65/72)

Only 100 of 137 total aligned reads displayed.

CGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGT  >  NC_000913/588339‑588631
                                                                                                                                                  |                                                                                                                                                  
cGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAg                                                                                                                                                 <  1:722581/150‑1 (MQ=255)
 gggTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGa                                                                                                                                                <  1:934942/150‑1 (MQ=255)
 gggTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGAc                                                                                                                                               <  2:1399468/151‑1 (MQ=255)
    tGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCgg                                                                                                                                            >  1:76277/1‑151 (MQ=255)
     gCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCgg                                                                                                                                            <  3:295423/150‑1 (MQ=255)
     gCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGt                                                                                                                                           <  2:353358/151‑1 (MQ=255)
      ccATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTg                                                                                                                                          >  2:761411/1‑151 (MQ=255)
      ccATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTg                                                                                                                                          >  4:232096/1‑151 (MQ=255)
       cATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGc                                                                                                                                         >  1:1385083/1‑151 (MQ=255)
        aTTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGCTCGCGCAGCGGCTTCCAGCACAGACCGGTGc                                                                                                                                         >  1:1685767/1‑150 (MQ=255)
        aTTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGcc                                                                                                                                        >  1:1522677/1‑151 (MQ=255)
         ttCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCg                                                                                                                                       >  1:51131/1‑151 (MQ=255)
         tgcgtcgCTGTATTTGCCGCCATTAAAGAATGAAGCGCTGGCGCGCAGCATGGTAACGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGAAGCGGCTTCCAGCACAGACCGGTGcc                                                                                                                                        <  1:300379/148‑1 (MQ=255)
           agtcgCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCgg                                                                                                                                      >  2:945047/2‑150 (MQ=255)
            gtcgCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGAc                                                                                                                                    >  2:94462/1‑151 (MQ=255)
            gtcgCTGTATTCGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCGGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCgg                                                                                                                                      <  1:1055992/149‑1 (MQ=255)
               gCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATc                                                                                                                                 >  1:1260721/1‑151 (MQ=255)
               gCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATc                                                                                                                                 <  2:1647428/151‑1 (MQ=255)
                  gTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCgg                                                                                                                               >  1:1135890/1‑150 (MQ=255)
                    aTTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGAtt                                                                                                                            >  2:1083535/1‑151 (MQ=255)
                       tGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGt                                                                                                                                           <  1:330375/133‑1 (MQ=255)
                       tGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGAttttt                                                                                                                         >  3:306764/1‑151 (MQ=255)
                       tGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGAttttt                                                                                                                         >  2:956035/1‑151 (MQ=255)
                       tGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGAtttt                                                                                                                          <  1:1117740/150‑1 (MQ=255)
                          cgccATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCAc                                                                                                                      >  2:366245/1‑151 (MQ=255)
                              aTTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGc                                                                                                                  >  2:1166116/1‑151 (MQ=255)
                               ttAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGc                                                                                                                  >  1:854212/1‑150 (MQ=255)
                               ttAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGc                                                                                                                  >  1:1072040/1‑150 (MQ=255)
                                     aaTGAGGCGCTGGCGCGCATCATGGTATCGCATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGGGCCGGACATCGGATTTTTCACCGGcatcat                                                                                                             >  1:1283565/1‑149 (MQ=255)
                                     aaTGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCatcatc                                                                                                            >  2:715276/1‑150 (MQ=255)
                                     aaTGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCAtcatca                                                                                                           <  3:518578/151‑1 (MQ=255)
                                           gcgcTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGt                                                                                                     <  2:713916/151‑1 (MQ=255)
                                           gcgcTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGt                                                                                                     >  1:662391/1‑151 (MQ=255)
                                            cgcTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGtt                                                                                                    <  2:1667832/151‑1 (MQ=255)
                                             gcTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGttt                                                                                                   >  2:1069472/1‑151 (MQ=255)
                                             gcTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGttt                                                                                                   >  1:446067/1‑151 (MQ=255)
                                                ggCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTtc                                                                                                 <  2:317489/150‑1 (MQ=255)
                                                   gcgcAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCg                                                                                             <  2:1126187/151‑1 (MQ=255)
                                                   gcgcAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCg                                                                                             <  2:241909/151‑1 (MQ=255)
                                                    cgcAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCgg                                                                                            >  2:250/1‑151 (MQ=255)
                                                    cgcAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCgg                                                                                            >  2:1384198/1‑151 (MQ=255)
                                                       agcaTGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTAt                                                                                         <  2:1479613/151‑1 (MQ=255)
                                                          aTGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCgg                                                                                      <  2:1647899/151‑1 (MQ=255)
                                                              tATCGGATACGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCaa                                                                                   >  1:1713939/1‑150 (MQ=255)
                                                               aTCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCaaa                                                                                  <  1:308469/150‑1 (MQ=255)
                                                               aTCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCaa                                                                                   >  2:992491/1‑149 (MQ=255)
                                                                 cGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGAc                                                                               <  2:764001/151‑1 (MQ=255)
                                                                 cGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGAc                                                                               <  2:764717/151‑1 (MQ=255)
                                                                     tGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCgg                                                                           >  1:1160605/1‑151 (MQ=255)
                                                                      gcgcCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCgg                                                                           <  2:474566/150‑1 (MQ=255)
                                                                       cgcCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCt                                                                         >  2:43095/1‑151 (MQ=255)
                                                                           aTTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCtata                                                                      >  1:914980/1‑150 (MQ=255)
                                                                           aTTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATaa                                                                     <  3:97227/151‑1 (MQ=255)
                                                                             tttGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATaaa                                                                    <  1:234012/150‑1 (MQ=255)
                                                                              ttGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACt                                                                  >  1:1280558/1‑151 (MQ=255)
                                                                              ttGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACt                                                                  >  1:646923/1‑151 (MQ=255)
                                                                              ttGGCCGTTAAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCGCCGGTATCGGCAAAGACGCGGCTATAAAc                                                                   >  1:432533/1‑150 (MQ=255)
                                                                               tGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACt                                                                  <  2:531607/150‑1 (MQ=255)
                                                                                 gCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTgag                                                               <  2:1311489/151‑1 (MQ=255)
                                                                                   cGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTgag                                                               >  1:1135/1‑149 (MQ=255)
                                                                                   cGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGagag                                                             <  1:1424586/151‑1 (MQ=255)
                                                                                          cGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCg                                                      >  4:551947/1‑151 (MQ=255)
                                                                                                  gctgTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTcc                                                <  2:1195786/149‑1 (MQ=255)
                                                                                                  gctgTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCcag                                              <  2:111940/151‑1 (MQ=255)
                                                                                                  gctgTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCcag                                              <  2:1410487/151‑1 (MQ=255)
                                                                                                  gctgTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCca                                               <  1:257617/150‑1 (MQ=255)
                                                                                                    tgTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCcagc                                             >  1:1168714/1‑150 (MQ=255)
                                                                                                     gTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCcagca                                            <  2:446141/150‑1 (MQ=255)
                                                                                                     gTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCcagc                                             >  2:1542393/1‑149 (MQ=255)
                                                                                                      ttcgggagCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCagcagc                                          >  1:305683/8‑151 (MQ=255)
                                                                                                        cGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAgcagca                                         <  2:156430/150‑1 (MQ=255)
                                                                                                         gACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAgcagca                                         <  2:966679/149‑1 (MQ=255)
                                                                                                         gACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCata                                       <  2:665849/151‑1 (MQ=255)
                                                                                                         gACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCata                                       <  1:1284345/151‑1 (MQ=255)
                                                                                                          aCGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCatat                                      >  1:1181087/1‑151 (MQ=255)
                                                                                                          aCGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCatat                                      <  2:288022/151‑1 (MQ=255)
                                                                                                          aCGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCatat                                      <  2:837433/151‑1 (MQ=255)
                                                                                                            ggCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTg                                    <  2:1028331/151‑1 (MQ=255)
                                                                                                                  gAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTcc                               <  2:835759/150‑1 (MQ=255)
                                                                                                                  gAAAAAGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCg                              >  2:42716/1‑151 (MQ=255)
                                                                                                                     aaaGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATg                           >  2:1717418/1‑151 (MQ=255)
                                                                                                                     aaaGATCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATg                           >  2:366384/1‑151 (MQ=255)
                                                                                                                          tCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGt                       >  1:1046426/1‑150 (MQ=255)
                                                                                                                          tCTGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGt                       >  2:1375860/1‑150 (MQ=255)
                                                                                                                            tGATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACt                    >  4:378850/1‑151 (MQ=255)
                                                                                                                             gATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGCCATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGCAGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGAGTCCGATGCGGTACTc                   >  2:1271073/1‑151 (MQ=255)
                                                                                                                             gATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACt                    <  2:963196/150‑1 (MQ=255)
                                                                                                                             gATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTc                   >  4:409089/1‑151 (MQ=255)
                                                                                                                             gATCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTc                   <  1:1240981/151‑1 (MQ=255)
                                                                                                                               tCGCGCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCgg                 >  4:414938/1‑151 (MQ=255)
                                                                                                                                cgcgCAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGc                <  2:862778/151‑1 (MQ=255)
                                                                                                                                 gcgcAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTTCGTCGGATGCGGTAcgaggc                <  1:286084/150‑6 (MQ=255)
                                                                                                                                  cgcAGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCtt              <  3:162123/151‑1 (MQ=255)
                                                                                                                                     aGCGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCa           <  1:738399/151‑1 (MQ=255)
                                                                                                                                       cGGCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGt         <  1:955914/151‑1 (MQ=255)
                                                                                                                                        ggCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTGCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGtt        <  1:1608174/151‑1 (MQ=255)
                                                                                                                                         gCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTg       >  3:372096/1‑151 (MQ=255)
                                                                                                                                         gCTTCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTg       <  4:548598/151‑1 (MQ=255)
                                                                                                                                           ttCCAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTc     >  1:816048/1‑151 (MQ=255)
                                                                                                                                              cAGCACAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGt  <  1:316818/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGT  >  NC_000913/588339‑588631

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