Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,334,354 T→C I36V (ATT→GTT)  yrbA ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,334,3540TC97.4% 598.5 / 3.1 194I36V (ATT→GTT) yrbAhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (5/0);  new base C (106/83);  total (111/83)
Fisher's exact test for biased strand distribution p-value = 7.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

Only 100 of 203 total aligned reads displayed.

GGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCC  >  NC_000913/3334264‑3334442
                                                                                          |                                                                                        
 gTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCaacaac                                                                                       <  2:340715/93‑1 (MQ=255)
 gTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCaacaac                                                                                       <  2:3516386/93‑1 (MQ=255)
  ttATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAAcc                                                                                      >  1:617476/1‑93 (MQ=255)
  ttATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAAcc                                                                                      >  1:1796108/1‑93 (MQ=255)
   tATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCt                                                                                     >  2:3928172/1‑93 (MQ=255)
   tATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCt                                                                                     <  2:3541286/93‑1 (MQ=255)
   tATCCCCAATATATTCCTTCAGCGGACGATAGACCGTCTGCTTTTTTTTAACCCCGCTCATGCCGTCAAACAACTCACCCGCGGCAACAACCt                                                                                     <  2:1605269/93‑1 (MQ=255)
    aTCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTg                                                                                    >  2:2935908/1‑93 (MQ=255)
     tCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGa                                                                                   <  1:989374/93‑1 (MQ=255)
     tCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGa                                                                                   >  2:1308442/1‑93 (MQ=255)
       cGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGaaa                                                                                 >  2:1874394/1‑93 (MQ=255)
        gCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAg                                                                                >  1:2110961/1‑93 (MQ=255)
        gCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAg                                                                                >  1:3473581/1‑93 (MQ=255)
         cAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGt                                                                               <  2:4287925/93‑1 (MQ=255)
         cAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGt                                                                               <  1:2330702/93‑1 (MQ=255)
            tatatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGc                                                                            <  2:3928882/93‑1 (MQ=255)
            tatatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGc                                                                            <  1:4694249/93‑1 (MQ=255)
            tatatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGc                                                                            >  1:709776/1‑93 (MQ=255)
              tatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTg                                                                          >  2:3739093/1‑93 (MQ=255)
              tatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTg                                                                          >  2:4581006/1‑93 (MQ=255)
               atTCCATCAGCGGACCATAGACCGTCTGTTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGc                                                                         <  1:3443659/93‑1 (MQ=255)
               atTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGc                                                                         <  2:2189153/93‑1 (MQ=255)
               atTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGc                                                                         <  1:3419958/93‑1 (MQ=255)
                ttCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGcc                                                                        >  2:1261363/1‑93 (MQ=255)
                ttCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGcc                                                                        >  1:4238354/1‑93 (MQ=255)
                 tCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCa                                                                       >  2:1885037/1‑93 (MQ=255)
                  ccATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCAt                                                                      <  1:4528955/93‑1 (MQ=255)
                   cATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATc                                                                     >  2:1073413/1‑93 (MQ=255)
                   cATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATc                                                                     <  1:20407/93‑1 (MQ=255)
                    aTCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGGCATCg                                                                    >  2:4317735/1‑93 (MQ=255)
                    aTCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCg                                                                    >  2:866873/1‑93 (MQ=255)
                     tCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGc                                                                   >  2:1745541/1‑93 (MQ=255)
                      cAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGcc                                                                  <  1:320911/93‑1 (MQ=255)
                          ggACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGaa                                                              >  2:200393/1‑93 (MQ=255)
                          ggACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGaa                                                              <  2:3249014/93‑1 (MQ=255)
                           gACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGaaa                                                             <  1:4043320/93‑1 (MQ=255)
                              cATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACgg                                                          >  1:1071153/1‑92 (MQ=255)
                               aTAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTg                                                         >  2:1372079/1‑93 (MQ=255)
                               aTAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTg                                                         >  1:1100580/1‑93 (MQ=255)
                                tAGACCGTCTGCTGTTTTTTAAGCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTgg                                                        >  2:2992830/1‑93 (MQ=255)
                                tAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTgg                                                        >  2:2355535/1‑93 (MQ=255)
                                 aGACCGTCTGCTGTTTTTTCACCCGACTCATGCCGGCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGACGGAAACGTGGa                                                       >  2:2035507/1‑93 (MQ=255)
                                 aGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGa                                                       >  2:4699641/1‑93 (MQ=255)
                                 aGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGa                                                       >  2:3345755/1‑93 (MQ=255)
                                  gACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGAc                                                      >  2:2691295/1‑93 (MQ=255)
                                  gACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGAc                                                      >  2:1477899/1‑93 (MQ=255)
                                     cGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTc                                                   <  2:820170/93‑1 (MQ=255)
                                     cGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTc                                                   >  1:2863754/1‑93 (MQ=255)
                                      gTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTcc                                                  >  1:3726746/1‑93 (MQ=255)
                                      gTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTcc                                                  <  2:1793271/93‑1 (MQ=255)
                                       tCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCt                                                 >  1:20020/1‑93 (MQ=255)
                                          gctgTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTgga                                              <  2:4244871/93‑1 (MQ=255)
                                           ctgTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTggag                                             <  1:2761935/93‑1 (MQ=255)
                                             gTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCCTCGCCGGAAACGTGGACTTCCTggaggg                                           >  2:1946202/1‑93 (MQ=255)
                                               tttttAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTggagggag                                         <  1:1330768/93‑1 (MQ=255)
                                                ttttAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTggagggagg                                        >  2:1673813/1‑92 (MQ=255)
                                                 tttAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGgagag                                       >  1:3765377/1‑93 (MQ=255)
                                                  ttAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGgagaga                                      <  1:1226926/93‑1 (MQ=255)
                                                   tAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGGGGCTGCCATCGCCGGCAACGTGGACTTCCTGGAGGgagagcg                                     >  1:2286708/1‑91 (MQ=255)
                                                    aaCCCCACTCATGCCGTCAAACAACTCGCCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGc                                    <  1:4381522/93‑1 (MQ=255)
                                                     aCCCGACTCATGCCGACAAACAACTCACCCACCGCAACAAACTGTAAGTGGCTGCCATAGACCGAAACGAGGACTTCCTAGAGGGAGAGAGCg                                   >  2:650737/1‑93 (MQ=255)
                                                      cccGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGt                                  <  2:4060413/93‑1 (MQ=255)
                                                       ccGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGtt                                 <  1:1088730/93‑1 (MQ=255)
                                                        cGACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGGtc                                >  2:343415/1‑93 (MQ=255)
                                                         gACTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtca                               >  2:2303290/1‑93 (MQ=255)
                                                          aCTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtcat                              >  1:4177891/1‑93 (MQ=255)
                                                           cTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGGCTTCCTGGAGGGAGAGCGCGTtcatc                             >  2:4617657/1‑93 (MQ=255)
                                                           cTCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGACGGAGCGTtgatc                             >  2:3745316/1‑93 (MQ=255)
                                                            tCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtcatca                            >  2:386465/1‑93 (MQ=255)
                                                            tCATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtcataa                            >  1:1560206/1‑91 (MQ=255)
                                                             cATGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCaa                           >  1:3370248/1‑93 (MQ=255)
                                                              aTGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAac                          <  1:1415935/93‑1 (MQ=255)
                                                               tGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAaca                         <  2:2131313/93‑1 (MQ=255)
                                                               tGCCGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAaca                         <  2:4553597/93‑1 (MQ=255)
                                                               agCCGGCAAACACCTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCACCAaca                         <  2:2524810/92‑1 (MQ=255)
                                                                  cGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGc                      <  2:723346/93‑1 (MQ=255)
                                                                  cGTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGc                      <  2:317506/93‑1 (MQ=255)
                                                                   gTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGct                     >  1:2053737/1‑93 (MQ=255)
                                                                   gTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGct                     >  1:309050/1‑93 (MQ=255)
                                                                   gTCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGct                     <  2:3275536/93‑1 (MQ=255)
                                                                    tCAAACAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAacaccccc                    >  2:3992254/1‑89 (MQ=255)
                                                                      aaaCAACTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTg                  >  1:2623150/1‑93 (MQ=255)
                                                                        acaacTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGaa                <  2:1173976/93‑1 (MQ=255)
                                                                        acaacTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGaa                <  1:2791407/93‑1 (MQ=255)
                                                                          accTCACCCACGGCAACACCCTGAAAGGGGCTGCCATCGCCGGAAACGGGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAAtt              <  2:2647298/91‑1 (MQ=255)
                                                                          aacTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGGGAGCGTTCATCAAAACGCTCTGAAtt              >  1:3505603/1‑93 (MQ=255)
                                                                          aacTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAAtt              <  1:1384013/93‑1 (MQ=255)
                                                                          aacTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAAtt              >  2:1561122/1‑93 (MQ=255)
                                                                          aacTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAAtt              <  1:1915356/93‑1 (MQ=255)
                                                                          aacTCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAAtt              <  2:2742114/93‑1 (MQ=255)
                                                                             tCACCCACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCa           >  1:794766/1‑93 (MQ=255)
                                                                               ccccACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCatt         <  1:2569170/92‑1 (MQ=255)
                                                                               ccccACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCatt         <  2:525160/92‑1 (MQ=255)
                                                                                cccACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACTTGGACTTCCTGGAGTGAGAGCGCGTTCATCGACACGCTCTGGATTTCatta        >  2:4165630/1‑93 (MQ=255)
                                                                                cccACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCatta        >  1:3116845/1‑93 (MQ=255)
                                                                                 ccACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCattat       >  1:915719/1‑93 (MQ=255)
                                                                                 ccACGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCattat       <  1:4728791/93‑1 (MQ=255)
                                                                                   aCGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTAttt     >  1:1120802/1‑93 (MQ=255)
                                                                                    cGGCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTAtttt    <  1:1179256/93‑1 (MQ=255)
                                                                                     ggCAACAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTc   <  1:3588140/93‑1 (MQ=255)
                                                                                          |                                                                                        
GGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCC  >  NC_000913/3334264‑3334442

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