Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,749,710 T→C L30P (CTT→CCT)  yiaU → putative transcriptional regulator LYSR‑type

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,749,7100TC98.3% 379.0 / ‑2.3 119L30P (CTT→CCT) yiaUputative transcriptional regulator LYSR‑type
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (113/4);  total (115/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.99e-01

Only 100 of 122 total aligned reads displayed.

GATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGC  >  NC_000913/3749621‑3749799
                                                                                         |                                                                                         
gATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTgg                                                                                        >  2:3672186/1‑93 (MQ=255)
 aTGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGc                                                                                       >  2:668792/1‑93 (MQ=255)
 aTGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGc                                                                                       >  1:1989415/1‑93 (MQ=255)
  tGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCa                                                                                      >  2:4639785/1‑93 (MQ=255)
   gACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCAt                                                                                     >  1:3499471/1‑93 (MQ=255)
    aCAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATc                                                                                    >  2:3180662/1‑93 (MQ=255)
    aCAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCACGACTGTACCTGGCAtg                                                                                    >  1:528784/1‑92 (MQ=255)
     cAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCg                                                                                   >  2:4391257/1‑93 (MQ=255)
      aaaaTTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGc                                                                                  >  1:2529275/1‑93 (MQ=255)
      aaaaTTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGc                                                                                  >  2:2499525/1‑93 (MQ=255)
       aaaTTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGca                                                                                 >  1:616164/1‑93 (MQ=255)
       aaaTTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGca                                                                                 >  2:820140/1‑93 (MQ=255)
             cAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGcc                                                                           >  2:108381/1‑93 (MQ=255)
                   aaGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTc                                                                     >  1:2735632/1‑93 (MQ=255)
                         cGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCaaa                                                               >  2:801493/1‑93 (MQ=255)
                         cGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCaaa                                                               >  1:87158/1‑93 (MQ=255)
                          gggAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAat                                                              >  2:2648432/1‑93 (MQ=255)
                          gggAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAat                                                              >  2:2648430/1‑93 (MQ=255)
                          gggAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAat                                                              >  2:4270784/1‑93 (MQ=255)
                          gggAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAat                                                              >  2:1802074/1‑93 (MQ=255)
                          gggAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAat                                                              >  1:2108025/1‑93 (MQ=255)
                          gggAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAat                                                              >  1:2049879/1‑93 (MQ=255)
                           ggAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAata                                                             >  2:1360193/1‑93 (MQ=255)
                            gAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAatat                                                            >  2:2514872/1‑93 (MQ=255)
                              gTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCt                                                          >  2:1967813/1‑93 (MQ=255)
                               ttAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCtt                                                         >  1:3791053/1‑93 (MQ=255)
                               ttAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCtt                                                         >  2:3282731/1‑93 (MQ=255)
                               ttAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCtt                                                         >  1:4172838/1‑93 (MQ=255)
                                 aaaaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGc                                                       >  2:1984143/1‑93 (MQ=255)
                                  aaaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCCACGTCAGCAAATATCTTGCt                                                      >  1:4615173/1‑93 (MQ=255)
                                  aaaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCt                                                      >  1:432963/1‑93 (MQ=255)
                                    aaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAACTATCTTGCTGa                                                    >  1:3023522/1‑93 (MQ=255)
                                    aaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGa                                                    >  1:3611192/1‑93 (MQ=255)
                                    aaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGa                                                    <  2:3669625/93‑1 (MQ=255)
                                    aaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGa                                                    >  1:1620202/1‑93 (MQ=255)
                                     aTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAt                                                   >  1:2244464/1‑93 (MQ=255)
                                     aTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAt                                                   >  1:1199893/1‑93 (MQ=255)
                                     aTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAt                                                   >  2:1030955/1‑93 (MQ=255)
                                       tATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAttt                                                 >  1:2518196/1‑93 (MQ=255)
                                        aTCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAtttt                                                >  2:3295379/1‑93 (MQ=255)
                                        aTCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAtttt                                                >  1:3428675/1‑93 (MQ=255)
                                         tctcGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTg                                               >  2:2855521/1‑93 (MQ=255)
                                         tctcGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTg                                               >  1:1110706/1‑93 (MQ=255)
                                          ctcGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGa                                              >  1:2327751/1‑93 (MQ=255)
                                          ctcGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGa                                              >  2:840632/1‑93 (MQ=255)
                                           tcGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGCTTTTGaa                                             >  1:3792535/1‑93 (MQ=255)
                                           tcGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGaa                                             >  2:3515394/1‑93 (MQ=255)
                                           tcGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGaa                                             >  1:1193861/1‑93 (MQ=255)
                                           tcGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGaa                                             >  1:923964/1‑93 (MQ=255)
                                            cGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAAt                                            >  1:4416058/1‑93 (MQ=255)
                                               tAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCaa                                         >  1:2029656/1‑93 (MQ=255)
                                                aaTCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCaaa                                        >  1:2966548/1‑93 (MQ=255)
                                                aaTCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCaaa                                        >  1:4269458/1‑93 (MQ=255)
                                                aaTCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCaaa                                        >  2:3632765/1‑93 (MQ=255)
                                                 aTCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCaaaa                                       >  2:4053164/1‑93 (MQ=255)
                                                 aTCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCaaaa                                       >  2:914333/1‑93 (MQ=255)
                                                   cGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGt                                     >  1:1784714/1‑93 (MQ=255)
                                                      tGCCAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTggg                                  >  2:3248403/1‑93 (MQ=255)
                                                         cAGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGttt                               >  2:4592808/1‑93 (MQ=255)
                                                          aGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTa                              >  1:2422504/1‑93 (MQ=255)
                                                          aGTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTa                              >  1:2956806/1‑93 (MQ=255)
                                                           gTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTaa                             >  2:1947361/1‑93 (MQ=255)
                                                           gTGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTaa                             >  1:4232390/1‑93 (MQ=255)
                                                            tGAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTaaa                            <  2:4449660/93‑1 (MQ=255)
                                                             gAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGACTCAAGAGTGGGTTTaaaa                           >  2:3055827/1‑93 (MQ=255)
                                                             gAAAATATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTaaaa                           >  1:1628508/1‑93 (MQ=255)
                                                               aaaTATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGt                         >  2:2556142/1‑93 (MQ=255)
                                                                aaTATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTc                        <  2:3334207/93‑1 (MQ=255)
                                                                aaTATCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGCTTAAAAGTc                        >  1:475110/1‑93 (MQ=255)
                                                                 atatCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCt                       >  2:522367/1‑93 (MQ=255)
                                                                 atatCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCt                       >  1:845149/1‑93 (MQ=255)
                                                                  tatCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCtt                      >  1:3028372/1‑93 (MQ=255)
                                                                  tatCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCtt                      >  2:3845438/1‑93 (MQ=255)
                                                                  tatCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCtt                      >  1:817284/1‑93 (MQ=255)
                                                                  tatCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCtt                      >  2:2781069/1‑93 (MQ=255)
                                                                   atCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCttt                     >  2:2953208/1‑93 (MQ=255)
                                                                    tCAGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTg                    >  2:674036/1‑93 (MQ=255)
                                                                      aGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGAc                  >  1:1643239/1‑93 (MQ=255)
                                                                      aGCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGAc                  >  2:1638604/1‑93 (MQ=255)
                                                                       gCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGAcc                 >  1:3166502/1‑93 (MQ=255)
                                                                       gCCATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTGTGAATCAAAAGTGGGTTTACAAGTCTTTGAcc                 >  1:1743078/1‑93 (MQ=255)
                                                                        ccATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCg                >  2:3246967/1‑93 (MQ=255)
                                                                         cATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGa               >  2:1918345/1‑93 (MQ=255)
                                                                         cATGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGa               >  1:1085176/1‑93 (MQ=255)
                                                                          aTGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGaa              >  1:4528743/1‑93 (MQ=255)
                                                                          aTGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGaa              >  2:681400/1‑93 (MQ=255)
                                                                           tGCCGCGACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAac             >  2:255911/1‑93 (MQ=255)
                                                                            gCCGCGACTGTACCTGGCATCGCACAGGCCACCGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCgcaac            >  1:1004797/1‑89 (MQ=255)
                                                                               gcgACTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGGTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAaccacc         >  2:4152536/1‑93 (MQ=255)
                                                                                 gACTGTACCTGGCATCGCACAGGCCAACGACAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGGTTAAAAGTCTTTGACCGAACCACCCg       >  2:36901/1‑93 (MQ=255)
                                                                                  agtGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCgg      >  2:1206597/3‑93 (MQ=255)
                                                                                  aCTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCgg      >  1:4395890/1‑93 (MQ=255)
                                                                                  aCTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCgg      >  1:469841/1‑93 (MQ=255)
                                                                                  aCTGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCgc      >  1:3281072/1‑92 (MQ=255)
                                                                                    tGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCa    >  2:3896635/1‑93 (MQ=255)
                                                                                    tGTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCa    >  1:3335566/1‑93 (MQ=255)
                                                                                     gTACTTGGCCTCGCACAGGCCAACGTAAGCAAATATCTTTCTGATTTTGAATCAAAAGGGGGCTTAAAAGTCTTTTTCTGCACCACCCGGGAg   >  2:3431089/1‑93 (MQ=255)
                                                                                     gTACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAg   >  1:4345460/1‑93 (MQ=255)
                                                                                      tACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGc  >  2:4446801/1‑93 (MQ=255)
                                                                                      tACCTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGc  >  1:4569250/1‑93 (MQ=255)
                                                                                         |                                                                                         
GATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGC  >  NC_000913/3749621‑3749799

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