Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,593,254 2 bp→GC coding (59‑60/807 nt) yjjM ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,593,2540CG100.0% 781.7 / NA 239A20G (GCG→GGCyjjMhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (134/105);  total (134/105)
*NC_0009134,593,2550GC99.2% 778.9 / ‑5.3 239A20G (GCG→GGC) yjjMhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (132/105);  minor base A (2/0);  total (134/105)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 244 total aligned reads displayed.

ACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACCGCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGAAG  >  NC_000913/4593167‑4593339
                                                                                       ||                                                                                    
aCTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCa                                                                                  >  2:1313142/1‑93 (MQ=255)
 cTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCAt                                                                                 >  1:2654705/1‑93 (MQ=255)
  tCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCAtt                                                                                >  1:3033614/1‑93 (MQ=255)
  tCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCAtt                                                                                <  1:2332985/93‑1 (MQ=255)
  tCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGACGCAtt                                                                                >  2:851532/1‑93 (MQ=255)
    aTTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTcg                                                                              >  1:1010807/1‑93 (MQ=255)
     ttGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTcgc                                                                             <  1:1181411/93‑1 (MQ=255)
      tGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTcgcg                                                                            <  1:676198/93‑1 (MQ=255)
       gAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGt                                                                           >  2:2836160/1‑93 (MQ=255)
        aaGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCATCCTGCGTCAAGACGCCGCATTCGCGTa                                                                          <  2:4019723/93‑1 (MQ=255)
         aGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTaa                                                                         <  1:1877169/93‑1 (MQ=255)
          gCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAg                                                                        <  1:4166043/93‑1 (MQ=255)
          gCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAg                                                                        >  1:2713376/1‑93 (MQ=255)
             gtagtaCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTgg                                                                     >  1:4261344/1‑93 (MQ=255)
               agtaCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCt                                                                   >  1:748213/1‑93 (MQ=255)
               agtaCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACCCCGCATTCGCGTAAGTGGcg                                                                   >  1:2475266/1‑92 (MQ=255)
                    aGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGaa                                                              <  2:45606/93‑1 (MQ=255)
                    aGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGaa                                                              <  1:1440969/93‑1 (MQ=255)
                     ggCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAat                                                             <  1:1213802/93‑1 (MQ=255)
                     ggCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAat                                                             >  2:2799081/1‑93 (MQ=255)
                      gctaaaCCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAata                                                            <  1:3301784/90‑1 (MQ=255)
                      gCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAata                                                            >  2:4246288/1‑93 (MQ=255)
                       cGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAatat                                                           >  2:3763923/1‑93 (MQ=255)
                        gAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATg                                                          <  2:1143086/93‑1 (MQ=255)
                         aaaCCGTCATCATGGTCAGGTTACGGCGAAAAAACAAAGTCGTTTCAAACTTCCGCCAAAACGCCGCATTGCCGTAAGGGGCTGAGAatatac                                                         >  2:1243488/1‑91 (MQ=255)
                          aaCCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGAc                                                        <  2:461304/93‑1 (MQ=255)
                          aaCCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGAc                                                        <  2:3623773/93‑1 (MQ=255)
                            ccGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGc                                                      <  2:724937/93‑1 (MQ=255)
                            ccGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGc                                                      <  1:1506249/93‑1 (MQ=255)
                              gTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCAc                                                    >  1:2865506/1‑93 (MQ=255)
                              ctcatcaTGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCAc                                                    <  1:1287638/92‑1 (MQ=255)
                               tcatcaTGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACt                                                   >  1:943893/1‑93 (MQ=255)
                                catcatGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACtg                                                  >  1:312333/1‑93 (MQ=255)
                                 atcatGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACtgt                                                 >  1:1834574/1‑93 (MQ=255)
                                 atcatGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACtgt                                                 <  1:217915/93‑1 (MQ=255)
                                  tcatGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGAGGCTGAGAATATGACGCACtcgg                                                >  2:3241041/1‑90 (MQ=255)
                                    atGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtggt                                              >  2:240422/1‑93 (MQ=255)
                                    atGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtggt                                              >  2:3396587/1‑93 (MQ=255)
                                    atGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtggt                                              >  2:1924743/1‑93 (MQ=255)
                                     tGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtggtg                                             >  1:1462611/1‑93 (MQ=255)
                                       gTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCg                                           <  1:228281/93‑1 (MQ=255)
                                         cAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGc                                         <  2:1675110/93‑1 (MQ=255)
                                          aGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCt                                        >  2:156188/1‑93 (MQ=255)
                                          aGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCt                                        <  2:4449852/93‑1 (MQ=255)
                                          aGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGCGCGGCt                                        <  2:4289506/93‑1 (MQ=255)
                                           ggTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTg                                       <  1:4560318/93‑1 (MQ=255)
                                           ggTGTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTg                                       <  2:1803598/93‑1 (MQ=255)
                                            gTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCGAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTgt                                      >  2:3672881/1‑92 (MQ=255)
                                            gTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGa                                      >  2:2403262/1‑93 (MQ=255)
                                             tttACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAt                                     <  2:267745/93‑1 (MQ=255)
                                              ttACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtg                                    <  1:819096/93‑1 (MQ=255)
                                                aCGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgtt                                  >  1:57189/1‑93 (MQ=255)
                                                  gcgcAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgttgt                                <  2:1034534/93‑1 (MQ=255)
                                                    gcAATAACATAGTCGTTTCCAACGTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTAc                              <  2:2627759/93‑1 (MQ=255)
                                                    gcAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTAc                              >  1:628016/1‑93 (MQ=255)
                                                    gcAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTAc                              >  2:3041559/1‑93 (MQ=255)
                                                     cAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACa                             <  2:370138/93‑1 (MQ=255)
                                                     cAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACa                             >  1:3078107/1‑93 (MQ=255)
                                                      aaTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAt                            >  2:1795668/1‑93 (MQ=255)
                                                      aaTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAt                            >  1:2876429/1‑93 (MQ=255)
                                                       aTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAtt                           <  1:355699/93‑1 (MQ=255)
                                                        tAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCGGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACActt                          >  1:1928049/1‑93 (MQ=255)
                                                         aaCATAGTCTGTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTc                         <  1:1170595/93‑1 (MQ=255)
                                                         aaCATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTc                         <  2:3489651/93‑1 (MQ=255)
                                                          aCATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTcc                        >  2:1065245/1‑93 (MQ=255)
                                                          aCATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTcc                        <  2:3905557/93‑1 (MQ=255)
                                                            aTAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGc                      >  2:3319720/1‑93 (MQ=255)
                                                            aTAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGc                      >  2:2170348/1‑93 (MQ=255)
                                                            aTAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGc                      >  1:3708831/1‑93 (MQ=255)
                                                             tAGTCGTTTCCACCCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGcc                     <  2:3510606/93‑1 (MQ=255)
                                                              aGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCa                    <  2:3049417/93‑1 (MQ=255)
                                                              aGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCa                    <  2:4150765/93‑1 (MQ=255)
                                                              aGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCa                    <  1:2059500/93‑1 (MQ=255)
                                                               gTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAg                   <  2:2733845/93‑1 (MQ=255)
                                                                tCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTCCATTTCCGCCAgc                  >  2:1335758/1‑93 (MQ=255)
                                                                tCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAgc                  >  1:3567964/1‑93 (MQ=255)
                                                                tCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAgc                  >  2:2925051/1‑93 (MQ=255)
                                                                tCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAgc                  >  2:1447632/1‑93 (MQ=255)
                                                                 cGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAgcg                 >  1:3756190/1‑93 (MQ=255)
                                                                  gTTTCCTACCTGCGTCAAGACGCCGCATTCGCGTTAGTGGCAGAGAATACGACGCACTGTGGTGCGGCTGTTGGTGTACATTTCCGCCAgcgc                >  2:4333752/1‑93 (MQ=255)
                                                                   tttCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTGCGCCAGCGCa               >  1:1143873/1‑93 (MQ=255)
                                                                   tttCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCCGCGCa               >  2:336679/1‑93 (MQ=255)
                                                                   tttCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCa               >  1:1715461/1‑93 (MQ=255)
                                                                     tCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACt             <  2:4427115/93‑1 (MQ=255)
                                                                     tCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACt             <  1:3927602/93‑1 (MQ=255)
                                                                      ccAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACtt            >  2:2295515/1‑93 (MQ=255)
                                                                       cAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACttt           >  2:2554556/1‑93 (MQ=255)
                                                                       cAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACttt           <  1:926492/93‑1 (MQ=255)
                                                                       cAACCGGCGGCAAGACGCCGCATTCTCGTAAGTGGCGGAGAATATGACGCACTGTGTTGCGGCTGATGTTGTACATTTCCGCCAGCGCAcgtg           >  2:3579907/1‑90 (MQ=255)
                                                                        aaCCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTg          >  1:843043/1‑93 (MQ=255)
                                                                         aCCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGc         >  1:1264618/1‑93 (MQ=255)
                                                                           cTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGa       <  1:3005095/93‑1 (MQ=255)
                                                                            tGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGAt      >  2:1214733/1‑93 (MQ=255)
                                                                             gCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATg     <  2:1949699/93‑1 (MQ=255)
                                                                             gCGTCAAGACGCCGCATTCGCGTAAGAGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTCCATTTCCGCCAGCGCACTTTGCGATg     >  2:4080871/1‑93 (MQ=255)
                                                                               gTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgga   <  2:3191162/93‑3 (MQ=255)
                                                                               gTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgga   <  1:3843321/93‑3 (MQ=255)
                                                                                tCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATggaa  >  1:4398539/1‑90 (MQ=255)
                                                                                tCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATggaa  >  1:4732412/1‑90 (MQ=255)
                                                                                tCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATggaa  <  1:504959/93‑4 (MQ=255)
                                                                                       ||                                                                                    
ACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACCGCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGAAG  >  NC_000913/4593167‑4593339

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