Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,765,916 A→G Y38C (TAT→TGT)  yibG → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,765,9160AG98.9% 911.3 / ‑2.2 277Y38C (TAT→TGT) yibGhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (129/145);  total (132/145)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 284 total aligned reads displayed.

TTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTATGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTATATCAATATTACT  >  NC_000913/3765826‑3766005
                                                                                          |                                                                                         
ttttttATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCtgtg                                                                                         >  2:164508/1‑93 (MQ=255)
ttttttATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCtgtg                                                                                         <  1:2212156/93‑1 (MQ=255)
ttttttATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCtgtg                                                                                         <  1:2103749/93‑1 (MQ=255)
ttttttATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCtgtg                                                                                         <  2:472197/93‑1 (MQ=255)
 tttttATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGa                                                                                        >  1:2429015/1‑93 (MQ=255)
 tttttATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGa                                                                                        <  1:725750/93‑1 (MQ=255)
  ttttATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGAt                                                                                       <  1:576468/93‑1 (MQ=255)
  ttttATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGAt                                                                                       >  1:3319961/1‑93 (MQ=255)
  ttttATTTCTCTTTTATTTGTCAATTGAATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGAt                                                                                       <  1:2064704/93‑1 (MQ=255)
   tttATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGAtt                                                                                      <  2:4489540/93‑1 (MQ=255)
   tttATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGAtt                                                                                      <  2:3678371/93‑1 (MQ=255)
    ttATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGAttt                                                                                     >  1:4070143/1‑93 (MQ=255)
      aTTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTaa                                                                                   >  1:1707867/1‑93 (MQ=255)
       tttCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAg                                                                                  >  1:1689359/1‑93 (MQ=255)
       tttCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAg                                                                                  <  2:2982284/93‑1 (MQ=255)
         tctcTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTg                                                                                <  1:2882905/93‑1 (MQ=255)
          ctctTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTga                                                                               <  2:110997/93‑1 (MQ=255)
           tctTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTgaa                                                                              <  2:4063936/93‑1 (MQ=255)
           tctTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTgaa                                                                              <  1:3071254/93‑1 (MQ=255)
           tctTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTgaa                                                                              >  1:3674133/1‑93 (MQ=255)
             ttttATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTgaaga                                                                            <  2:3926168/93‑1 (MQ=255)
              tttATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTgaagaa                                                                           <  2:3800429/93‑1 (MQ=255)
              tttATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTgaagaa                                                                           >  1:1252092/1‑93 (MQ=255)
              tttATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTgaagaa                                                                           >  2:2389369/1‑93 (MQ=255)
               ttATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGaaa                                                                          <  2:577860/93‑1 (MQ=255)
               ttATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGaaa                                                                          >  2:196140/1‑93 (MQ=255)
                tATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGaaaa                                                                         <  2:3638266/93‑1 (MQ=255)
                 aTTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGaaaat                                                                        >  1:4082232/1‑92 (MQ=255)
                 aTTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAg                                                                        <  1:1151345/93‑1 (MQ=255)
                  tttGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGt                                                                       >  2:3400655/1‑93 (MQ=255)
                  tttGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGt                                                                       >  2:829380/1‑93 (MQ=255)
                     gTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTtga                                                                    >  2:633895/1‑93 (MQ=255)
                      tCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTtgat                                                                   <  1:2359153/93‑1 (MQ=255)
                      tCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGCATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTtgat                                                                   >  1:4022487/1‑93 (MQ=255)
                         aTTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGaa                                                                >  1:150093/1‑93 (MQ=255)
                         aTTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGaa                                                                <  1:860569/93‑1 (MQ=255)
                         aTTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGaa                                                                <  1:474570/93‑1 (MQ=255)
                         aTTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGaa                                                                <  1:1016924/93‑1 (MQ=255)
                          ttGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGaaa                                                               <  1:1036824/93‑1 (MQ=255)
                          ttGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGaaa                                                               >  1:4261033/1‑93 (MQ=255)
                           tGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAAc                                                              <  1:3477094/93‑1 (MQ=255)
                           tGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAAc                                                              >  1:3392706/1‑93 (MQ=255)
                             cATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTg                                                            >  2:1512058/1‑93 (MQ=255)
                              aTGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGa                                                           >  1:2697577/1‑93 (MQ=255)
                               tGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGaa                                                          <  2:2098556/93‑1 (MQ=255)
                                  tGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGta                                                       <  2:1475412/93‑1 (MQ=255)
                                     tGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATaa                                                    >  2:668492/1‑93 (MQ=255)
                                      gATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATaaa                                                   >  1:4642032/1‑93 (MQ=255)
                                      gATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATaaa                                                   <  1:2823755/93‑1 (MQ=255)
                                        agtgAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAACCTGAAGTATAAATc                                                 <  1:2250119/92‑1 (MQ=255)
                                         gtgAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATca                                                >  2:10432/1‑93 (MQ=255)
                                          tgAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATcac                                               >  2:359539/1‑93 (MQ=255)
                                          tgAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATcac                                               >  1:98579/1‑93 (MQ=255)
                                            aaaaTAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCacaa                                             <  1:63409/93‑1 (MQ=255)
                                            aaaaTAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCacaa                                             <  1:1265948/93‑1 (MQ=255)
                                             aaaTAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCacaac                                            <  1:1105225/93‑1 (MQ=255)
                                             aaaTAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCacaac                                            >  2:3889459/1‑93 (MQ=255)
                                               aTAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACAt                                          <  2:3722490/93‑1 (MQ=255)
                                               aTAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTATGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACAt                                          >  2:3991773/1‑93 (MQ=255)
                                                 aaaacaaaacGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATgg                                        >  2:3255691/1‑93 (MQ=255)
                                                  aaacaaaacGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGt                                       >  2:1224682/1‑93 (MQ=255)
                                                      aaaacGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGAtt                                   >  2:192834/1‑93 (MQ=255)
                                                       aaacGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTc                                  >  2:2319747/1‑93 (MQ=255)
                                                       aaacGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTc                                  <  2:749866/93‑1 (MQ=255)
                                                       aaacGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTc                                  >  1:1528294/1‑93 (MQ=255)
                                                        aacGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCa                                 >  1:1334219/1‑93 (MQ=255)
                                                         acGAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCag                                <  2:2054120/93‑1 (MQ=255)
                                                           gAAAGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCagag                              >  2:4442256/1‑93 (MQ=255)
                                                             aaGTATGATGGGGTCTACAGCAATGACCTGTGATTTGAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGc                            >  1:2651256/1‑93 (MQ=255)
                                                              aGTATGATGGGGTCTACAGCAATGNCCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGca                           <  2:4594834/93‑1 (MQ=255)
                                                              aGTATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGca                           <  2:1509978/93‑1 (MQ=255)
                                                                tATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGcatc                         >  1:1777095/1‑93 (MQ=255)
                                                                tATGATGGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGcatc                         <  1:4327157/93‑1 (MQ=255)
                                                                 atgatgGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGcatca                        >  2:1589753/1‑93 (MQ=255)
                                                                    atgGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCAttt                     <  1:2979097/93‑1 (MQ=255)
                                                                     tgGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTc                    <  2:2684679/93‑1 (MQ=255)
                                                                     tgGGGTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTc                    <  1:2952486/93‑1 (MQ=255)
                                                                       gggTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGa                  >  2:1012505/1‑93 (MQ=255)
                                                                       gggTCTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGa                  >  1:1610015/1‑93 (MQ=255)
                                                                           cTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACtat              >  2:3114695/1‑93 (MQ=255)
                                                                           cTACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACtat              >  1:3757886/1‑93 (MQ=255)
                                                                            tACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACtata             <  2:2535224/93‑1 (MQ=255)
                                                                            tACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACtata             >  2:3423287/1‑93 (MQ=255)
                                                                            tACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACtata             >  2:2897528/1‑93 (MQ=255)
                                                                            tACAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACtata             >  2:1971345/1‑93 (MQ=255)
                                                                             aCAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGACGTATAAATCACAACATGATGATTCAGAGGCATCATTTCGACTatat            >  1:129161/1‑93 (MQ=255)
                                                                             aCAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTatat            >  2:274610/1‑93 (MQ=255)
                                                                              cAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTATATc           <  2:2700031/93‑1 (MQ=255)
                                                                              cAGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTATATc           <  2:4172755/93‑1 (MQ=255)
                                                                               aGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTATATCa          >  2:3502274/1‑93 (MQ=255)
                                                                               aGCAATGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTATATCa          <  1:4415125/93‑1 (MQ=255)
                                                                                acAAGGACCTTTGATTTAAGGGAAAAAAGGTTGATGAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTATATCAat        <  1:1176858/92‑1 (MQ=255)
                                                                                    tGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTATATCAATAtt     >  2:1791985/1‑93 (MQ=255)
                                                                                    tGACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTATATCAATAtt     <  1:3860018/93‑1 (MQ=255)
                                                                                     gACCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTATATCAATATTa    <  1:1094539/93‑1 (MQ=255)
                                                                                      aCCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTATATCAATATTAc   >  2:4275325/1‑93 (MQ=255)
                                                                                      aCCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTATATCAATATTAc   <  2:219510/93‑1 (MQ=255)
                                                                                      aCCTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTATATCAATATTAc   <  1:4663435/93‑1 (MQ=255)
                                                                                       ccTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTATATCAATATTACt  <  1:446584/93‑1 (MQ=255)
                                                                                       ccTGTGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTATATCAATATTACt  <  1:1870541/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTCTACAGCAATGACCTATGATTTAAGTGAAGAAAAGTTGATGAAACTGAAGTATAAATCACAACATGGTGATTCAGAGGCATCATTTCGACTATATCAATATTACT  >  NC_000913/3765826‑3766005

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