Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,722,918 T→A N247Y (AAC→TAC)  b1647 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,722,9180TA100.0% 759.8 / NA 197N247Y (AAC→TAC) b1647hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (112/85);  total (112/85)

Only 100 of 204 total aligned reads displayed.

TTTCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTTAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAAC  >  NC_000913/1722828‑1723006
                                                                                          |                                                                                        
 ttCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAac                                                                                       <  1:5758307/93‑1 (MQ=255)
   cTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAacca                                                                                     <  1:2042859/93‑1 (MQ=255)
   cTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAAGCGTAAacca                                                                                     >  2:5896107/1‑93 (MQ=255)
     gCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccacc                                                                                   <  1:5444466/93‑1 (MQ=255)
     gCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccacc                                                                                   <  1:5054474/93‑1 (MQ=255)
     gCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccacc                                                                                   <  2:3926839/93‑1 (MQ=255)
     gCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCCGCGATGGTAAACGTAATACCCAGGCGTAAaccacc                                                                                   >  1:2342416/1‑93 (MQ=255)
        tCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGt                                                                                <  2:4089078/93‑1 (MQ=255)
         cGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAGTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGtt                                                                               <  2:3391860/93‑1 (MQ=255)
          gACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTc                                                                              <  2:2374195/93‑1 (MQ=255)
          gACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTc                                                                              >  2:1275167/1‑93 (MQ=255)
          gACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTc                                                                              >  2:1320597/1‑93 (MQ=255)
          gACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTc                                                                              <  1:3672315/93‑1 (MQ=255)
           aCAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCa                                                                             >  1:3776996/1‑93 (MQ=255)
           aCAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCa                                                                             >  1:1942356/1‑93 (MQ=255)
             aGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAAt                                                                           <  2:5272306/93‑1 (MQ=255)
             aGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAAt                                                                           >  2:3621068/1‑93 (MQ=255)
             aGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAAt                                                                           >  2:6161539/1‑93 (MQ=255)
              gCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATc                                                                          >  2:4260899/1‑93 (MQ=255)
              gCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATc                                                                          <  2:5221697/93‑1 (MQ=255)
                tGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGa                                                                        <  2:2146804/93‑1 (MQ=255)
                 gCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAg                                                                       >  2:2794062/1‑93 (MQ=255)
                   ccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGcc                                                                     <  1:2964913/93‑1 (MQ=255)
                   ccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGcc                                                                     <  2:4883712/93‑1 (MQ=255)
                   ccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGcc                                                                     <  2:3563571/93‑1 (MQ=255)
                   ccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCCAGcc                                                                     >  2:257073/1‑93 (MQ=255)
                        cgcgCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGt                                                                <  1:3215250/93‑1 (MQ=255)
                        cgcgCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGt                                                                >  2:3716844/1‑93 (MQ=255)
                          cgcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGttt                                                              <  2:3073338/93‑1 (MQ=255)
                           gcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTa                                                             >  2:3252924/1‑93 (MQ=255)
                             tCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAAc                                                           >  1:1112947/1‑93 (MQ=255)
                             tCAATTTTACCTGAACCGATAATTGGCAGCCGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAAc                                                           >  1:5287256/1‑93 (MQ=255)
                                aTTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcc                                                        >  2:2684758/1‑93 (MQ=255)
                                 ttttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcct                                                       <  2:3931492/93‑1 (MQ=255)
                                    tACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTg                                                    >  2:3288454/1‑93 (MQ=255)
                                    tACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTg                                                    <  2:247701/93‑1 (MQ=255)
                                    tACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTg                                                    <  1:4954908/93‑1 (MQ=255)
                                     aCCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGc                                                   >  1:670389/1‑93 (MQ=255)
                                     aCCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGc                                                   >  2:5615086/1‑93 (MQ=255)
                                     aCCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAAACACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGc                                                   >  2:1735171/1‑93 (MQ=255)
                                       cTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGcca                                                 >  2:4838956/1‑93 (MQ=255)
                                        tGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccac                                                <  2:1564669/93‑1 (MQ=255)
                                        tGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccac                                                <  1:556361/93‑1 (MQ=255)
                                         gAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccacc                                               >  1:4476065/1‑93 (MQ=255)
                                          aaCCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccacca                                              <  1:53798/93‑1 (MQ=255)
                                           aCCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCcaccac                                             >  1:2107251/1‑93 (MQ=255)
                                           aCCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCcaccac                                             >  1:1041835/1‑93 (MQ=255)
                                            ccGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACcaca                                            <  2:2308595/93‑1 (MQ=255)
                                            ccGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACcaca                                            <  2:4332195/93‑1 (MQ=255)
                                              gATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGc                                          <  1:2500305/93‑1 (MQ=255)
                                               aTAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGcc                                         >  1:4701291/1‑93 (MQ=255)
                                                    tGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtc                                    >  1:5019551/1‑93 (MQ=255)
                                                     ggCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtca                                   <  2:5201636/93‑1 (MQ=255)
                                                     ggCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtca                                   >  2:4466259/1‑93 (MQ=255)
                                                     ggCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtca                                   >  2:1758698/1‑93 (MQ=255)
                                                     ggCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtca                                   >  1:3447844/1‑93 (MQ=255)
                                                      gCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcat                                  <  1:1802715/93‑1 (MQ=255)
                                                         gCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCAc                               <  2:6259247/93‑1 (MQ=255)
                                                         gCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCAc                               >  2:1962617/1‑93 (MQ=255)
                                                         gCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCAc                               >  1:4794961/1‑93 (MQ=255)
                                                           ggCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGc                             <  2:1684007/93‑1 (MQ=255)
                                                           ggCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGc                             >  1:3252886/1‑93 (MQ=255)
                                                            gCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCa                            >  2:6067682/1‑93 (MQ=255)
                                                            gCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCa                            >  2:1481855/1‑93 (MQ=255)
                                                             cTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAg                           >  1:560975/1‑93 (MQ=255)
                                                             cTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAg                           <  2:3993590/93‑1 (MQ=255)
                                                              tGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGc                          >  2:3727797/1‑93 (MQ=255)
                                                                  aTGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCt                      <  2:4429586/93‑1 (MQ=255)
                                                                  aTGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCt                      >  2:4566838/1‑93 (MQ=255)
                                                                  aTGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCt                      >  2:5960228/1‑93 (MQ=255)
                                                                  aTGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCt                      >  2:1507227/1‑93 (MQ=255)
                                                                   tGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg                     >  1:77949/1‑93 (MQ=255)
                                                                   tGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg                     >  2:107361/1‑93 (MQ=255)
                                                                   tGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg                     >  2:3231353/1‑93 (MQ=255)
                                                                    ggTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTgg                    >  2:4954602/1‑93 (MQ=255)
                                                                    ggTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTgg                    >  2:5638921/1‑93 (MQ=255)
                                                                       aatcGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCTCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGaaa                 <  2:4742153/90‑1 (MQ=255)
                                                                          cGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaataa              <  1:3586635/93‑1 (MQ=255)
                                                                          cGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaataa              <  2:6215652/93‑1 (MQ=255)
                                                                            tAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaataatg            <  1:4672046/93‑2 (MQ=255)
                                                                            tAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatc            >  2:2027087/1‑93 (MQ=255)
                                                                            tAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatc            <  2:473184/93‑1 (MQ=255)
                                                                            tAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatc            >  1:641161/1‑93 (MQ=255)
                                                                            tAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatc            <  2:3594224/93‑1 (MQ=255)
                                                                              aTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatcat          >  2:379891/1‑93 (MQ=255)
                                                                              aTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatcat          >  1:418134/1‑93 (MQ=255)
                                                                              aTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatcat          >  1:3990878/1‑93 (MQ=255)
                                                                              aTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatcat          <  2:5412623/93‑1 (MQ=255)
                                                                               tACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatcatc         <  1:2271461/93‑1 (MQ=255)
                                                                                aCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAAtcatca        <  1:5850328/93‑1 (MQ=255)
                                                                                  ccAGGCGTAAACCACCTGTTCAATCGAGGCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGCAATAATCATCAtt      >  1:4021016/1‑93 (MQ=255)
                                                                                  ccAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCAtt      >  1:5033132/1‑93 (MQ=255)
                                                                                  ccAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCAtt      >  2:3877705/1‑93 (MQ=255)
                                                                                   cAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTa     >  2:2974759/1‑93 (MQ=255)
                                                                                    aGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGGAATACTAATCATTaa    >  1:95980/1‑93 (MQ=255)
                                                                                    aGGAGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaa    <  2:5562735/93‑1 (MQ=255)
                                                                                     ggCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa   <  1:4835478/93‑1 (MQ=255)
                                                                                     ggCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa   >  2:5695875/1‑93 (MQ=255)
                                                                                     ggCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa   >  1:628099/1‑93 (MQ=255)
                                                                                      gCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAAc  >  1:1546618/1‑93 (MQ=255)
                                                                                          |                                                                                        
TTTCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTTAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAAC  >  NC_000913/1722828‑1723006

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