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% 2006.9 / NA 514A20G (GCG→GGCyjjMhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (253/261);  total (253/261)
*NC_0009134,593,2550GC99.8% 1977.2 / NA 514A20G (GCG→GGC) yjjMhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (251/261);  total (253/261)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 525 total aligned reads displayed.

TTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACCGCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGAAGTGGG  >  NC_000913/4593164‑4593343
                                                                                          ||                                                                                         
ttCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCa                                                                                       >  1:14107883/1‑96 (MQ=255)
 tCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCAt                                                                                      <  1:8126517/96‑1 (MQ=255)
 tCACTCATTGAAGCGGTAGTACAGGAGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCAt                                                                                      >  2:6290106/1‑96 (MQ=255)
  cACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCAtt                                                                                     <  1:6815584/96‑1 (MQ=255)
   aCTCATTGACGCGGTAGTACAGGCGAAACCGTCACCATGGTCAGGTTTACGCGCAATAACATAGTCCTTTCCAACCTGCGTCAAGACGCCGCATTc                                                                                    <  1:3255904/96‑1 (MQ=255)
   aCTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTc                                                                                    >  2:600698/1‑96 (MQ=255)
   aCTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTc                                                                                    <  2:4749947/96‑1 (MQ=255)
   aCTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTc                                                                                    <  2:12803396/96‑1 (MQ=255)
    cTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTcg                                                                                   <  1:3759355/96‑1 (MQ=255)
     tCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTcgc                                                                                  <  1:2135829/96‑1 (MQ=255)
       aTTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGt                                                                                <  1:3675005/96‑1 (MQ=255)
         tGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTaa                                                                              >  1:5935302/1‑96 (MQ=255)
         tGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTaa                                                                              <  1:5126492/96‑1 (MQ=255)
           aaGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGt                                                                            <  1:11632346/96‑1 (MQ=255)
             gCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTgg                                                                          <  2:5289939/96‑1 (MQ=255)
              cGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGc                                                                         >  1:7464444/1‑96 (MQ=255)
               ggTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCt                                                                        <  1:6201390/96‑1 (MQ=255)
                gtagtaCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTg                                                                       <  1:11003413/96‑1 (MQ=255)
                 tagtaCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTga                                                                      >  2:1350551/1‑96 (MQ=255)
                  agtaCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTgag                                                                     >  1:6643838/1‑96 (MQ=255)
                  agtaCAGGAGAAACCGTCATCATGGGCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTATGATGCTgag                                                                     >  1:5798087/1‑96 (MQ=255)
                   gtaCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTgaga                                                                    <  2:9681488/96‑1 (MQ=255)
                   gtaCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGGAAGTGGCTgaga                                                                    >  1:2428637/1‑96 (MQ=255)
                     aCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAat                                                                  <  2:8089876/96‑1 (MQ=255)
                      cAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAata                                                                 <  2:5488434/96‑1 (MQ=255)
                       aGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAatat                                                                >  1:2612582/1‑96 (MQ=255)
                       aGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAatat                                                                >  2:2780573/1‑96 (MQ=255)
                       aGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAatat                                                                >  2:5002650/1‑96 (MQ=255)
                       aGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAatag                                                                >  1:5311344/1‑95 (MQ=255)
                       aGGAGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAatat                                                                >  1:4691430/1‑96 (MQ=255)
                        ggCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATg                                                               >  2:6612636/1‑96 (MQ=255)
                          cGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGAc                                                             >  1:8632349/1‑96 (MQ=255)
                          cGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACAGAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGAc                                                             >  2:12615730/1‑96 (MQ=255)
                           gAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACg                                                            >  2:4192291/1‑96 (MQ=255)
                           gAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACg                                                            >  1:8856483/1‑96 (MQ=255)
                            aaaCCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGc                                                           >  2:8655774/1‑96 (MQ=255)
                             aaCCGTCATCATGCTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCa                                                          >  1:13441672/1‑96 (MQ=255)
                              aCCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGAGGCGGCATTCGCGTAAGTGGCTGAGAATATGACGCAc                                                         >  1:11548336/1‑96 (MQ=255)
                               ccGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACt                                                        <  2:380770/96‑1 (MQ=255)
                                cGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACtg                                                       <  2:1865421/96‑1 (MQ=255)
                                cGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACtg                                                       >  2:11623134/1‑96 (MQ=255)
                                 gTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACtgt                                                      >  1:634316/1‑96 (MQ=255)
                                 gTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACtgt                                                      >  1:9150893/1‑96 (MQ=255)
                                     tcatGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtggtg                                                  <  1:1308873/96‑1 (MQ=255)
                                      catGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGc                                                 <  2:11460131/96‑1 (MQ=255)
                                      catGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGc                                                 <  1:9740064/96‑1 (MQ=255)
                                       atGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCg                                                >  2:1008196/1‑96 (MQ=255)
                                        tGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCgg                                               <  1:7186535/96‑1 (MQ=255)
                                        tGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCAGCATTCGCGTAAGAGGCTGAGAATATTGCGCACTGTGGTGCgg                                               >  1:13938753/1‑96 (MQ=255)
                                         ggTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGc                                              >  2:7430824/1‑96 (MQ=255)
                                           tCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTg                                            >  1:5542894/1‑96 (MQ=255)
                                            cAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGa                                           <  2:4356266/96‑1 (MQ=255)
                                             aGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAt                                          >  1:10007008/1‑96 (MQ=255)
                                             aGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAt                                          <  2:10532733/96‑1 (MQ=255)
                                              ggTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtg                                         >  2:9006586/1‑96 (MQ=255)
                                                 ttACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgttg                                      <  1:11884648/96‑1 (MQ=255)
                                                 ttACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgttg                                      >  2:7629253/1‑96 (MQ=255)
                                                 ttACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgttg                                      >  2:7244936/1‑96 (MQ=255)
                                                 ttACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgttg                                      >  2:14336233/1‑96 (MQ=255)
                                                 ttACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgtag                                      >  1:10442882/1‑94 (MQ=255)
                                                  tccgcgCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgttgt                                     <  2:2247976/94‑1 (MQ=255)
                                                     gcgcAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACa                                  <  2:12919315/96‑1 (MQ=255)
                                                       gcAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAAAATGAGGACTTGGGGGGCGGCGGATGTTGTAAAtt                                >  2:315715/1‑96 (MQ=255)
                                                        cAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAttt                               <  2:11873809/96‑1 (MQ=255)
                                                         aaTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTTATGTTGTACATTTc                              >  1:14417799/1‑96 (MQ=255)
                                                         aaTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAttgc                              <  1:10785110/96‑3 (MQ=255)
                                                         aaTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTc                              >  1:12195307/1‑96 (MQ=255)
                                                          aTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCCCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTcc                             >  2:13266814/1‑96 (MQ=255)
                                                           tAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGATTATGACGCACTGTGGTGCGGATGATGTTGTACATTTCCg                            <  1:12615730/96‑1 (MQ=255)
                                                            aaCATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGc                           >  1:4937025/1‑96 (MQ=255)
                                                             aCATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGcc                          >  1:9946345/1‑96 (MQ=255)
                                                              cATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCa                         >  2:1211932/1‑96 (MQ=255)
                                                              cATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCa                         <  2:8849827/96‑1 (MQ=255)
                                                               aTAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAg                        >  1:7030633/1‑96 (MQ=255)
                                                                tAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAgc                       <  2:11224567/96‑1 (MQ=255)
                                                                tAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAgc                       <  2:8244299/96‑1 (MQ=255)
                                                                  gTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAgcgc                     <  2:10898832/96‑1 (MQ=255)
                                                                  gTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAgcgc                     <  1:3435921/96‑1 (MQ=255)
                                                                  gTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAgcgc                     >  2:4662240/1‑96 (MQ=255)
                                                                     gTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGGATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACt                  >  2:6867054/1‑96 (MQ=255)
                                                                     gTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACt                  >  2:9181620/1‑96 (MQ=255)
                                                                      tttCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACtt                 >  2:2083623/1‑96 (MQ=255)
                                                                        tCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTg               <  1:9486505/96‑1 (MQ=255)
                                                                         ccAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGc              >  1:2788997/1‑96 (MQ=255)
                                                                           aaCCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGa            <  1:4103867/96‑1 (MQ=255)
                                                                             ccTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGGACTTTGCGATg         >  1:7808396/1‑96 (MQ=255)
                                                                             ccTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATg         <  2:12274647/96‑1 (MQ=255)
                                                                              cTGNGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgg        >  2:5672318/1‑95 (MQ=255)
                                                                              cTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgg        <  2:6360923/96‑2 (MQ=255)
                                                                              cTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgg        >  2:346026/1‑95 (MQ=255)
                                                                               tGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgga       >  1:1915302/1‑94 (MQ=255)
                                                                                gCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATggaa      <  1:8338167/96‑4 (MQ=255)
                                                                                 cGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATggaag     <  2:13784202/96‑5 (MQ=255)
                                                                                   tCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTg    <  2:13999400/96‑1 (MQ=255)
                                                                                   tCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTg    <  1:13619345/96‑1 (MQ=255)
                                                                                   tCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTg    <  1:4503225/96‑1 (MQ=255)
                                                                                    cAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTCGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTgg   >  2:10393294/1‑96 (MQ=255)
                                                                                     aaGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTggg  >  1:4264411/1‑96 (MQ=255)
                                                                                     aaGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTggg  <  2:9398209/96‑1 (MQ=255)
                                                                                     aaGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTggg  >  2:2559150/1‑96 (MQ=255)
                                                                                          ||                                                                                         
TTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACCGCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGAAGTGGG  >  NC_000913/4593164‑4593343

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: