Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,587,161 T→C I289V (ATT→GTT)  b1505 ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,587,1610TC100.0% 365.7 / NA 114I289V (ATT→GTT) b1505putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (53/61);  total (53/61)

Only 100 of 117 total aligned reads displayed.

CAATTGCACCGAATGGAACGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAATAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAAGAGCAACACGGTTTT  >  NC_000913/1587071‑1587250
                                                                                          |                                                                                         
cAATTGCACCGAATGGAACGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAg                                                                                         <  2:190151/93‑1 (MQ=255)
  aTTGCACCGAATGGAACGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCa                                                                                       >  2:1342316/1‑93 (MQ=255)
  aTTGCACCGAATGGAACGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAAAAgaa                                                                                       >  2:316767/1‑91 (MQ=255)
    tGCACCGAATGGAACGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCAcc                                                                                     <  2:2852021/93‑1 (MQ=255)
     tccccGAATGGAACGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCg                                                                                    <  1:1256333/90‑1 (MQ=255)
       aCCGAATGGAACGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTg                                                                                  <  1:640639/93‑1 (MQ=255)
       aCCGAATGGAACGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTg                                                                                  >  1:3673349/1‑93 (MQ=255)
        ccGAATGGAACGCTCTTATTCCCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGa                                                                                 <  2:4675651/93‑1 (MQ=255)
          gAATGGTACGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGt                                                                               >  2:742911/1‑93 (MQ=255)
           aaTGGACCGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGtt                                                                              <  2:2149410/93‑1 (MQ=255)
           aaTGGAACGCTCTTATTCCCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGtt                                                                              <  2:1049015/93‑1 (MQ=255)
           aaTGGAACGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGtt                                                                              >  2:1555062/1‑93 (MQ=255)
           aaTGGAACGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGtt                                                                              >  2:768193/1‑93 (MQ=255)
             tGGAACGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTgg                                                                            >  2:145433/1‑93 (MQ=255)
             tGGAACGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTgg                                                                            >  1:3573862/1‑93 (MQ=255)
               gAACGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGa                                                                          >  2:3545568/1‑93 (MQ=255)
                aaCGCTCTTATTCCCGTACTTCAACGTCATTAATCCTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGAt                                                                         <  2:122872/93‑1 (MQ=255)
                 aCGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATg                                                                        >  2:2477759/1‑93 (MQ=255)
                  cGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGa                                                                       <  2:4206436/93‑1 (MQ=255)
                  cGCTCTTATTACCGTACTTCAACGTCATTAATACTTCCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGa                                                                       <  2:2719432/93‑1 (MQ=255)
                   gCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGAc                                                                      >  1:824890/1‑93 (MQ=255)
                   gCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGAc                                                                      >  1:3179676/1‑93 (MQ=255)
                    ctctTATTCCCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACa                                                                     <  2:1592868/93‑1 (MQ=255)
                    ctctTATACCCGTTCTACAACGTCTTTAATCCTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACa                                                                     <  2:3912342/93‑1 (MQ=255)
                     tctTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAg                                                                    >  1:506041/1‑93 (MQ=255)
                      tttattaCCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGt                                                                   <  2:4064061/92‑1 (MQ=255)
                        tattaCCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGa                                                                 >  2:3970579/1‑93 (MQ=255)
                         attaCCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGAc                                                                >  2:3295794/1‑93 (MQ=255)
                          ttaCCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGAcc                                                               >  1:2194933/1‑93 (MQ=255)
                           tcccGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCa                                                              <  1:159774/91‑1 (MQ=255)
                            aCCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCAc                                                             <  1:3741970/93‑1 (MQ=255)
                             ccGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACg                                                            >  1:4075885/1‑93 (MQ=255)
                               gTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGt                                                          >  1:3706355/1‑93 (MQ=255)
                               gTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGt                                                          <  1:3081892/93‑1 (MQ=255)
                               gTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGCGACCACGGt                                                          >  2:1335129/1‑93 (MQ=255)
                               gTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAATTGGGATGACAGTGACCACGGt                                                          >  2:315681/1‑93 (MQ=255)
                                tACTTCACCGTCATTAATACTTCCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGtt                                                         <  2:2442671/93‑1 (MQ=255)
                                tACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGtt                                                         >  1:3182888/1‑93 (MQ=255)
                                tACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGtt                                                         <  1:4172663/93‑1 (MQ=255)
                                  cTTCAACGTCATTAATCCTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTc                                                       <  2:1223117/93‑1 (MQ=255)
                                   ttCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCa                                                      >  2:327877/1‑93 (MQ=255)
                                   ttCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCa                                                      >  1:1139459/1‑93 (MQ=255)
                                   ttCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCa                                                      >  1:24047/1‑93 (MQ=255)
                                     cAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATc                                                    >  1:856289/1‑93 (MQ=255)
                                       aCGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCa                                                  <  2:2936070/93‑1 (MQ=255)
                                         gTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGt                                                >  2:1926104/1‑93 (MQ=255)
                                         gTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGt                                                <  1:218363/93‑1 (MQ=255)
                                          tCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGtt                                               >  1:2501150/1‑93 (MQ=255)
                                           cATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTc                                              <  2:4614605/93‑1 (MQ=255)
                                           cATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGCTc                                              >  2:3268081/1‑93 (MQ=255)
                                            aTTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCa                                             >  1:4132516/1‑93 (MQ=255)
                                             ttAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCaa                                            <  1:234999/93‑1 (MQ=255)
                                             ttAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCaa                                            >  1:3734815/1‑93 (MQ=255)
                                               aaTCCTTCCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACa                                          <  2:2872920/93‑1 (MQ=255)
                                               aaTACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACa                                          <  2:1445027/93‑1 (MQ=255)
                                                attcTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACAt                                         <  1:4362980/90‑1 (MQ=255)
                                                 tgcTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACAtt                                        <  2:4153234/91‑1 (MQ=255)
                                                 tACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACAtt                                        <  2:821772/93‑1 (MQ=255)
                                                 tACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACAtt                                        >  1:4524318/1‑93 (MQ=255)
                                                 tACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACAtt                                        >  1:4462388/1‑93 (MQ=255)
                                                  ccTTCCCCGCCGAGATGTGTTGTAAATGTTGCTCTGGCAACATTACCGTGGCTTGGGATGACAGTGACCGCGGTTTCATCCAGTTCAACATTa                                       <  2:4241370/92‑1 (MQ=255)
                                                  aCTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTa                                       >  2:648344/1‑93 (MQ=255)
                                                  aCTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTa                                       <  1:334603/93‑1 (MQ=255)
                                                    tttCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATc                                     <  1:1570389/93‑1 (MQ=255)
                                                    tttCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATc                                     >  1:2872181/1‑93 (MQ=255)
                                                     tccccGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCt                                    <  2:207938/91‑1 (MQ=255)
                                                     tccccGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCt                                    <  2:2101137/91‑1 (MQ=255)
                                                      tCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGATTCATCCAGTTCAACATTATCTg                                   >  1:3865733/1‑93 (MQ=255)
                                                         cgccgATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCaa                                <  2:3282589/93‑1 (MQ=255)
                                                         cgccgATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCaa                                <  1:1810562/93‑1 (MQ=255)
                                                           ccgATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAgg                              <  1:736212/93‑1 (MQ=255)
                                                           ccgATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAgg                              <  2:2123547/93‑1 (MQ=255)
                                                           ccgATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAgg                              >  1:3015454/1‑93 (MQ=255)
                                                             gATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGa                            <  2:4211219/93‑1 (MQ=255)
                                                             gATTTGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGa                            <  2:1913837/93‑1 (MQ=255)
                                                               tttGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAAt                          >  1:1507995/1‑93 (MQ=255)
                                                               tttGTGCATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAAt                          >  2:1607732/1‑93 (MQ=255)
                                                                     cATTAAATGTTGCTCTGGCAACATCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGt                    >  1:3767232/1‑93 (MQ=255)
                                                                     cATTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGt                    >  2:4363150/1‑93 (MQ=255)
                                                                      aTTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGtt                   <  1:4037005/93‑1 (MQ=255)
                                                                      aTTAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGtt                   >  1:3238331/1‑93 (MQ=255)
                                                                       ttAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTa                  >  1:4359566/1‑93 (MQ=255)
                                                                       ttAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTa                  <  1:869808/93‑1 (MQ=255)
                                                                        tAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTaa                 >  1:4044398/1‑93 (MQ=255)
                                                                        tAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTaa                 >  1:4002642/1‑93 (MQ=255)
                                                                        tAAATGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTaa                 <  1:4094820/93‑1 (MQ=255)
                                                                         aaaTGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAag                <  2:1154997/93‑1 (MQ=255)
                                                                          aaTGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAaga               <  1:2820301/93‑1 (MQ=255)
                                                                          aaTGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAaga               >  2:3429740/1‑93 (MQ=255)
                                                                            tGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAAGAGc             <  2:1886869/93‑1 (MQ=255)
                                                                            tGTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAAGAGc             >  2:4021760/1‑93 (MQ=255)
                                                                             gTTGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAAGAGCa            >  2:2146364/1‑93 (MQ=255)
                                                                              ttGCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAAGAGCaa           <  2:4433365/93‑1 (MQ=255)
                                                                                gCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAAGAGCAaca         >  2:3791547/1‑93 (MQ=255)
                                                                                gCTCTGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAAGAGCAaca         >  1:3845361/1‑93 (MQ=255)
                                                                                    tGGCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAAGAGCAACACGGt     <  1:3440448/93‑1 (MQ=255)
                                                                                      gCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAAGAGCAACACGGttt   <  1:2132440/93‑1 (MQ=255)
                                                                                      gCAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAAGAGCAACACGGttt   <  1:4211234/93‑1 (MQ=255)
                                                                                       cAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAAGAGCAACACGGtttt  <  2:1667419/93‑1 (MQ=255)
                                                                                       cAACAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAAGAGCAACACGGtttt  >  1:3469337/1‑93 (MQ=255)
                                                                                          |                                                                                         
CAATTGCACCGAATGGAACGCTCTTATTACCGTACTTCAACGTCATTAATACTTTCCCGCCGATTTGTGCATTAAATGTTGCTCTGGCAATAGCACCGTGAGTTGGGATGACAGTGACCACGGTTTCATCCAGTTCAACATTATCTGCAAGGGAATTCGCGTTAAGAGCAACACGGTTTT  >  NC_000913/1587071‑1587250

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