Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,341,542 G→A A329T (GCG→ACG)  lapB → lipopolysaccharide assembly protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,341,5420GA97.6% 130.5 / NA 41A329T (GCG→ACG) lapBlipopolysaccharide assembly protein B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/24);  minor base C (1/0);  total (17/24)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GTCAGCGAAACGCTGGAAATGTTGCAAACCTGCTATCAGCAGTTGGGTAAAACTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGC  >  NC_000913/1341299‑1341785
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
gTCAGCGAAACGCTGGAAATGTTGCAAACCTGCTATCAGCAGTTGGGTAAAACTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGga                                                                                                                                                                                                                                                 >  1:122441/1‑248 (MQ=255)
     cGAAACGCTGGAAATGTTGCAAACCTGCTATCAGCAGTTGGGTAAAACTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGc                                                                                                                                                                                                                                         >  1:338703/1‑251 (MQ=255)
          cGCTGGAAATGTTGCAAACCTGCTATCAGCAGTGGGGTAAAACGGACGAATGGGCAGAATTCCTGCAGCGGGCGGTGGAAGAGAACACCGGGGCCGATGCGGAATGGAGGCTGGCTGATATCAGCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTAGATTACGCGCCAGCGTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGct                                                                                                                                                                                                                                    <  2:186577/251‑2 (MQ=255)
                          aaCCTGCTATCAGCAGTTGGGTAAAACTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGAGGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTAAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTCCCACTTAAATGAAACGGAAGAAACGAGTGCCAAAGAGAGCCTGATgg                                                                                                                                                                                                                    >  2:455757/1‑251 (MQ=255)
                                 tatcagacgTTGGGCAAAACTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACAGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGCAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCg                                                                                                                                                                                                             >  2:268952/9‑251 (MQ=255)
                                              gTAAAACTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGCTTGGc                                                                                                                                                                                                <  2:91143/251‑1 (MQ=255)
                                                  aaCTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGATATATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCAATTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGGGCTGCGTGACATGGTTGGCg                                                                                                                                                                                               >  1:626831/1‑248 (MQ=255)
                                                          aaTGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGGACgt                                                                                                                                                                                    <  2:184113/251‑1 (MQ=255)
                                                                cAGAATTCCTGCAGCGTGCGGTGGCAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACgtagta                                                                                                                                                                                >  2:519205/1‑249 (MQ=255)
                                                                cAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACATGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAGGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAg                                                                                                                                                                              <  2:229705/251‑1 (MQ=255)
                                                                 aGAATTCCTGCAGCGTTCGGTGGAAGAGAACACCGTTGCCGATGCTGAATTGATGCGTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACATGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGGGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGGTGGCGAGAAGTTACgtagtatgc                                                                                                                                                                             <  2:386882/251‑4 (MQ=255)
                                                                              cGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCATAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTTGTAAGCCTCGTTATCGttt                                                                                                                                                                <  2:526764/251‑2 (MQ=255)
                                                                                   ggTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGa                                                                                                                                                           <  1:160589/251‑1 (MQ=255)
                                                                                       gAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGAACGTAGTAAGCATCGTTATCGTTGCCAGAAAtg                                                                                                                                                       >  1:355090/1‑251 (MQ=255)
                                                                                        aaGAGAACACCCGTGCCGATGCTGAATTGATGCTTGCTGATATCAACGAAGGACGGGACGGTAGTGACGCCGCACAGGTCAATATTACGCGCCAGCTTCAGCGACATCCGACCATGCGTGTGTTACATAAGTTAATGGATTACAACTTAAATGAACCGGAAGAAAGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGCTGGCGAGAAGGACCGTCGTAAGCCTCGTTATCGTTGCAAGCAAtgg                                                                                                                                                      >  2:29025/1‑250 (MQ=255)
                                                                                            gaACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGtt                                                                                                                                                  <  1:105812/251‑1 (MQ=255)
                                                                                                     tGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGATGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCAt                                                                                                                                         >  2:477387/1‑251 (MQ=255)
                                                                                                               gAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCctct                                                                                                                                <  1:567969/250‑1 (MQ=255)
                                                                                                               gAATGGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATGACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCTTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCctc                                                                                                                                 <  2:311087/249‑1 (MQ=255)
                                                                                                                aaTTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTAc                                                                                                                              >  2:200351/1‑251 (MQ=255)
                                                                                                                 aTTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACt                                                                                                                             <  1:77922/251‑1 (MQ=255)
                                                                                                                       gcttgctGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCAt                                                                                                                        <  1:114489/250‑1 (MQ=255)
                                                                                                                              gATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGGACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCg                                                                                                                  <  1:553734/249‑1 (MQ=255)
                                                                                                                              gATATCATCGAACCGCGCGACGGTAGTGAGGCCGCACAGGCCTATATTACGCGCCCGCTACCGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGCAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGACGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTc                                                                                                                >  2:613687/1‑251 (MQ=255)
                                                                                                                                 atcatcGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTg                                                                                                             >  1:508688/1‑251 (MQ=255)
                                                                                                                                       gAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGc                                                                                                       <  2:508677/251‑1 (MQ=255)
                                                                                                                                                         gAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAAcc                                                                                     <  1:552867/251‑1 (MQ=255)
                                                                                                                                                              cGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTc                                                                                >  1:504682/1‑251 (MQ=255)
                                                                                                                                                                  caGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACTATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACAGGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCGCTACTGGCATTGTCCGTCTTGTCGTGCCGGGTCAACCATTAAACCGATTCGCgg                                                                            <  1:221198/251‑1 (MQ=255)
                                                                                                                                                                        tataTTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTgt                                                                      >  2:166814/1‑250 (MQ=255)
                                                                                                                                                                         atatTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCtt                                                                        <  1:575519/248‑1 (MQ=255)
                                                                                                                                                                                 ggcCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGGCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTGGATGGCCTGTa                                                             <  2:31446/250‑1 (MQ=255)
                                                                                                                                                                                                          tGCGTGTGTTCCATAAGTTCATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACttt                                      >  1:204271/1‑249 (MQ=255)
                                                                                                                                                                                                                     cATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTTAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAAGTTTTAAAAAAAATCCGACTTTAGTTACAACATAc                         <  2:122441/251‑1 (MQ=255)
                                                                                                                                                                                                                        aaGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTGTTAAAAAAAATCCGACGTTAGTTACAACATACTaa                      <  2:232737/251‑1 (MQ=255)
                                                                                                                                                                                                                         aGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGTCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGTTGTTACCTCATACACCCTCTACTGGCATTGTCCGGCTTGTCTGGCGTGGTCAACCATTAAACCGATTCGCGGTCTTGATTGCCTGTAATTGTTAAAAAAGATCCGAGTTGAGTTTCAACTTTCTAat                     <  2:326114/251‑1 (MQ=255)
                                                                                                                                                                                                                           ttAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAattt                   >  1:180270/1‑250 (MQ=255)
                                                                                                                                                                                                                           ttAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTCGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAatt                    >  1:183822/1‑250 (MQ=255)
                                                                                                                                                                                                                                      ccACTTAAATGAAACGGAAGAAGGGCTTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGTTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCa        <  1:463003/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ccACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCa        <  1:268952/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aaaTGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGc  <  1:477387/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
GTCAGCGAAACGCTGGAAATGTTGCAAACCTGCTATCAGCAGTTGGGTAAAACTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGC  >  NC_000913/1341299‑1341785

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