Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,803,395 A→G 5.7% L248S (TTG→TCG)  waaB ← lipopolysaccharide 1,6‑galactosyltransferase; UDP‑D‑galactose:(glucosyl)lipopolysaccharide‑1, 6‑D‑galactosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,803,3950AG 5.7% 233.6 / 4.1 72L248S (TTG→TCG) waaBlipopolysaccharide 1,6‑galactosyltransferase; UDP‑D‑galactose:(glucosyl)lipopolysaccharide‑1, 6‑D‑galactosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base G (2/2);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TTGAATAATGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCA  >  NC_000913/3803250‑3803541
                                                                                                                                                 |                                                                                                                                                  
ttGAATAATGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTc                                                                                                                                                <  1:1021558/150‑1 (MQ=255)
ttGAATAATGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTc                                                                                                                                                >  2:1080149/1‑150 (MQ=255)
  gAATAATGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTccc                                                                                                                                              >  2:1433271/1‑150 (MQ=255)
   aataatGTCAGCAGGCCCGGAGACACAGTCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGCCCAATTCCCa                                                                                                                                             <  1:538522/150‑1 (MQ=255)
   aataatGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCa                                                                                                                                             <  1:704613/150‑1 (MQ=255)
    ataatGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAg                                                                                                                                            >  2:21778/1‑150 (MQ=255)
    ataatGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAg                                                                                                                                            <  2:546966/150‑1 (MQ=255)
    ataatGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAg                                                                                                                                            >  1:884581/1‑150 (MQ=255)
    ataatGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAg                                                                                                                                            >  2:1224950/1‑150 (MQ=255)
    ataatGTCAGCAGGCCCGGAGACACAAACGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTACCAg                                                                                                                                            >  2:295383/1‑150 (MQ=255)
      aatGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAggg                                                                                                                                          >  1:302713/1‑150 (MQ=255)
      aatGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAggg                                                                                                                                          >  1:1239575/1‑150 (MQ=255)
        tGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGat                                                                                                                                        >  1:439405/1‑150 (MQ=255)
        tGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGat                                                                                                                                        >  1:832719/1‑150 (MQ=255)
              cagGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATAttgtt                                                                                                                                  <  2:1209645/150‑1 (MQ=255)
                ggCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGc                                                                                                                                <  2:572188/150‑1 (MQ=255)
                    cGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAAc                                                                                                                            <  2:1509879/150‑1 (MQ=255)
                    cGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAAc                                                                                                                            <  1:605601/150‑1 (MQ=255)
                         acacAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATAc                                                                                                                       >  1:523471/1‑150 (MQ=255)
                          cacaATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATAcc                                                                                                                      >  1:1161654/1‑150 (MQ=255)
                          cacaATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATAcc                                                                                                                      >  1:1161628/1‑150 (MQ=255)
                            caATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAg                                                                                                                    >  2:1247829/1‑150 (MQ=255)
                             aaTCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGa                                                                                                                   >  1:738279/1‑150 (MQ=255)
                                   cACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCGATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTc                                                                                                             >  2:521274/1‑150 (MQ=255)
                                   cACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCGATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTc                                                                                                             >  2:1382584/1‑150 (MQ=255)
                                   cACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTc                                                                                                             <  1:1422859/150‑1 (MQ=255)
                                    aCTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCg                                                                                                            >  1:620570/1‑150 (MQ=255)
                                         tACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGatcat                                                                                                       <  1:1066676/150‑1 (MQ=255)
                                               gAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAatat                                                                                                 <  2:528276/150‑1 (MQ=255)
                                                 aTCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATAttt                                                                                               <  1:1158936/150‑1 (MQ=255)
                                                     cccACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAAtt                                                                                           >  2:1318219/1‑150 (MQ=255)
                                                       cACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATtct                                                                                         <  2:83373/150‑1 (MQ=255)
                                                       cACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATtct                                                                                         <  1:1074490/150‑1 (MQ=255)
                                                       cACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATtct                                                                                         <  2:1415507/150‑1 (MQ=255)
                                                       cACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATtct                                                                                         <  2:1256501/150‑1 (MQ=255)
                                                        aCGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATtctc                                                                                        <  1:669116/150‑1 (MQ=255)
                                                           aCAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTAc                                                                                     <  1:1349663/150‑1 (MQ=255)
                                                             aaGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAANAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCa                                                                                   >  2:37127/1‑150 (MQ=255)
                                                                 cTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAg                                                                               >  1:1228602/1‑150 (MQ=255)
                                                                 cTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAg                                                                               <  2:1045381/150‑1 (MQ=255)
                                                                 cTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTCAATTCTCTACCATAAg                                                                               >  1:1228594/1‑150 (MQ=255)
                                                                    cAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCt                                                                            <  1:462909/150‑1 (MQ=255)
                                                                    cAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCt                                                                            <  2:1451782/150‑1 (MQ=255)
                                                                      aGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTgg                                                                          <  1:407119/150‑1 (MQ=255)
                                                                      aGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTgg                                                                          <  1:1080149/150‑1 (MQ=255)
                                                                       gTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGc                                                                         <  2:144346/150‑1 (MQ=255)
                                                                           ggTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCAttt                                                                     >  1:1455465/1‑150 (MQ=255)
                                                                             tCATTGGGAAACCTTCAAAAGATGACGTGAGCGATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCAttttt                                                                   <  1:1129544/150‑1 (MQ=255)
                                                                               aTTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCg                                                                 >  2:269435/1‑150 (MQ=255)
                                                                                  gggAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGaaa                                                              >  2:38435/1‑150 (MQ=255)
                                                                                     aaaCCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCg                                                           >  2:543422/1‑150 (MQ=255)
                                                                                       aCCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGa                                                         <  2:832719/150‑1 (MQ=255)
                                                                                       aCCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGa                                                         >  1:1660150/1‑150 (MQ=255)
                                                                                          ttCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGcc                                                      >  1:1115042/1‑150 (MQ=255)
                                                                                           tCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCa                                                     >  1:26927/1‑150 (MQ=255)
                                                                                                  aTGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCGATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCa                                              <  1:521274/150‑1 (MQ=255)
                                                                                                  aTGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCGATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCa                                              <  1:1382584/150‑1 (MQ=255)
                                                                                                     aCGTGAGCAATAACGCAGTCAGTTTTTCAATATCTTGCTGGACCTATTCGCAGGGATTTTGTTGCCAAGCATACGAGACATTTCGATCATCAATATTTATTTCTCTACGATAAGCCTGGCATTTTTCGAAATCGGAGCCATGACCCAAtc                                           <  2:1015310/150‑2 (MQ=255)
                                                                                                     aCGTGAGCAATAACGCACTCACTTTTTCAATCTCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATa                                           <  2:1015368/150‑1 (MQ=255)
                                                                                                            cAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAg                                    >  2:1226880/1‑150 (MQ=255)
                                                                                                                  cGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCa                              <  2:866458/150‑1 (MQ=255)
                                                                                                                   gCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAg                             >  1:659299/1‑150 (MQ=255)
                                                                                                                     actcactTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGtt                           <  2:1117985/150‑1 (MQ=255)
                                                                                                                         actTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCt                       >  1:434961/1‑150 (MQ=255)
                                                                                                                             tttCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCg                   >  2:1653611/1‑150 (MQ=255)
                                                                                                                             tttCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCg                   <  2:594793/150‑1 (MQ=255)
                                                                                                                             tttCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCg                   >  1:1262637/1‑150 (MQ=255)
                                                                                                                                cAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCtt                <  1:780105/150‑1 (MQ=255)
                                                                                                                                cAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATAACAGACAATTCGATCATCAATATTTAAGTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCtt                <  1:780103/150‑1 (MQ=255)
                                                                                                                                  atatCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTtgt              <  1:1171863/150‑1 (MQ=255)
                                                                                                                                           tGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCa     >  1:98266/1‑150 (MQ=255)
                                                                                                                                              aCCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCa  >  2:735749/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TTGAATAATGTCAGCAGGCCCGGAGACACAATCGGCACTGATACAAGGAATCCCCCACGACAAGGCTTCAAGTAGGGTCATTGGGAAACCTTCAAAAGATGACGTGAGCAATAACGCACTCACTTTTTCAATATCTTGCTGGACCAATTCCCAGGGATATTGTTGCCAACCATACCAGACAATTCGATCATCAATATTTAATTCTCTACCATAAGCCTGGCATTTTTCGAAATCGGAGCCATCACCCAATACATGGAGTTTCCAGTTACCTTTCGCTTGTGATAAGCCATCA  >  NC_000913/3803250‑3803541

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