Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,808,078 A→G intergenic (+6/‑157) yniC → / → b1728 putative phosphatase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,808,0780AG98.5% 411.7 / ‑2.9 135intergenic (+6/‑157)yniC/b1728putative phosphatase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (67/66);  total (69/66)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

Only 100 of 142 total aligned reads displayed.

AGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACAC  >  NC_000913/1807988‑1808166
                                                                                          |                                                                                        
aGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCggg                                                                                        <  1:3988391/93‑1 (MQ=255)
 gcgcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCgggg                                                                                       >  2:2137729/1‑93 (MQ=255)
 acgcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCgggg                                                                                       <  2:2279291/92‑1 (MQ=255)
  cgcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGgggg                                                                                      <  1:935248/93‑1 (MQ=255)
    cAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCa                                                                                    >  2:3574309/1‑93 (MQ=255)
    cAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCa                                                                                    <  1:1322661/93‑1 (MQ=255)
     aaaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAg                                                                                   >  2:1323610/1‑93 (MQ=255)
     aaaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAg                                                                                   >  1:3877077/1‑93 (MQ=255)
      aaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGt                                                                                  >  1:851359/1‑93 (MQ=255)
      aaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGt                                                                                  <  1:532441/93‑1 (MQ=255)
        aTGATCCGCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCGGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGa                                                                                <  2:4096780/93‑1 (MQ=255)
        aTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGa                                                                                >  1:1041713/1‑93 (MQ=255)
        aTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTGCTCGGTTAATGACCGGGGGCAGTGa                                                                                <  2:1906579/93‑1 (MQ=255)
         tGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAt                                                                               <  1:3485205/93‑1 (MQ=255)
          gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtc                                                                              >  1:313907/1‑93 (MQ=255)
          gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtc                                                                              >  2:4275030/1‑93 (MQ=255)
           aTCCGCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtct                                                                             <  2:3412515/93‑1 (MQ=255)
            tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtctc                                                                            >  1:926554/1‑93 (MQ=255)
              cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGc                                                                          <  1:2129615/93‑1 (MQ=255)
                   tttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGccc                                                                     >  2:3239845/1‑93 (MQ=255)
                   tttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGccc                                                                     >  2:4501819/1‑93 (MQ=255)
                    ttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGAACGGGGGCAGTGATCTCGCTGcccc                                                                    <  2:4519755/93‑1 (MQ=255)
                     tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCt                                                                   >  2:597875/1‑93 (MQ=255)
                     tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCt                                                                   >  2:2215114/1‑93 (MQ=255)
                      gTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTg                                                                  <  1:1727405/93‑1 (MQ=255)
                       tATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTgg                                                                 >  1:1397413/1‑93 (MQ=255)
                       tATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTgg                                                                 >  2:3533858/1‑93 (MQ=255)
                        aTTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGt                                                                <  1:2612795/93‑1 (MQ=255)
                        aTTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGt                                                                >  2:2620306/1‑93 (MQ=255)
                        aTTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGt                                                                >  2:3279696/1‑93 (MQ=255)
                         ttAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGtt                                                               >  1:2719384/1‑93 (MQ=255)
                          tAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTc                                                              <  2:1858240/93‑1 (MQ=255)
                          tAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTc                                                              <  1:2619987/93‑1 (MQ=255)
                            gCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTATTGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                            <  2:4613057/93‑1 (MQ=255)
                            gCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                            >  1:730242/1‑93 (MQ=255)
                              aGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTa                                                          >  2:394094/1‑93 (MQ=255)
                              aGACGTAAAACTGTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTa                                                          <  1:2146688/93‑1 (MQ=255)
                               gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                         >  2:1479161/1‑93 (MQ=255)
                               gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                         <  2:735523/93‑1 (MQ=255)
                               gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                         >  2:35094/1‑93 (MQ=255)
                               gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                         <  1:3841654/93‑1 (MQ=255)
                               gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGGTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                         >  2:4357832/1‑93 (MQ=255)
                               aaCGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                         <  1:4352349/92‑1 (MQ=255)
                                aCGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATc                                                        <  1:1747742/93‑1 (MQ=255)
                                      aaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGCAAt                                                  >  2:2903457/1‑93 (MQ=255)
                                      aaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                  >  2:2997602/1‑93 (MQ=255)
                                      aaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGCGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                  >  1:2047076/1‑93 (MQ=255)
                                       aCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAtt                                                 <  1:536151/93‑1 (MQ=255)
                                       aCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAtt                                                 >  1:4270319/1‑93 (MQ=255)
                                          ttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTCCCCCTGGTTCTTTCTCTGAAATTGCa                                              >  1:252736/1‑93 (MQ=255)
                                             aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTc                                           >  1:1211719/1‑93 (MQ=255)
                                             aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCCTTATCTGAAATTGCATTc                                           >  2:778112/1‑93 (MQ=255)
                                              tCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCa                                          >  1:4548209/1‑93 (MQ=255)
                                              tCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCa                                          >  1:3106207/1‑93 (MQ=255)
                                                 cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACt                                       >  1:2489628/1‑93 (MQ=255)
                                                 cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACt                                       <  1:2330381/93‑1 (MQ=255)
                                                  tGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTg                                      >  2:218688/1‑93 (MQ=255)
                                                  tGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTACACTg                                      >  2:387970/1‑93 (MQ=255)
                                                       gAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGa                                 <  2:1656555/93‑1 (MQ=255)
                                                       gAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGa                                 >  2:715045/1‑93 (MQ=255)
                                                         aCTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAtt                               >  2:1586163/1‑93 (MQ=255)
                                                         aCTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAtt                               <  2:90970/93‑1 (MQ=255)
                                                         aCTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAtt                               <  2:2887611/93‑1 (MQ=255)
                                                         aCTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAtt                               <  2:3049841/93‑1 (MQ=255)
                                                            cACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAAt                            <  2:616992/93‑1 (MQ=255)
                                                            cACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTACATTCAACTGACGGATTAAt                            >  1:3946849/1‑93 (MQ=255)
                                                             aCCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                           >  1:295501/1‑93 (MQ=255)
                                                             aCCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                           <  1:1533242/93‑1 (MQ=255)
                                                             aCCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                           >  1:1932331/1‑93 (MQ=255)
                                                              ccGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCg                          <  1:390403/93‑1 (MQ=255)
                                                              ccGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCg                          <  1:3832978/93‑1 (MQ=255)
                                                                gCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                        <  2:3419366/93‑1 (MQ=255)
                                                                gCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                        >  1:873329/1‑93 (MQ=255)
                                                                 cAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                       <  1:4353188/93‑1 (MQ=255)
                                                                 cAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                       >  1:670704/1‑93 (MQ=255)
                                                                  aaaaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCaa                      >  1:4166616/1‑93 (MQ=255)
                                                                  aaaaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCaa                      >  1:3170898/1‑93 (MQ=255)
                                                                  aaaaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCaa                      >  1:3152330/1‑93 (MQ=255)
                                                                  aaaaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCaa                      <  2:3785596/93‑1 (MQ=255)
                                                                   aaaGACCTTCTCGGTTAATGACCGGGGGCAGTCATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAt                     <  2:2974406/93‑1 (MQ=255)
                                                                    aaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAtt                    <  2:595678/93‑1 (MQ=255)
                                                                    aaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAtt                    <  1:4672824/93‑1 (MQ=255)
                                                                      gACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTa                  >  2:1055336/1‑93 (MQ=255)
                                                                       aCCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTaa                 >  1:555622/1‑93 (MQ=255)
                                                                          ttCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAagag              >  2:1565552/1‑93 (MQ=255)
                                                                          ttCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAagag              >  1:2752012/1‑93 (MQ=255)
                                                                           tctcGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAgaga             >  2:2579943/1‑93 (MQ=255)
                                                                           tctcGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAgaga             <  1:2711290/93‑1 (MQ=255)
                                                                            ctcGGTTAATGACCGGGTGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaa            <  1:2366626/93‑1 (MQ=255)
                                                                             tcGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaaa           <  1:974379/93‑1 (MQ=255)
                                                                             tcGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaaa           >  2:351102/1‑93 (MQ=255)
                                                                             tcGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaaa           >  2:1756206/1‑93 (MQ=255)
                                                                               ggTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga         <  1:3793244/93‑1 (MQ=255)
                                                                               ggTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga         <  2:2019421/93‑1 (MQ=255)
                                                                               ggTCAATGACCGGGGGCAGTGATCTCGCTGCCCCTGTTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga         <  2:3571506/93‑1 (MQ=255)
                                                                                gTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag        >  1:1810761/1‑93 (MQ=255)
                                                                                gTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag        >  2:2574682/1‑93 (MQ=255)
                                                                                 ttAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGt       <  1:225826/93‑1 (MQ=255)
                                                                                  tAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGtt      <  1:4219687/93‑1 (MQ=255)
                                                                                     tGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTaca   >  1:3503411/1‑93 (MQ=255)
                                                                                          |                                                                                        
AGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACAC  >  NC_000913/1807988‑1808166

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