Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,001,486 2 bp→GC coding (316‑317/903 nt) rarD ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,001,4860CG99.6% 1519.6 / NA 450R106A (CGG→GCG) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base G (239/209);  total (241/209)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,001,4870GC99.8% 1528.7 / NA 450R106A (CGG→GCG) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (240/209);  total (241/209)
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.00e+00

Only 100 of 451 total aligned reads displayed.

AGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTC  >  NC_000913/4001399‑4001578
                                                                                       ||                                                                                           
agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                                                                                         >  1:381887/1‑93 (MQ=255)
 gccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCa                                                                                        >  1:1001377/1‑93 (MQ=255)
 gccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCa                                                                                        >  2:2836163/1‑93 (MQ=255)
  ccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAg                                                                                       >  1:3996912/1‑93 (MQ=255)
  ccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGCAACCAAGGCTCGCTTCCAg                                                                                       <  1:6552151/93‑1 (MQ=255)
   caTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGc                                                                                      >  2:6813260/1‑93 (MQ=255)
   caTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGc                                                                                      >  2:5392030/1‑93 (MQ=255)
   caTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGc                                                                                      >  1:5866591/1‑93 (MQ=255)
    aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa                                                                                     <  1:5171042/93‑1 (MQ=255)
    aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa                                                                                     <  2:272203/93‑1 (MQ=255)
    aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa                                                                                     <  1:3089218/93‑1 (MQ=255)
     ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat                                                                                    >  1:3288611/1‑93 (MQ=255)
     ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat                                                                                    >  2:3821858/1‑93 (MQ=255)
     ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat                                                                                    >  1:4199375/1‑93 (MQ=255)
      tGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCata                                                                                   <  2:6676579/93‑1 (MQ=255)
      tGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGCAACCAAGGCTCGCTTCCAGCata                                                                                   <  1:924654/93‑1 (MQ=255)
       gCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCatat                                                                                  <  1:4353448/93‑1 (MQ=255)
       gCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGGTAATAAAGTAACCAAGGCTCGCTTCCAGCatat                                                                                  >  2:6013406/1‑93 (MQ=255)
        cATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATg                                                                                 <  2:2048841/93‑1 (MQ=255)
          tCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGt                                                                               >  1:7382675/1‑93 (MQ=255)
             ggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAt                                                                            <  1:3610775/93‑1 (MQ=255)
              gcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAtt                                                                           <  1:2972617/93‑1 (MQ=255)
                ggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGt                                                                         >  1:3404586/1‑93 (MQ=255)
                     cgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAcc                                                                    <  1:7477634/90‑1 (MQ=255)
                     cgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAcc                                                                    <  2:6116786/90‑1 (MQ=255)
                       ctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGc                                                                  >  2:243678/2‑93 (MQ=255)
                       ctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGc                                                                  >  1:5670669/2‑93 (MQ=255)
                        tCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGcc                                                                 <  2:1003476/93‑1 (MQ=255)
                        tCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGcc                                                                 >  2:5889843/1‑93 (MQ=255)
                        tCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGcc                                                                 >  2:4313186/1‑93 (MQ=255)
                           ccGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGCAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAg                                                              <  1:5829797/93‑1 (MQ=255)
                            cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGa                                                             >  1:6230365/1‑93 (MQ=255)
                            cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGa                                                             <  2:1918526/93‑1 (MQ=255)
                            cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGa                                                             >  2:2621651/1‑93 (MQ=255)
                            cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGa                                                             <  1:1946921/93‑1 (MQ=255)
                            cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGCCTCCCTTCCACCATATGGTGATTGTTCACCGCCCAGa                                                             >  2:1135048/1‑93 (MQ=255)
                             gAGGAAAATCATCCCCAGCACAATGTTCACCAACGGGGTAATAAAGCAACCAAGGGTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGAt                                                            >  1:2696661/1‑93 (MQ=255)
                               ggAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCGGCATATGGTGATTGTTCACCGCCCAGATaa                                                          >  2:7164226/1‑93 (MQ=255)
                               ggAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaa                                                          >  1:4726376/1‑93 (MQ=255)
                               ggAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaa                                                          >  2:5476394/1‑93 (MQ=255)
                               ggAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaa                                                          <  1:2163102/93‑1 (MQ=255)
                                gAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaaa                                                         >  2:524416/1‑93 (MQ=255)
                                  aaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACa                                                       <  1:819313/93‑1 (MQ=255)
                                   aaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAg                                                      <  1:1678476/93‑1 (MQ=255)
                                   aaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAg                                                      >  1:3030641/1‑93 (MQ=255)
                                   aaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGCAACCAAGGCTCGGTTCCCGCATATGGTGATTGTTGACCGCCCAGATAAAAAg                                                      >  1:1592214/1‑93 (MQ=255)
                                        tcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGcc                                                 >  1:2068808/1‑93 (MQ=255)
                                        tcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGcc                                                 >  2:1730325/1‑93 (MQ=255)
                                        tcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGCAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGcc                                                 >  2:2877225/1‑93 (MQ=255)
                                          cccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAg                                               <  1:5315618/93‑1 (MQ=255)
                                           ccAGCACAATGTTCACCAGCGGGTCAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGt                                              >  1:4989761/1‑93 (MQ=255)
                                            cAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGtt                                             <  1:6448988/93‑1 (MQ=255)
                                            cAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGtt                                             >  1:1431581/1‑93 (MQ=255)
                                                acaATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGcca                                         >  2:401943/1‑93 (MQ=255)
                                                acaATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGcca                                         >  2:305376/1‑93 (MQ=255)
                                                acaATGTTCACCAGCGGGTTAATAAAGCAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGcca                                         >  2:3163655/1‑93 (MQ=255)
                                                 caATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccac                                        >  2:2384472/1‑93 (MQ=255)
                                                  aaTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacc                                       >  2:5159943/1‑93 (MQ=255)
                                                  aaTGTTCACCAGCGGGTTAATAAAGCAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacc                                       >  1:287070/1‑93 (MQ=255)
                                                   aTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacca                                      >  2:743737/1‑93 (MQ=255)
                                                   aTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacca                                      >  1:7148656/1‑93 (MQ=255)
                                                   aTGTTCACCAGCGGGTTAATAAAGCAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacca                                      <  1:5544859/93‑1 (MQ=255)
                                                    tGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCaa                                     >  2:3460256/1‑93 (MQ=255)
                                                    tGTTCACCAGCGGGTTAATAAAGCAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCaa                                     >  2:128751/1‑93 (MQ=255)
                                                     gTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAAt                                    <  2:6990326/93‑1 (MQ=255)
                                                      ttCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATc                                   >  1:6608700/1‑93 (MQ=255)
                                                       tCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCa                                  <  1:6300070/93‑1 (MQ=255)
                                                       tCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCa                                  <  2:3754860/93‑1 (MQ=255)
                                                       tCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCa                                  <  2:7499406/93‑1 (MQ=255)
                                                         aCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGc                                <  1:4600468/93‑1 (MQ=255)
                                                         aCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGc                                >  1:7269351/1‑93 (MQ=255)
                                                         aCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCAAGTTGCCACCAATCAGc                                >  1:4426597/1‑93 (MQ=255)
                                                          ccAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCa                               <  1:673085/93‑1 (MQ=255)
                                                          ccAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCa                               <  2:6232677/93‑1 (MQ=255)
                                                          ccAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCa                               <  2:6485636/93‑1 (MQ=255)
                                                            aGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACg                             >  1:1137409/1‑93 (MQ=255)
                                                            aGCGGGTTAATAAAGAAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACAAATCAGCACg                             >  1:7015275/1‑93 (MQ=255)
                                                             gCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACgg                            <  1:6083948/93‑1 (MQ=255)
                                                              cGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGc                           <  2:3906659/93‑1 (MQ=255)
                                                              cGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGc                           >  1:317072/1‑93 (MQ=255)
                                                               gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa                          <  2:36861/93‑1 (MQ=255)
                                                               gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa                          <  2:1575040/93‑1 (MQ=255)
                                                               gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa                          <  1:2361651/93‑1 (MQ=255)
                                                                ggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCag                         >  2:4860414/1‑93 (MQ=255)
                                                                ggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCag                         >  2:2191573/1‑93 (MQ=255)
                                                                ggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCag                         >  2:6308420/1‑93 (MQ=255)
                                                                 gTTAATAAAGCAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCaga                        <  2:790933/93‑1 (MQ=255)
                                                                  ttAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCagag                       >  1:6765362/1‑93 (MQ=255)
                                                                  ttAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCagag                       <  2:4799878/93‑1 (MQ=255)
                                                                  ttAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCagag                       <  1:5974980/93‑1 (MQ=255)
                                                                    aataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCGCCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGAc                     <  2:7317225/93‑1 (MQ=255)
                                                                     ataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACt                    >  1:5430027/1‑93 (MQ=255)
                                                                       aaaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGc                  <  1:1251804/93‑1 (MQ=255)
                                                                        aaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGcc                 >  1:4802962/1‑93 (MQ=255)
                                                                           tAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGc              >  2:1437912/1‑93 (MQ=255)
                                                                            aaCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCa             >  2:656014/1‑93 (MQ=255)
                                                                            aaCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCa             <  2:4222494/93‑1 (MQ=255)
                                                                            aaCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCa             >  2:3327010/1‑93 (MQ=255)
                                                                              ccAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATa           <  1:3949350/93‑1 (MQ=255)
                                                                               cAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATaa          >  2:451741/1‑93 (MQ=255)
                                                                                       ||                                                                                           
AGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTC  >  NC_000913/4001399‑4001578

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