Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,594,634 T→C I42V (ATA→GTA)  b2475 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,594,6340TC84.8% 562.8 / 90.6 243I42V (ATA→GTA) b2475hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (16/21);  new base C (105/101);  total (121/122)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 244 total aligned reads displayed.

TTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACAT  >  NC_000913/2594544‑2594726
                                                                                          |                                                                                            
  gACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGt                                                                                          >  1:7289882/1‑93 (MQ=255)
     cgcgTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTaaa                                                                                       >  1:901872/1‑93 (MQ=255)
     cgcgTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTaaa                                                                                       <  2:5338024/93‑1 (MQ=255)
     cgcgTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTaaa                                                                                       <  2:4675076/93‑1 (MQ=255)
      gcgTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTaaaa                                                                                      <  2:5134096/93‑1 (MQ=255)
        gTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATc                                                                                    >  1:3526702/1‑93 (MQ=255)
         ttGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCa                                                                                   <  1:6589218/93‑1 (MQ=255)
         ttGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCa                                                                                   <  2:4573179/93‑1 (MQ=255)
          tGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAg                                                                                  >  2:3840720/1‑93 (MQ=255)
          tGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAg                                                                                  >  2:5456744/1‑93 (MQ=255)
           gATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGc                                                                                 >  2:406631/1‑93 (MQ=255)
            aTTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGcc                                                                                <  2:1839101/93‑1 (MQ=255)
            aTTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGcc                                                                                <  1:3643034/93‑1 (MQ=255)
            aTTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGcc                                                                                <  1:4005309/93‑1 (MQ=255)
             ttgttgttgGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGccc                                                                               <  2:5626134/93‑1 (MQ=255)
               gttgttgGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGc                                                                             >  1:2549586/1‑93 (MQ=255)
                ttgttgGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGcc                                                                            >  2:5772828/1‑93 (MQ=255)
                ttgttgGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGcc                                                                            >  1:1034407/1‑93 (MQ=255)
                  gttgGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCtt                                                                          <  2:5740394/93‑1 (MQ=255)
                   ttgGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCttt                                                                         >  1:3376440/1‑93 (MQ=255)
                    tgGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTa                                                                        <  2:4843455/93‑1 (MQ=255)
                    tgGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTa                                                                        <  2:6089538/93‑1 (MQ=255)
                     gggAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTAc                                                                       <  1:5767303/93‑1 (MQ=255)
                      ggAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTAcc                                                                      >  1:586912/1‑93 (MQ=255)
                       gAACCCGGCTGCCCGGTCATCAACCCGGTAAAATCAACACCATAGTAGCCTGAAACCAGCACGACTCCCAGTAAAATCACCCCGCCTTTACCg                                                                     <  1:5263278/93‑1 (MQ=255)
                         aaCCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCt                                                                   <  1:3146774/93‑1 (MQ=255)
                          aCCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCtt                                                                  >  2:6617927/1‑93 (MQ=255)
                          aCCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCtt                                                                  <  2:808666/93‑1 (MQ=255)
                          aCCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCtt                                                                  >  1:4316923/1‑93 (MQ=255)
                          aCCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCtt                                                                  <  1:3771243/93‑1 (MQ=255)
                          aCCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCtt                                                                  >  1:1288748/1‑93 (MQ=255)
                           ccGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTg                                                                 <  2:4099457/93‑1 (MQ=255)
                            cGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTgg                                                                >  2:4384865/1‑93 (MQ=255)
                             ggCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGc                                                               >  2:4145023/1‑93 (MQ=255)
                              gCTGCCCGGTCATCACCCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCa                                                              <  1:1269101/93‑1 (MQ=255)
                              gCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCa                                                              <  1:4484609/93‑1 (MQ=255)
                               cTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAg                                                             <  1:2207965/93‑1 (MQ=255)
                               cTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAg                                                             >  2:113649/1‑93 (MQ=255)
                                 gCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGc                                                           >  1:6454163/1‑93 (MQ=255)
                                  cccGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCg                                                          <  2:6898422/93‑1 (MQ=255)
                                    cGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCgga                                                        >  1:1204075/1‑91 (MQ=255)
                                       tcatcaACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAc                                                     >  2:3113016/1‑93 (MQ=255)
                                       tcatcaACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAc                                                     <  1:748061/93‑1 (MQ=255)
                                         atcaACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCg                                                   <  2:176219/93‑1 (MQ=255)
                                          tcaACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCgg                                                  >  2:1276110/1‑93 (MQ=255)
                                          tcaACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCgg                                                  <  2:2657545/93‑1 (MQ=255)
                                           caACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCggg                                                 <  2:6835547/93‑1 (MQ=255)
                                             aCCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGAc                                               <  1:4613246/93‑1 (MQ=255)
                                             aCCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGAc                                               <  1:1184107/93‑1 (MQ=255)
                                             aCCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGAc                                               <  1:3243355/93‑1 (MQ=255)
                                             aCCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGAc                                               >  1:596130/1‑93 (MQ=255)
                                             aCCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGAc                                               >  2:3235984/1‑93 (MQ=255)
                                              cccGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGAcc                                              <  2:5929879/93‑1 (MQ=255)
                                               ccGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCg                                             <  1:6743354/93‑1 (MQ=255)
                                               ccGGTTAAATCAACACCATAGTAGCCTCCACCCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCg                                             <  1:2342783/93‑1 (MQ=255)
                                                cGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGc                                            >  1:5419657/1‑93 (MQ=255)
                                                cGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGc                                            <  2:3364764/93‑1 (MQ=255)
                                                cGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGc                                            >  1:5070190/1‑93 (MQ=255)
                                                  gTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGccc                                          >  2:3020062/1‑93 (MQ=255)
                                                   ttAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCa                                         >  2:5223118/1‑93 (MQ=255)
                                                     aaaTCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATa                                       <  1:7070193/93‑1 (MQ=255)
                                                      aaTCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAg                                      >  2:371512/1‑93 (MQ=255)
                                                       aTCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGa                                     <  1:2296650/93‑1 (MQ=255)
                                                       aTCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGa                                     <  2:395377/93‑1 (MQ=255)
                                                       aTCAACACCAAAGTAGCCTGCACCCAGCACAACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGa                                     <  1:3435507/93‑1 (MQ=255)
                                                         cAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATg                                   >  1:1719925/1‑93 (MQ=255)
                                                          aaCACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATgg                                  >  1:4670865/1‑93 (MQ=255)
                                                          aaCACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATgg                                  >  2:3842844/1‑93 (MQ=255)
                                                           acacCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGa                                 >  2:4900000/1‑93 (MQ=255)
                                                            cacCATAGTAGCCTGCAACCAGCATGACTACCAGTAAAATCAGCCCACCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGac                                <  2:4782988/93‑1 (MQ=255)
                                                            cacCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGac                                <  2:304796/93‑1 (MQ=255)
                                                              ccATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGacca                              >  2:4293066/1‑93 (MQ=255)
                                                               cATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGaccac                             <  1:3005036/93‑1 (MQ=255)
                                                                aTAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGaccacc                            >  1:6830749/1‑93 (MQ=255)
                                                                 tagtagCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGAccacca                           <  1:96621/93‑1 (MQ=255)
                                                                  agtagCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCag                          <  1:2974989/93‑1 (MQ=255)
                                                                   gtagCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCaga                         <  1:5608803/93‑1 (MQ=255)
                                                                    tagCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCagag                        >  2:549602/1‑93 (MQ=255)
                                                                      gCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCt                      >  2:2885372/1‑93 (MQ=255)
                                                                        cTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGt                    <  1:2580769/93‑1 (MQ=255)
                                                                        cTGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGt                    <  2:2060944/93‑1 (MQ=255)
                                                                         tGCTACCAGCACGACTACCAGTAAAATCAGCTCGCCTTTGCCGCGTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGtt                   <  2:1500870/93‑1 (MQ=255)
                                                                         tGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGtt                   >  1:5589585/1‑93 (MQ=255)
                                                                         tGCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGtt                   >  2:6214573/1‑93 (MQ=255)
                                                                          gCAACCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTg                  >  2:380317/1‑93 (MQ=255)
                                                                            aaCCAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTgcg                <  2:754357/93‑1 (MQ=255)
                                                                              ccAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGcc              >  2:7204025/1‑93 (MQ=255)
                                                                              ccAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGcc              <  2:5984706/93‑1 (MQ=255)
                                                                              ccAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGcc              >  1:7367887/1‑93 (MQ=255)
                                                                              ccAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGcc              <  2:268922/93‑1 (MQ=255)
                                                                               cAGCACGACTATCAATAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCt             >  2:6444604/1‑93 (MQ=255)
                                                                               cAGCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCt             >  2:5320319/1‑93 (MQ=255)
                                                                                 gCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCATCAGAGCTGTTGCGCCTGt           >  2:2138141/1‑93 (MQ=255)
                                                                                 gCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGt           <  2:2613779/93‑1 (MQ=255)
                                                                                 gCACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGt           >  1:1823176/1‑93 (MQ=255)
                                                                                  cACGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTc          <  2:1815544/93‑1 (MQ=255)
                                                                                    cGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCtt        <  1:4106787/93‑1 (MQ=255)
                                                                                    cGACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCtt        <  2:2411344/93‑1 (MQ=255)
                                                                                       cTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCaa     >  1:7249683/1‑93 (MQ=255)
                                                                                         aTCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACa   >  1:7106307/1‑93 (MQ=255)
                                                                                          |                                                                                            
TTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACAT  >  NC_000913/2594544‑2594726

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