Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,725,518 (G)6→5 coding (836/1098 nt) nemA → N‑ethylmaleimide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,725,5130G.98.9% 806.9 / ‑0.2 179coding (831/1098 nt)nemAN‑ethylmaleimide reductase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (95/82);  total (97/82)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 211 total aligned reads displayed.

GGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTA  >  NC_000913/1725421‑1725600
                                                                                            |                                                                                       
 gCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgg                                                                                        >  1:1883357/1‑93 (MQ=255)
 gCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgg                                                                                        >  2:1145783/1‑93 (MQ=255)
 gCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgg                                                                                        >  2:1733801/1‑93 (MQ=255)
  cccGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCggg                                                                                       >  1:1160914/1‑93 (MQ=255)
  cccGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCggg                                                                                       <  1:6092066/93‑1 (MQ=255)
   ccGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggg                                                                                      <  1:3668575/93‑1 (MQ=255)
   cagAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggg                                                                                      <  2:3610410/91‑1 (MQ=255)
    cGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCggggg                                                                                     <  1:1035490/93‑1 (MQ=255)
    cGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCggggg                                                                                     <  1:5321821/93‑1 (MQ=255)
     gAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggggt                                                                                    >  2:2919863/1‑92 (MQ=255)
      aaTGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggggtg                                                                                   >  2:721617/1‑91 (MQ=255)
      aaTGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggggtg                                                                                   <  1:5935610/93‑3 (MQ=255)
       aTGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggggtga                                                                                  >  2:2580675/1‑90 (MQ=255)
        tGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggggtgaa                                                                                 >  1:3140724/1‑89 (MQ=255)
        tGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggggtgaa                                                                                 >  2:3418554/1‑89 (MQ=255)
         gaagaaGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAAc                                                                               >  2:3523970/1‑93 (MQ=255)
         gaagaaGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAAc                                                                               >  1:3039686/1‑93 (MQ=255)
          aagaagCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAAcc                                                                              >  2:3200281/1‑93 (MQ=255)
             aagCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGta                                                                           <  1:1260909/93‑1 (MQ=255)
             aagCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGta                                                                           <  2:5780923/93‑1 (MQ=255)
             aagCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGta                                                                           >  1:4989506/1‑93 (MQ=255)
                ccGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATAc                                                                        <  2:2127079/93‑1 (MQ=255)
                ccGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATAc                                                                        <  2:5943786/93‑1 (MQ=255)
                 cGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACt                                                                       <  2:3262805/93‑1 (MQ=255)
                 cGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACt                                                                       <  2:22020/93‑1 (MQ=255)
                   aTGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGa                                                                     <  1:4612128/93‑1 (MQ=255)
                    tGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGAt                                                                    >  1:3564963/1‑93 (MQ=255)
                     gCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATg                                                                   <  1:5453452/93‑1 (MQ=255)
                     gCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGAGAACCGTATACTGATg                                                                   >  2:5251836/1‑93 (MQ=255)
                       aCTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCg                                                                 >  1:4210514/1‑93 (MQ=255)
                       aCTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCg                                                                 >  1:5402201/1‑93 (MQ=255)
                        cTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGt                                                                >  1:5408420/1‑93 (MQ=255)
                           tATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTcc                                                             <  1:106171/93‑1 (MQ=255)
                            aTCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGGGAACCGTATACTGATGCGTTCcg                                                            >  2:551925/1‑93 (MQ=255)
                            aTCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCcg                                                            >  1:6250771/1‑93 (MQ=255)
                             tCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCcgc                                                           <  1:1597620/93‑1 (MQ=255)
                               tGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGa                                                         >  2:652864/1‑93 (MQ=255)
                               tGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGa                                                         >  1:3166836/1‑93 (MQ=255)
                                gATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGaa                                                        >  2:5725629/1‑93 (MQ=255)
                                 aTTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGaaa                                                       <  1:5902838/93‑1 (MQ=255)
                                  ttGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGaaaa                                                      >  2:5762014/1‑93 (MQ=255)
                                   tGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGaaaaa                                                     <  1:5836504/93‑1 (MQ=255)
                                   tGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGaaaaa                                                     <  2:2035209/93‑1 (MQ=255)
                                    gAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAg                                                    >  2:6017360/1‑93 (MQ=255)
                                         acTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTAcgc                                               >  1:5098039/1‑93 (MQ=255)
                                          cTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTAcgcg                                              <  2:1647712/93‑1 (MQ=255)
                                           tGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACgcgc                                             <  2:237920/93‑1 (MQ=255)
                                            gggTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGcc                                            <  2:3075089/93‑1 (MQ=255)
                                              gTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCg                                          <  2:5105075/93‑1 (MQ=255)
                                               tAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGt                                         <  1:5718592/93‑1 (MQ=255)
                                               tAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGt                                         <  1:3887991/93‑1 (MQ=255)
                                                aaaCGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGtt                                        >  1:4864167/1‑93 (MQ=255)
                                                 aaCGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGttt                                       >  1:96444/1‑93 (MQ=255)
                                                 aaCGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGttt                                       >  2:2892570/1‑93 (MQ=255)
                                                  aCGCTGCATATCTTATCTGCCTATGTCAGAACCATAATGGGCGGGGGTGAACCGTATACTGACGCGTTCCCCGAAAAAGTACGCGCCCGTTTc                                      <  2:4859943/93‑1 (MQ=255)
                                                  aCGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTc                                      >  2:3943371/1‑93 (MQ=255)
                                                   cgcgGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTcc                                     <  2:119371/93‑1 (MQ=255)
                                                   cgcgGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTcc                                     <  2:4227153/93‑1 (MQ=255)
                                                    gcgGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCa                                    >  2:3885705/1‑93 (MQ=255)
                                                    gcgGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCa                                    >  2:2625547/1‑93 (MQ=255)
                                                      ggCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACg                                  <  2:2759544/93‑1 (MQ=255)
                                                      ggCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACg                                  >  1:5597445/1‑93 (MQ=255)
                                                      ggCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACg                                  >  1:850261/1‑93 (MQ=255)
                                                      ggCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGCGAACCGCATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTTCACg                                  >  2:3728470/1‑93 (MQ=255)
                                                        cATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGt                                <  2:4790569/93‑1 (MQ=255)
                                                         aTTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTc                               <  1:3997922/93‑1 (MQ=255)
                                                          ttGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTcc                              >  2:5466014/1‑93 (MQ=255)
                                                            gCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGa                            >  2:6175622/1‑93 (MQ=255)
                                                             cTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGAt                           <  2:5137773/93‑1 (MQ=255)
                                                              ttATCTGCATATGTCAGAACCATATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGAtt                          <  1:6112537/93‑1 (MQ=255)
                                                              ttATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGAtt                          <  1:3963295/93‑1 (MQ=255)
                                                               tATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTa                         >  2:1307604/1‑93 (MQ=255)
                                                               tATCTGCATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTa                         >  2:5579516/1‑93 (MQ=255)
                                                                     cATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGgcg                   >  1:1370960/1‑93 (MQ=255)
                                                                     cATATGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGgcg                   <  2:1428946/93‑1 (MQ=255)
                                                                      atatGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGgcgc                  <  2:4529258/93‑1 (MQ=255)
                                                                      atatGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGgcgc                  >  1:1820358/1‑93 (MQ=255)
                                                                       tatGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCa                 >  1:4508176/1‑93 (MQ=255)
                                                                         tGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAgg               >  2:1076246/1‑93 (MQ=255)
                                                                         tGTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAgg               <  2:5283527/93‑1 (MQ=255)
                                                                          gTCAGAACCAGATTGGGCGGGGGGGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAggg              >  1:1944434/1‑92 (MQ=255)
                                                                          gTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAggg              >  1:257497/1‑92 (MQ=255)
                                                                          gTCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGt              >  2:2833426/1‑93 (MQ=255)
                                                                           tCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTg             >  1:742792/1‑93 (MQ=255)
                                                                           tCAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTg             >  1:2743314/1‑93 (MQ=255)
                                                                            cAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGc            <  1:640765/93‑1 (MQ=255)
                                                                            cAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGc            <  1:3200281/93‑1 (MQ=255)
                                                                            cAGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGc            <  1:3595458/93‑1 (MQ=255)
                                                                             aGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCa           >  2:2477057/1‑93 (MQ=255)
                                                                             aGAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCa           <  1:626027/93‑1 (MQ=255)
                                                                              gAACCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCaa          >  2:1838474/1‑92 (MQ=255)
                                                                               aaCCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATa         <  2:5798649/93‑1 (MQ=255)
                                                                                aCCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACTCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATac        >  2:4373695/1‑93 (MQ=255)
                                                                                aCCAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATac        >  1:2487995/1‑93 (MQ=255)
                                                                                  cAGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATacac      >  1:2393474/1‑93 (MQ=255)
                                                                                   aGATTGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATAcaca     >  2:5672170/1‑93 (MQ=255)
                                                                                   aGATTGGGCGGGGGCGAACCGTTTAGTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTTCACGGTCCGGTTATGGGCGTTGGTGCGTacagt     >  2:4502409/1‑91 (MQ=255)
                                                                                      ttGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTa  >  1:309138/1‑93 (MQ=255)
                                                                                      ttGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTa  <  2:4036828/93‑1 (MQ=255)
                                                                                      ttGGGCGGGGGTGAACCGTATACTGATGCGTTCCGCCAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTa  >  1:2137775/1‑93 (MQ=255)
                                                                                            |                                                                                       
GGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTA  >  NC_000913/1725421‑1725600

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: