Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,798,234 A→T 12.5% intergenic (+27/+5) waaL → / ← waaU O‑antigen ligase/lipopolysaccharide core biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,798,2340AT12.5% 577.8 / 63.8 184intergenic (+27/+5)waaL/waaUO‑antigen ligase/lipopolysaccharide core biosynthesis
Reads supporting (aligned to +/- strand):  ref base A (80/81);  new base T (12/11);  total (92/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

Only 100 of 199 total aligned reads displayed.

TGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCG  >  NC_000913/3798087‑3798381
                                                                                                                                                   |                                                                                                                                                   
tGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaa                                                                                                                                                    >  1:1511610/1‑149 (MQ=255)
  gCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaaa                                                                                                                                                  <  1:1179370/149‑1 (MQ=255)
  gCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaaa                                                                                                                                                  <  2:626440/149‑1 (MQ=255)
   cATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAAc                                                                                                                                                 <  1:1772902/149‑1 (MQ=255)
            gATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtctt                                                                                                                                        <  2:1226529/149‑1 (MQ=255)
                    tCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaaa                                                                                                                                                  <  4:888855/131‑1 (MQ=255)
                    tCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCActtc                                                                                                                                          >  4:1430627/1‑139 (MQ=255)
                    tCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCActtc                                                                                                                                          <  3:1430627/139‑1 (MQ=255)
                      aGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa                                                                                                                                              <  2:1685626/133‑1 (MQ=255)
                      aGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa                                                                                                                                              >  1:1685626/1‑133 (MQ=255)
                       gtgatgtgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAAc                                                                                                                             >  4:1071093/1‑149 (MQ=255)
                          atgtgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTa                                                                                                                          <  3:1071093/149‑1 (MQ=255)
                             tgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTaa                                                                                                                         >  3:89489/1‑147 (MQ=255)
                              gatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAtct                                                                                                                      >  4:917653/1‑149 (MQ=255)
                               atCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTa                                                                                                                                       <  3:565468/131‑1 (MQ=255)
                               atCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTa                                                                                                                                       >  4:565468/1‑131 (MQ=255)
                                  atatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                                                                  >  2:1034471/1‑149 (MQ=255)
                                        gCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCAc                                                                                                                                             <  4:689208/116‑1 (MQ=255)
                                          aCGCAGCATTCCAATTATCATTATATCCGCTATAGTTCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAct                                                                                                                                            <  4:1363277/115‑1 (MQ=255)
                                          aCGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAct                                                                                                                                            >  3:1363277/1‑115 (MQ=255)
                                            gcagcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTCACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa                                                                                                                                              >  3:1838371/1‑111 (MQ=255)
                                            gcagcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTACCAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCAc                                                                                                                                             <  1:396487/112‑1 (MQ=255)
                                             cagcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGt                                                                                                       <  3:692385/149‑1 (MQ=255)
                                                         tATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACt                                                                                                                            >  2:948220/1‑116 (MQ=255)
                                                            cATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAc                                                                                                                                             >  4:407241/1‑96 (MQ=255)
                                                             aTTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtc                                                                                                                                          <  3:387454/98‑1 (MQ=255)
                                                              ttATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa                                                                                                                                              <  1:1473203/93‑1 (MQ=255)
                                                              ttATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa                                                                                                                                              >  3:989359/1‑93 (MQ=255)
                                                              ttATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGa                                                                                                             <  1:815073/126‑1 (MQ=255)
                                                                 tatCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtc                                                                                                                                          <  2:1526176/94‑1 (MQ=255)
                                                                   tCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCtt                                                                                                                   <  1:10025/115‑1 (MQ=255)
                                                                   tCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCtt                                                                                                                   >  2:10025/1‑115 (MQ=255)
                                                                   tCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGat                                                                                     <  2:1072958/145‑1 (MQ=255)
                                                                     cGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAtct                                                                                                                      >  4:31948/1‑110 (MQ=255)
                                                                      gCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAt                                                                                                                        >  2:973188/1‑107 (MQ=255)
                                                                           aGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTAc                                                                                                                <  3:367122/110‑1 (MQ=255)
                                                                            gTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaaa                                                                                                                                                  <  1:1007351/75‑1 (MQ=255)
                                                                            gTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaaa                                                                                                   <  2:396121/122‑1 (MQ=255)
                                                                               ctctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACt                                                                                                               >  1:1297638/1‑107 (MQ=255)
                                                                               ctctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGat                                                                                     >  2:726097/1‑133 (MQ=255)
                                                                                tctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAc                                                                                                                                             <  3:190317/76‑1 (MQ=255)
                                                                                tctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                                                                              <  1:1827861/123‑1 (MQ=255)
                                                                                tctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                                                                              >  2:1827861/1‑123 (MQ=255)
                                                                                  ttACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTAt                                                                                                                                      <  2:494028/81‑1 (MQ=255)
                                                                                  ttACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                                                                              >  4:783547/1‑121 (MQ=255)
                                                                                  ttACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                                                                              <  3:783547/121‑1 (MQ=255)
                                                                                      tCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCtata                                                                                                                                 <  2:1480571/82‑1 (MQ=255)
                                                                                      tCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATctct                                                                                                                    <  1:956050/95‑1 (MQ=255)
                                                                                      tCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGt                                                                                                       >  1:1085268/1‑108 (MQ=255)
                                                                                         tCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                                                                 <  4:1613886/143‑1 (MQ=255)
                                                                                         tCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATATTg                                                                 >  3:1613886/1‑143 (MQ=255)
                                                                                           aTTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTa                                                                                                                                       <  3:905289/71‑1 (MQ=255)
                                                                                           aTTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTatat                                                                     >  3:847847/1‑137 (MQ=255)
                                                                                              aataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGat                                                                                                            >  3:1076289/1‑94 (MQ=255)
                                                                                              aataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGaa                                                                                                            <  4:1076289/95‑1 (MQ=255)
                                                                                              aataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                                                                              <  1:135285/109‑1 (MQ=255)
                                                                                              aataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                                                                              >  2:135285/1‑109 (MQ=255)
                                                                                               ataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCtt                                                                                                                   >  4:666588/1‑87 (MQ=255)
                                                                                                   tCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCa                                                                                             <  2:1168127/105‑1 (MQ=255)
                                                                                                   tCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTa                                                                             >  3:455023/1‑121 (MQ=255)
                                                                                                   tCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg                                                               >  3:1499991/1‑135 (MQ=255)
                                                                                                   tCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg                                                               <  4:1499991/135‑1 (MQ=255)
                                                                                                     gTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaaa                                                                                                   <  1:1232228/97‑1 (MQ=255)
                                                                                                       aaCAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATAc                                                          >  3:1568142/1‑136 (MQ=255)
                                                                                                       aaCAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATAc                                                          <  4:1568142/136‑1 (MQ=255)
                                                                                                                attaattaAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATa                                                           >  3:1167518/1‑126 (MQ=255)
                                                                                                                attaattaAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATa                                                           <  4:1167518/126‑1 (MQ=255)
                                                                                                                attaattaAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt                                                     >  3:1703819/1‑132 (MQ=255)
                                                                                                                attaattaAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt                                                     <  4:1703819/132‑1 (MQ=255)
                                                                                                                     ttaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGtt                                                                                                >  3:1778692/1‑84 (MQ=255)
                                                                                                                     ttaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGtt                                                                                                <  4:1778692/84‑1 (MQ=255)
                                                                                                                     ttaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCa                                                                                             <  2:1299117/87‑1 (MQ=255)
                                                                                                                      taaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCggg                              <  4:896010/149‑1 (MQ=255)
                                                                                                                         gAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGtt                                                                  >  1:692217/1‑110 (MQ=255)
                                                                                                                          aaTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCt                                                                                 <  3:631742/94‑1 (MQ=255)
                                                                                                                          aaTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCt                                                                                 >  4:631742/1‑94 (MQ=255)
                                                                                                                          aaTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg                                                               >  1:1142512/1‑112 (MQ=255)
                                                                                                                          aaTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                                                                 >  3:681980/1‑110 (MQ=255)
                                                                                                                             aaaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                                                                 <  1:1438643/107‑1 (MQ=255)
                                                                                                                             aaaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                                                                 >  2:1438643/1‑107 (MQ=255)
                                                                                                                             aaaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                                                                 <  3:756816/107‑1 (MQ=255)
                                                                                                                             aaaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                                                                 >  4:756816/1‑107 (MQ=255)
                                                                                                                              aaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGaa                                                                                                            >  3:302741/1‑63 (MQ=255)
                                                                                                                               aCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgg                                                                <  4:1863629/106‑1 (MQ=255)
                                                                                                                                    tttAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAg                                                                                                  >  1:1617736/1‑67 (MQ=255)
                                                                                                                                    tttAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCgg                               <  4:803023/134‑1 (MQ=255)
                                                                                                                                      tAAGAAGTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTg                                                    >  1:1063339/1‑111 (MQ=255)
                                                                                                                                         gaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg                                                               <  4:1107188/97‑1 (MQ=255)
                                                                                                                                         gaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg                                                               >  3:1107188/1‑97 (MQ=255)
                                                                                                                                          aagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaaa                                                                                                   >  4:376172/1‑60 (MQ=37)
                                                                                                                                          aagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAAt                 >  3:1321624/1‑142 (MQ=255)
                                                                                                                                          aagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAAt                 <  4:1321624/142‑1 (MQ=255)
                                                                                                                                           agTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGa         >  2:1726129/1‑149 (MQ=255)
                                                                                                                                            gTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGAttt                                           >  2:1265735/1‑114 (MQ=255)
                                                                                                                                             tGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGata       >  1:300371/1‑149 (MQ=255)
                                                                                                                                                 ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAAt                 <  1:87764/135‑1 (MQ=255)
                                                                                                                                                 ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAg            <  2:1088501/140‑1 (MQ=255)
                                                                                                                                                 ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTcc   >  4:1813362/1‑149 (MQ=255)
                                                                                                                                                  tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCg  <  4:246561/149‑1 (MQ=255)
                                                                                                                                                   aaaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAgg           >  2:865376/1‑139 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCG  >  NC_000913/3798087‑3798381

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