Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,952,380 G→A 35.6% L60L (CTG→CTA ilvM → acetolactate synthase 2 small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,952,3800GA35.6% 105.4 / 80.3 95L60L (CTG→CTAilvMacetolactate synthase 2 small subunit
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base A (17/17);  total (47/48)
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

GTTTCCACGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAGGAATA  >  NC_000913/3952283‑3952478
                                                                                                 |                                                                                                  
gTTTCCACGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTAtt                                                                                                  >  1:811794/1‑100 (MQ=255)
 tttCCACGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGttt                                                                                                 <  1:839258/100‑1 (MQ=255)
 tttCCACGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGttt                                                                                                 <  1:1954926/100‑1 (MQ=255)
  ttCCACGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTa                                                                                                >  1:1823652/1‑100 (MQ=255)
    ccACGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGt                                                                                              <  1:2074285/100‑1 (MQ=255)
      aCGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCa                                                                                            >  1:1412379/1‑100 (MQ=255)
        gTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGt                                                                                          >  2:1373987/1‑100 (MQ=255)
         tCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGtt                                                                                         >  1:214296/1‑100 (MQ=255)
          cTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTta                                                                                        <  1:1675260/100‑1 (MQ=255)
          cTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTta                                                                                        <  2:728536/100‑1 (MQ=255)
           tGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTtaa                                                                                       <  1:2212334/100‑1 (MQ=255)
           tGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTtaa                                                                                       >  1:734561/1‑100 (MQ=255)
           tGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTtaa                                                                                       <  1:1156659/100‑1 (MQ=255)
            gCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTtaaa                                                                                      >  1:1556769/1‑100 (MQ=255)
            gCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTtaaa                                                                                      >  1:1592786/1‑100 (MQ=255)
            gCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTtaaa                                                                                      <  2:2040126/100‑1 (MQ=255)
             cTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTtaaat                                                                                     <  1:1049308/100‑1 (MQ=255)
             cTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTtaaa                                                                                      <  2:1485840/99‑1 (MQ=255)
             cTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTtaaat                                                                                     >  2:160745/1‑100 (MQ=255)
                 aTGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAAc                                                                                 <  2:286686/100‑1 (MQ=255)
                  tGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACt                                                                                <  1:1553391/100‑1 (MQ=255)
                      tatGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGTCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACtggtg                                                                            <  2:369082/100‑1 (MQ=255)
                      tatGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACttgtg                                                                            >  1:1425811/1‑100 (MQ=255)
                      tatGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACtggtg                                                                            <  1:677480/100‑1 (MQ=255)
                      tatGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACtggtg                                                                            >  2:337518/1‑100 (MQ=255)
                       atGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACtggtgg                                                                           <  1:1621384/100‑1 (MQ=255)
                       atGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACtggtgg                                                                           <  1:726015/100‑1 (MQ=255)
                        tGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGa                                                                          <  1:1041736/100‑1 (MQ=255)
                        tGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGa                                                                          >  1:2157887/1‑100 (MQ=255)
                         ggCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGAc                                                                         >  1:364555/1‑100 (MQ=255)
                          gccgccgccAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACg                                                                        >  1:1414689/1‑100 (MQ=255)
                          gccgccgccAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACg                                                                        >  1:516282/1‑100 (MQ=255)
                            cgccgccAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTc                                                                      >  2:849375/1‑100 (MQ=255)
                             gccgccAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCg                                                                     <  2:1110877/100‑1 (MQ=255)
                              ccgccAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGc                                                                    <  1:132347/100‑1 (MQ=255)
                              ccgccAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGc                                                                    >  1:429556/1‑100 (MQ=255)
                               cgccAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGca                                                                   <  1:1546908/100‑1 (MQ=255)
                                      gATGCACAAAATATAAATATCGAATTGAGCGGTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTg                                                            <  2:1821850/100‑1 (MQ=255)
                                      gATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTg                                                            <  1:1032484/100‑1 (MQ=255)
                                      gATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTg                                                            <  2:964186/100‑1 (MQ=255)
                                       aTGCACAAAATATAAATATCGATTTTACCGTTGCTAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGc                                                           <  2:207011/100‑1 (MQ=255)
                                       aTGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGc                                                           <  1:2126449/100‑1 (MQ=255)
                                         gCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCa                                                         >  1:1747755/1‑100 (MQ=255)
                                          cacaAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCAt                                                        >  2:1861618/1‑100 (MQ=255)
                                          cacaAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCAt                                                        >  1:1760257/1‑100 (MQ=255)
                                           acaAAATATAAATATCGAATTCACCGTTGCCAGCCCACGGTCGGTCGACTTACTGCTTACACAGTTAAATAAACCGGTGGACGTCGCACACGTTGCCATc                                                       >  2:196854/1‑100 (MQ=255)
                                            caAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCt                                                      >  1:1489534/1‑100 (MQ=255)
                                             aaaaTATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTg                                                     <  1:1861630/100‑1 (MQ=255)
                                              aaatatatatatCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTg                                                     <  1:149135/99‑1 (MQ=255)
                                                     aatatCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGca                                             >  2:1098920/1‑100 (MQ=255)
                                                     aatatCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGca                                             <  2:2065082/100‑1 (MQ=255)
                                                       tatCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGcaca                                           >  2:764647/1‑100 (MQ=255)
                                                         tCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAAc                                         >  1:2199446/1‑100 (MQ=255)
                                                         tCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAAc                                         >  1:1551960/1‑100 (MQ=255)
                                                         tCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAAc                                         >  1:1557038/1‑100 (MQ=255)
                                                          cGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAAcc                                        <  1:1513875/100‑1 (MQ=255)
                                                          cGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAAcc                                        <  1:1012197/100‑1 (MQ=255)
                                                          cGAATGGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAAcc                                        <  2:193705/100‑1 (MQ=255)
                                                           gAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACca                                       >  1:2245195/1‑100 (MQ=255)
                                                            aaTTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACcac                                      <  1:1014069/100‑1 (MQ=255)
                                                               tGACCGTTGCCAGCCCATGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCTATCTGCCAGAGCACAACCACATc                                   <  2:161686/100‑1 (MQ=255)
                                                                 aCCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTACATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATcac                                 <  1:535193/100‑1 (MQ=255)
                                                                  ccGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATcaca                                >  1:1263848/1‑100 (MQ=255)
                                                                  ccGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATcaca                                >  1:65633/1‑100 (MQ=255)
                                                                  ccGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATcaca                                >  2:134315/1‑100 (MQ=255)
                                                                  ccGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATcaca                                >  2:1762306/1‑100 (MQ=255)
                                                                  ccGTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATcac                                 >  1:2021382/1‑99 (MQ=255)
                                                                  ccGTTGCCAGCACACGGTTGGTCGATTTACTATTTAGTCAGTTAAATATTCTGGTTGACCTCTCTCACGTTGCGATCTGCCAGAGCACAACCACATcac                                 <  2:550244/99‑1 (MQ=255)
                                                                   cGTTGCCAGCCCACGGTCGGACGACTTAGTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACCCGTTGCCATCTGCCAGAGGACAACTACAGcaca                                >  2:1104899/1‑99 (MQ=255)
                                                                    gTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGTCATCTGCCAGAGCACAACCACATCacaac                              <  2:1957834/100‑1 (MQ=255)
                                                                    gTTGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCacaac                              >  2:667411/1‑100 (MQ=255)
                                                                    gTTGCCAGCCCACGGTCGGTCGACTTACTATGTAGTCAGTTAAATAAACTGGTGTACGTCGCACACGTTGTTATCTGCGAGATCACAACCACATCacaac                              <  2:1957787/100‑1 (MQ=255)
                                                                     ttGCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTTAATAGACTGGTGGACGTCGCACACGTTGCGATCTGCCAGAGCACAACCACAGCacaaca                             <  2:434559/100‑1 (MQ=255)
                                                                       gCCAGCCCACGGTCGGTCGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACaaa                           >  1:1680618/1‑100 (MQ=255)
                                                                             ccACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCcgcg                     <  1:1067476/100‑1 (MQ=255)
                                                                             ccACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCcgcg                     >  1:668775/1‑100 (MQ=255)
                                                                             ccACGGTCGGTCGACGTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCcgcg                     <  2:1926974/100‑1 (MQ=255)
                                                                                 ggtcggtcGACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTg                 >  1:575472/1‑100 (MQ=255)
                                                                                   tcggtcgACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAg               <  1:1112444/100‑1 (MQ=255)
                                                                                     ggtcgACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAgcg             >  2:2052341/1‑100 (MQ=255)
                                                                                     ggtcgACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAgcg             <  1:1075637/100‑1 (MQ=255)
                                                                                     ggtcgACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAgcg             >  2:1637837/1‑100 (MQ=255)
                                                                                     ggtcgACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAgcg             >  1:1564207/1‑100 (MQ=255)
                                                                                       tcgACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCa           >  2:1777383/1‑100 (MQ=255)
                                                                                       tcgACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCa           <  1:1331332/100‑1 (MQ=255)
                                                                                        cgACTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGAACAACCACATCACAACAAATCCGCGCCTGAGCGCa           >  2:926678/1‑99 (MQ=255)
                                                                                           cTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAg       >  1:550502/1‑100 (MQ=255)
                                                                                           cTTACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAg       <  2:905467/100‑1 (MQ=255)
                                                                                            ttACTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAg       <  2:1530726/99‑1 (MQ=255)
                                                                                             caCTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAGGa     >  1:1676818/2‑100 (MQ=255)
                                                                                              aCTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAGGaa    <  1:764647/100‑1 (MQ=255)
                                                                                              aCTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAGGaa    >  2:109801/1‑100 (MQ=255)
                                                                                              aCTATTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAGGaa    >  2:1166096/1‑100 (MQ=255)
                                                                                               cTGTTTAGTCAGTTAAATAAACTGGTGGTCGTCGCACACCTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAGGaa    <  2:1823640/99‑1 (MQ=255)
                                                                                                tGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAGGAata  <  1:1095058/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
GTTTCCACGTCTGCTCAATGAATATGGCCGCCGCCAGCGATGCACAAAATATAAATATCGAATTGACCGTTGCCAGCCCACGGTCGGTCGACTTACTGTTTAGTCAGTTAAATAAACTGGTGGACGTCGCACACGTTGCCATCTGCCAGAGCACAACCACATCACAACAAATCCGCGCCTGAGCGCAAAAGGAATA  >  NC_000913/3952283‑3952478

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