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% 325.0 / NA 92N247Y (AAC→TAC) b1647hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (48/44);  total (48/44)

TTCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTTAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACA  >  NC_000913/1722829‑1723007
                                                                                         |                                                                                         
ttCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAac                                                                                        <  1:4622745/93‑1 (MQ=255)
ttCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAac                                                                                        >  2:855163/1‑93 (MQ=255)
 tCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAacc                                                                                       <  1:2958784/93‑1 (MQ=255)
   tGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccac                                                                                     >  2:744101/1‑93 (MQ=255)
   tGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccac                                                                                     <  1:3126566/93‑1 (MQ=255)
   tGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccac                                                                                     <  2:2241222/93‑1 (MQ=255)
    gCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccacc                                                                                    <  2:853757/93‑1 (MQ=255)
     cTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCt                                                                                   <  1:3053945/93‑1 (MQ=255)
     cTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCt                                                                                   >  2:258005/1‑93 (MQ=255)
      ttCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTg                                                                                  <  1:1041949/93‑1 (MQ=255)
       tCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGt                                                                                 <  2:3950014/93‑1 (MQ=255)
       tCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGt                                                                                 <  1:1683132/93‑1 (MQ=255)
         caCAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTc                                                                               <  1:3408262/92‑1 (MQ=255)
          aCAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCa                                                                              <  2:692729/93‑1 (MQ=255)
          aCAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCa                                                                              >  2:3000135/1‑93 (MQ=255)
           cAGCTGCCCGTACGCGCTCAATTTTCCCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCaa                                                                             <  2:1429139/93‑1 (MQ=255)
           cAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCaa                                                                             >  2:3901557/1‑93 (MQ=255)
            aGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGCTCAAt                                                                            >  2:4395442/1‑93 (MQ=255)
             gCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATc                                                                           >  2:857795/1‑93 (MQ=255)
             gCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATc                                                                           >  2:3389788/1‑93 (MQ=255)
               tGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGa                                                                         <  2:3842828/93‑1 (MQ=255)
               tGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGa                                                                         <  1:4583382/93‑1 (MQ=255)
                  ccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGcc                                                                      <  2:1831138/93‑1 (MQ=255)
                  ccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGcc                                                                      <  1:1389196/93‑1 (MQ=255)
                   cGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGccc                                                                     >  1:1095893/1‑93 (MQ=255)
                   cGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGccc                                                                     >  1:1430537/1‑93 (MQ=255)
                    gTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcg                                                                    >  1:698007/1‑93 (MQ=255)
                    gTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcg                                                                    >  2:1068078/1‑93 (MQ=255)
                     tACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcgc                                                                   >  2:33892/1‑93 (MQ=255)
                     tACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcgc                                                                   >  2:942937/1‑93 (MQ=255)
                      aCGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcgcg                                                                  <  2:3867155/93‑1 (MQ=255)
                        gcgcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGtt                                                                >  2:4570877/1‑93 (MQ=255)
                          gcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTa                                                              >  1:4617623/1‑93 (MQ=255)
                          gcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTa                                                              <  1:562286/93‑1 (MQ=255)
                          gcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTa                                                              >  2:4366461/1‑93 (MQ=255)
                           cTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTaa                                                             <  2:3024624/93‑1 (MQ=255)
                           cTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTaa                                                             <  1:891498/93‑1 (MQ=255)
                            tCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAAc                                                            >  1:2105457/1‑93 (MQ=255)
                            tCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAAc                                                            >  1:2643179/1‑93 (MQ=255)
                            tCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAAc                                                            <  2:569705/93‑1 (MQ=255)
                             cAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAACCCGCGTTTaagt                                                           >  1:523366/1‑91 (MQ=255)
                              aaTTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActc                                                          <  1:3633016/93‑1 (MQ=255)
                               aTTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcc                                                         >  2:1785598/1‑93 (MQ=255)
                               aTTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcc                                                         >  1:697417/1‑93 (MQ=255)
                                 tttACCTGAACCGATAATTGGCGGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcctc                                                       <  2:1788362/93‑1 (MQ=255)
                                 tttACCTGAACCGATAATTGGCAGCGGCTGCGATTGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcctc                                                       <  2:2754961/93‑1 (MQ=255)
                                  ttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCctct                                                      >  2:2998359/1‑93 (MQ=255)
                                    aCCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGc                                                    >  1:742971/1‑93 (MQ=255)
                                    aCCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGc                                                    >  1:3273814/1‑93 (MQ=255)
                                           ccGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACcaca                                             <  2:3084134/93‑1 (MQ=255)
                                            cGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAg                                            >  1:3549332/1‑93 (MQ=255)
                                             gATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGc                                           <  2:3185926/93‑1 (MQ=255)
                                               tAATTGGCAGCGGCTGCCATGGTAACCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTATCTCCTCTGCCACCACAGCCa                                         <  1:3805489/93‑1 (MQ=255)
                                                 aTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGt                                       >  1:2337736/1‑93 (MQ=255)
                                                  ttGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGtt                                      >  2:623479/1‑93 (MQ=255)
                                                  ttGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGtt                                      >  1:102493/1‑93 (MQ=255)
                                                  ttGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGtt                                      >  2:639632/1‑93 (MQ=255)
                                                    ggCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtca                                    <  2:4034358/93‑1 (MQ=255)
                                                     gCAGCGGCTGCGATGGTAAACGTAATCCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcat                                   <  1:1809886/93‑1 (MQ=255)
                                                      cAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatc                                  <  2:2787518/93‑1 (MQ=255)
                                                        gCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTACTATGCCACCACATCCAGTTCATCAc                                >  2:947507/1‑93 (MQ=255)
                                                             tGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGc                           <  1:378201/93‑1 (MQ=255)
                                                             tGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGc                           >  2:4082916/1‑93 (MQ=255)
                                                             tGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGc                           <  1:696161/93‑1 (MQ=255)
                                                                gATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGc                        >  2:4393923/1‑93 (MQ=255)
                                                                  tGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg                      >  1:2419002/1‑93 (MQ=255)
                                                                  tGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg                      <  2:1940462/93‑1 (MQ=255)
                                                                  tGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg                      >  2:4364388/1‑93 (MQ=255)
                                                                   ggTAAACGTAATCCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTgg                     <  2:3768133/93‑1 (MQ=255)
                                                                     tAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGaa                   >  2:3379002/1‑93 (MQ=255)
                                                                     tAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGaa                   >  2:382854/1‑93 (MQ=255)
                                                                     tAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGaa                   >  2:4011878/1‑93 (MQ=255)
                                                                      aaaCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGaaa                  >  1:1202859/1‑93 (MQ=255)
                                                                      aaaCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGaaa                  >  2:4257461/1‑93 (MQ=255)
                                                                      aaaCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGaaa                  >  2:1145410/1‑93 (MQ=255)
                                                                       aaCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaat                 >  1:3695848/1‑93 (MQ=255)
                                                                       aaCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaat                 >  1:3840787/1‑93 (MQ=255)
                                                                       aaCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaat                 <  1:4483147/93‑1 (MQ=255)
                                                                        aCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaata                <  2:1871301/93‑1 (MQ=255)
                                                                        aCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaata                <  2:3800093/93‑1 (MQ=255)
                                                                         cGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaataa               <  1:1752027/93‑1 (MQ=255)
                                                                          gTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaataat              <  2:2639397/93‑1 (MQ=255)
                                                                          gTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaataat              <  2:2390097/93‑1 (MQ=255)
                                                                          gTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaataat              >  2:2272579/1‑93 (MQ=255)
                                                                          gTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGCTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaataat              >  1:1758796/1‑93 (MQ=255)
                                                                           tAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatc             <  2:900698/93‑1 (MQ=255)
                                                                            aaTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatca            >  1:4641477/1‑93 (MQ=255)
                                                                              atcccAGCCATAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAAAAatcatc          <  2:4454019/91‑1 (MQ=255)
                                                                               aCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAAtcatca         >  1:159560/1‑93 (MQ=255)
                                                                               aCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAAtcatca         <  1:865419/93‑1 (MQ=255)
                                                                                cccAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATcatcat        >  1:4081602/1‑93 (MQ=255)
                                                                                  cAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTa      >  1:3890184/1‑93 (MQ=255)
                                                                                   aGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaa     <  2:474569/93‑1 (MQ=255)
                                                                                     gCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAAc   >  1:524083/1‑93 (MQ=255)
                                                                                      cGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACa  <  2:4444115/93‑1 (MQ=255)
                                                                                      cGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACa  <  2:4319610/93‑1 (MQ=255)
                                                                                         |                                                                                         
TTCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTTAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACA  >  NC_000913/1722829‑1723007

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