Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,561,503 A→G intergenic (+366/‑94) b2447 → / → b2448 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,561,5030AG99.8% 1387.0 / NA 425intergenic (+366/‑94)b2447/b2448hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (232/192);  total (232/193)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 435 total aligned reads displayed.

CACCAAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTG  >  NC_000913/2561413‑2561592
                                                                                          |                                                                                         
   caAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTcc                                                                                      <  2:4583821/93‑1 (MQ=255)
    aaaGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTccc                                                                                     >  2:9314515/1‑93 (MQ=255)
     aaGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCca                                                                                    >  1:3238643/1‑93 (MQ=255)
          aacaaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaac                                                                               <  2:1835216/93‑1 (MQ=255)
          aacaaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaac                                                                               <  2:6121485/93‑1 (MQ=255)
           acacGATGGGAGGTGAATTTTGTGATGCAGATCCCTTTTTCCCTTCAGGCATTTTTGTTTAATAAAATCAACGAAAGGCGGTTCCcaacaaca                                                                              <  2:2926161/93‑1 (MQ=255)
           acaaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaaca                                                                              >  2:9188193/1‑93 (MQ=255)
           acaaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaaca                                                                              >  1:2690478/1‑93 (MQ=255)
            caaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaacaa                                                                             >  1:10010658/1‑93 (MQ=255)
            caaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaacaa                                                                             >  2:6432239/1‑93 (MQ=255)
             aaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACaaa                                                                            <  2:7414056/93‑1 (MQ=255)
              aGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACaaaa                                                                           <  1:8114798/93‑1 (MQ=255)
              aGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACaaaa                                                                           <  1:8539151/93‑1 (MQ=255)
                aTGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATa                                                                         >  1:7588772/1‑93 (MQ=255)
                 tGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATaa                                                                        >  1:1164119/1‑93 (MQ=255)
                  gggAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATaaa                                                                       <  2:6696389/93‑1 (MQ=255)
                  gggAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATaaa                                                                       >  2:8236587/1‑93 (MQ=255)
                  gggAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATaaa                                                                       >  2:9695179/1‑93 (MQ=255)
                    gagaTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGa                                                                     <  2:417561/93‑1 (MQ=255)
                    gAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAGCGAAAGGCGGTTCCCAACAACAAAATAAAGa                                                                     >  1:2352819/1‑93 (MQ=255)
                    gAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGa                                                                     <  1:818022/93‑1 (MQ=255)
                    gAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGa                                                                     <  1:5835569/93‑1 (MQ=255)
                    gAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGa                                                                     >  1:5332555/1‑93 (MQ=255)
                      ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATg                                                                   <  2:7125991/93‑1 (MQ=255)
                      ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATg                                                                   <  2:8086485/93‑1 (MQ=255)
                      ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATg                                                                   >  1:9418738/1‑93 (MQ=255)
                      ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATg                                                                   <  1:9051517/93‑1 (MQ=255)
                      ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATg                                                                   >  1:2567345/1‑93 (MQ=255)
                      ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATg                                                                   <  1:4353368/93‑1 (MQ=255)
                          aaTTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAg                                                               >  1:8669472/1‑93 (MQ=255)
                           aTTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAgg                                                              >  2:6631578/1‑93 (MQ=255)
                              ttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTaa                                                           >  1:646901/1‑93 (MQ=255)
                               tgtgATGCAGATCGCTTTTTCGCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTaaa                                                          >  2:9888658/1‑93 (MQ=255)
                                 tgATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaa                                                        <  2:4915507/93‑1 (MQ=255)
                                  gATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAt                                                       <  2:768331/93‑1 (MQ=255)
                                  gATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAt                                                       <  1:5149544/93‑1 (MQ=255)
                                    tGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGa                                                     <  2:1701803/93‑1 (MQ=255)
                                    tGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGa                                                     >  2:9015916/1‑93 (MQ=255)
                                    tGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGa                                                     >  1:4490261/1‑93 (MQ=255)
                                      cAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGaaa                                                   <  2:1467649/93‑1 (MQ=255)
                                       aGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAAc                                                  >  1:5229774/1‑93 (MQ=255)
                                        gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACg                                                 <  1:6347389/93‑1 (MQ=255)
                                         aTCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACgg                                                >  1:8817164/1‑93 (MQ=255)
                                          tCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGt                                               >  1:6206261/1‑93 (MQ=255)
                                           cGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTg                                              <  1:4100836/93‑1 (MQ=255)
                                             cTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTggg                                            >  1:1715408/1‑93 (MQ=255)
                                              tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGa                                           >  2:3358247/1‑93 (MQ=255)
                                              tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGa                                           <  1:673440/93‑1 (MQ=255)
                                               ttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGaa                                          <  1:4873194/93‑1 (MQ=255)
                                               ttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGaa                                          >  1:9162409/1‑93 (MQ=255)
                                                tttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAg                                         >  2:1137495/1‑93 (MQ=255)
                                                tttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAg                                         >  1:5514850/1‑93 (MQ=255)
                                                 ttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAgg                                        >  2:3575777/1‑92 (MQ=255)
                                                  tCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTg                                       <  2:1502815/93‑1 (MQ=255)
                                                   cccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTgg                                      >  1:3612019/1‑93 (MQ=255)
                                                   cccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTgg                                      >  1:6113627/1‑93 (MQ=255)
                                                   cccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTgg                                      >  1:5992095/1‑93 (MQ=255)
                                                    ccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGc                                     <  2:8683741/93‑1 (MQ=255)
                                                      ttCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCtt                                   >  1:3359629/1‑93 (MQ=255)
                                                      ttCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCtt                                   >  1:2039051/1‑93 (MQ=255)
                                                      ttCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCtt                                   <  2:3632508/93‑1 (MQ=255)
                                                       tCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACTGAGGGAAGTGGCttt                                  >  1:7148291/1‑93 (MQ=255)
                                                       tCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCttt                                  >  1:2909504/1‑93 (MQ=255)
                                                        cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTa                                 <  2:5316564/93‑1 (MQ=255)
                                                         aGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTAc                                >  2:6301024/1‑93 (MQ=255)
                                                          ggCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACt                               >  1:896007/1‑93 (MQ=255)
                                                          ggCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACt                               >  1:7953765/1‑93 (MQ=255)
                                                           gCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTg                              <  2:944246/93‑1 (MQ=255)
                                                            cATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTgg                             <  2:7806686/93‑1 (MQ=255)
                                                             aTTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGt                            >  1:2554143/1‑93 (MQ=255)
                                                               ttttGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGa                          <  2:3089576/93‑1 (MQ=255)
                                                                tttGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGCGGATTTACTGGGGag                         >  1:8731828/1‑92 (MQ=255)
                                                                  tGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATgg                       >  2:7182380/1‑93 (MQ=255)
                                                                   gTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGt                      >  2:6944942/1‑93 (MQ=255)
                                                                    tttAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTa                     >  1:6778526/1‑93 (MQ=255)
                                                                     ttAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTaa                    <  1:7074079/93‑1 (MQ=255)
                                                                       aaTTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATg                  >  2:2371440/1‑93 (MQ=255)
                                                                       aaTTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATg                  >  1:7496500/1‑93 (MQ=255)
                                                                        aTTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGt                 <  1:56707/93‑1 (MQ=255)
                                                                        aTTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGt                 <  2:5933535/93‑1 (MQ=255)
                                                                         ttCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                <  1:3876005/93‑1 (MQ=255)
                                                                         ttCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                <  1:311722/93‑1 (MQ=255)
                                                                          tcaatcaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTg               >  1:7563905/1‑93 (MQ=255)
                                                                           caatcaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTgg              <  1:4006086/93‑1 (MQ=255)
                                                                            aatcaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGt             >  1:3453755/1‑93 (MQ=255)
                                                                            aatcaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGt             <  1:9613868/93‑1 (MQ=255)
                                                                            aatcaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGt             >  2:3919658/1‑93 (MQ=255)
                                                                            aatcaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGt             <  1:3922759/93‑1 (MQ=255)
                                                                             atcaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtt            <  1:2788218/93‑1 (MQ=255)
                                                                               caaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTg          >  2:5013535/1‑93 (MQ=255)
                                                                               caaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTg          >  1:9893232/1‑93 (MQ=255)
                                                                                aaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGc         >  2:4750251/1‑93 (MQ=255)
                                                                                 aCGAAATGCGGTTCCCAACAACAAAGTAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCt        <  2:9617671/93‑1 (MQ=255)
                                                                                 aCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCt        <  2:5126599/93‑1 (MQ=255)
                                                                                   gAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGt      >  2:2749331/1‑93 (MQ=255)
                                                                                    aaaGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTa     <  2:2209541/93‑1 (MQ=255)
                                                                                    aaaGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTa     <  1:2183921/93‑1 (MQ=255)
                                                                                    aaaGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTa     >  2:4338385/1‑93 (MQ=255)
                                                                                      aGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAtt   <  2:7606499/93‑1 (MQ=255)
                                                                                       ggCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTg  <  1:4146257/93‑1 (MQ=255)
                                                                                          |                                                                                         
CACCAAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTG  >  NC_000913/2561413‑2561592

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