Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,839,484 A→G intergenic (‑51/‑30) b1758 ← / → b1759 putative cytochrome oxidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,839,4840AG99.3% 952.0 / ‑3.5 296intergenic (‑51/‑30)b1758/b1759putative cytochrome oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (155/139);  total (156/140)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 298 total aligned reads displayed.

GTTTAATCCGGGGATGAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCATTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTCGCG  >  NC_000913/1839394‑1839573
                                                                                          |                                                                                         
gTTTAATCCGGGGATGAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGtt                                                                                         >  1:8595329/1‑93 (MQ=255)
 tttAATCCGGGGATGAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGttt                                                                                        >  1:2754278/1‑93 (MQ=255)
    aaTCCGGGGATGAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTc                                                                                     <  2:7110787/93‑1 (MQ=255)
     aTCCGTGGATGAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCg                                                                                    <  1:3248336/93‑1 (MQ=255)
     aTCCGGGGATGAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCg                                                                                    >  2:8430645/1‑93 (MQ=255)
      tCCGGGGATGAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGc                                                                                   <  2:1652251/93‑1 (MQ=255)
        cGGGGATGAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCtt                                                                                 <  2:4299375/93‑1 (MQ=255)
         ggggATGAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTa                                                                                <  1:1495567/93‑1 (MQ=255)
         ggggATGAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTa                                                                                <  1:9221209/93‑1 (MQ=255)
            gATGAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCt                                                                             <  1:4507873/93‑1 (MQ=255)
             aTGAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTg                                                                            >  2:5801042/1‑93 (MQ=255)
              tGAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGCTATTATGGCGTTTTTCGCTTATCTGc                                                                           >  1:6688824/1‑93 (MQ=255)
              tGAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGc                                                                           >  2:8565147/1‑93 (MQ=255)
               gAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGcc                                                                          >  2:412520/1‑93 (MQ=255)
               gAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGcc                                                                          <  2:6084253/93‑1 (MQ=255)
                  gATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCgt                                                                       >  2:5892062/1‑93 (MQ=255)
                   aTGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCgtg                                                                      <  1:7078420/93‑1 (MQ=255)
                    tGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCgtgt                                                                     <  1:9146012/93‑1 (MQ=255)
                    tGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCATTTTTCGCTTATCTGCCCgtgt                                                                     >  1:9795097/1‑93 (MQ=255)
                         gTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAAt                                                                >  1:7214176/1‑93 (MQ=255)
                         gTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAAt                                                                >  1:4824158/1‑93 (MQ=255)
                            tAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTa                                                             >  1:7058131/1‑93 (MQ=255)
                             aGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTAt                                                            >  2:7549894/1‑93 (MQ=255)
                              gCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATg                                                           >  2:8867562/1‑93 (MQ=255)
                               cacaGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGa                                                          <  2:10027759/93‑1 (MQ=255)
                               cacaGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGa                                                          <  1:8193632/93‑1 (MQ=255)
                               cacaGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGa                                                          <  2:1888155/93‑1 (MQ=255)
                                acaGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGaa                                                         <  1:8116316/93‑1 (MQ=255)
                                 caGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCgtg                                                                      <  1:6196255/79‑1 (MQ=255)
                                 caGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCgtg                                                                      >  2:6196254/1‑79 (MQ=255)
                                 caGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGaaa                                                        <  1:585868/93‑1 (MQ=255)
                                 caGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGaaa                                                        >  2:4596754/1‑93 (MQ=255)
                                  aGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGaaaa                                                       <  1:5448365/93‑1 (MQ=255)
                                  aGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGaaaa                                                       <  1:5242654/93‑1 (MQ=255)
                                     cACTCCTTAAATATAACGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGa                                                    <  1:8972262/93‑1 (MQ=255)
                                     cACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGa                                                    <  2:2329910/93‑1 (MQ=255)
                                      aCTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGAt                                                   <  1:3980846/93‑1 (MQ=255)
                                       cTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCTTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGAtt                                                  <  2:6880749/93‑1 (MQ=255)
                                       cTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGAtt                                                  >  1:6360545/1‑93 (MQ=255)
                                       cTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGAtt                                                  <  2:4355452/93‑1 (MQ=255)
                                       cTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGAtt                                                  >  2:5375770/1‑93 (MQ=255)
                                        tCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTg                                                 >  2:9675074/1‑93 (MQ=255)
                                        tCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTg                                                 >  2:8387192/1‑93 (MQ=255)
                                        tCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTg                                                 <  1:4705506/93‑1 (MQ=255)
                                         ccTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGa                                                >  1:6566569/1‑93 (MQ=255)
                                          cTTAAATATAAAGCCTTTCTGATTTAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGaa                                               <  1:2587659/93‑1 (MQ=255)
                                          cTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGaa                                               <  1:10039549/93‑1 (MQ=255)
                                          cTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGaa                                               <  2:3539131/93‑1 (MQ=255)
                                           ttAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAg                                              <  1:7287727/93‑1 (MQ=255)
                                            tAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAgt                                             >  2:6339770/1‑93 (MQ=255)
                                            tAAACATAGAGCCTTTCTGATTTAGGAACCGTGCGGATATTATGGCTTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAgt                                             <  1:2139986/93‑1 (MQ=255)
                                             taaTATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAgtt                                            <  1:7355566/92‑1 (MQ=255)
                                               atatAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAgttgt                                          <  1:1816120/93‑1 (MQ=255)
                                                tataAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAgttgtt                                         >  2:7238320/1‑93 (MQ=255)
                                                tataAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAgttgtt                                         >  2:7184916/1‑93 (MQ=255)
                                                aataAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAgttgtt                                         <  1:6933913/92‑1 (MQ=255)
                                                 ataAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGttgttg                                        >  2:8673269/1‑93 (MQ=255)
                                                 ataAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGttgttg                                        >  2:2692361/1‑93 (MQ=255)
                                                  taAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTgc                                       <  2:8210432/93‑1 (MQ=255)
                                                  taAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTgc                                       >  2:1337683/1‑93 (MQ=255)
                                                  taAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTgc                                       <  2:9052608/93‑1 (MQ=255)
                                                   aaaGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTgcc                                      <  2:7739324/93‑1 (MQ=255)
                                                   aaaGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTgcc                                      <  2:5231387/93‑1 (MQ=255)
                                                     aGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTgccgc                                    <  1:7973778/93‑1 (MQ=255)
                                                       ccTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCca                                  <  1:1934501/93‑1 (MQ=255)
                                                        cTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCcat                                 <  2:1181547/93‑1 (MQ=255)
                                                        cTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCcat                                 >  2:2064701/1‑93 (MQ=255)
                                                        cTTTCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGGTGTTGCCGCcat                                 >  2:2764645/1‑93 (MQ=255)
                                                         tttCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCcatc                                >  1:7411435/1‑93 (MQ=255)
                                                         tttCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCcatc                                >  2:2160959/1‑93 (MQ=255)
                                                          ttCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCcatca                               <  2:9236830/93‑1 (MQ=255)
                                                          ttCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCcatca                               <  2:8779642/93‑1 (MQ=255)
                                                          ttCTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCcatca                               <  2:4829623/93‑1 (MQ=255)
                                                            cTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCAtt                             >  2:3849080/1‑93 (MQ=255)
                                                            cTGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCAtt                             >  1:6439600/1‑93 (MQ=255)
                                                             tGATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTg                            >  1:5045298/1‑93 (MQ=255)
                                                              gATTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCATGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGa                           >  2:3063513/1‑93 (MQ=255)
                                                               aTTGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGaa                          <  1:8401481/93‑1 (MQ=255)
                                                                ttGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAAc                         >  1:455273/1‑93 (MQ=255)
                                                                ttGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAAc                         >  1:1638307/1‑93 (MQ=255)
                                                                ttGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAAc                         >  1:7480926/1‑93 (MQ=255)
                                                                ttGAGCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAAc                         >  1:5126749/1‑93 (MQ=255)
                                                                    gCAACAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGa                     <  2:741690/93‑1 (MQ=255)
                                                                       aCAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATgg                  <  1:9344668/93‑1 (MQ=255)
                                                                       aCAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATgg                  >  2:4904426/1‑93 (MQ=255)
                                                                       aCAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATgg                  >  2:7985057/1‑93 (MQ=255)
                                                                        cAGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGc                 >  1:7731777/1‑93 (MQ=255)
                                                                         aGTGCGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCa                >  1:564788/1‑93 (MQ=255)
                                                                             cGGATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAAt            <  1:4249398/93‑1 (MQ=255)
                                                                               gATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAAttt          >  2:1942841/1‑93 (MQ=255)
                                                                               gATATTATGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAAttt          >  1:6752513/1‑93 (MQ=255)
                                                                                  attatGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTTCCGCCATCATTGAACGTGATGGCAAAATTTTAc       <  2:5602503/93‑1 (MQ=255)
                                                                                  attatGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTAc       >  2:2526749/1‑93 (MQ=255)
                                                                                   ttatGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACt      <  2:3254653/93‑1 (MQ=255)
                                                                                   ttatGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGAAAATTTTAcct     >  1:3341161/1‑91 (MQ=255)
                                                                                    tatGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTc     >  1:2387272/1‑93 (MQ=255)
                                                                                     atGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTcg    <  2:3662872/93‑1 (MQ=255)
                                                                                      tGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTcgc   <  2:1375146/93‑1 (MQ=255)
                                                                                      tGGCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTcgc   >  2:147694/1‑93 (MQ=255)
                                                                                       ggCGTTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTcgcg  >  1:3013417/1‑93 (MQ=255)
                                                                                          |                                                                                         
GTTTAATCCGGGGATGAAGATGGCGGTCTAGCACAGGCACTCCTTAAATATAAAGCCTTTCTGATTGAGCAACAGTGCGGATATTATGGCATTTTTCGCTTATCTGCCCGTGTGTAATTTATGAAAATGATTGAAGTTGTTGCCGCCATCATTGAACGTGATGGCAAAATTTTACTCGCG  >  NC_000913/1839394‑1839573

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