Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,807,555 G→T 5.4% Q182K (CAG→AAG)  waaQ ← lipopolysaccharide core biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,807,5550GT 5.4% 410.7 / 10.0 112Q182K (CAG→AAG) waaQlipopolysaccharide core biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (46/60);  new base T (2/4);  total (48/64)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 112 total aligned reads displayed.

TTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAA  >  NC_000913/3807409‑3807703
                                                                                                                                                  |                                                                                                                                                    
ttttGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGAcg                                                                                                                                                   >  2:1352253/1‑150 (MQ=255)
       tCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTc                                                                                                                                            >  2:1544558/1‑150 (MQ=255)
          cAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTCGAAAATTTATCGTTATCCCAACACTGAAATATCTGACGCGCTGTCGGt                                                                                                                                         >  2:384722/1‑150 (MQ=255)
              cTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGa                                                                                                                                     <  1:799526/150‑1 (MQ=255)
               tCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGAt                                                                                                                                    >  2:866477/1‑150 (MQ=255)
               tCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGAt                                                                                                                                    >  2:866451/1‑150 (MQ=255)
               tCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGAt                                                                                                                                    >  2:1093508/1‑150 (MQ=255)
               tCGTGCAATCTCATCCACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGAt                                                                                                                                    <  1:1553559/150‑1 (MQ=255)
                  tGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGac                                                                                                                                 >  2:220838/1‑150 (MQ=255)
                     aaTCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGATAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGacaac                                                                                                                              <  1:546066/150‑1 (MQ=255)
                          cATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAAt                                                                                                                         <  1:1082174/150‑1 (MQ=255)
                           aTCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATg                                                                                                                        >  1:875666/1‑150 (MQ=255)
                           aTCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATg                                                                                                                        <  2:1101506/150‑1 (MQ=255)
                            tCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCTATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAAGATCTGACGAGCTGTCGGTTGGATGACACCATAATGa                                                                                                                       >  2:1306069/1‑150 (MQ=255)
                            tCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGa                                                                                                                       <  1:467898/150‑1 (MQ=255)
                               acacAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCt                                                                                                                    <  1:30352/150‑1 (MQ=255)
                                 acaAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCttt                                                                                                                  <  2:289844/150‑1 (MQ=255)
                                   aaGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTa                                                                                                                >  1:1469932/1‑150 (MQ=255)
                                     gCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACg                                                                                                              >  1:1138368/1‑150 (MQ=255)
                                      cGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGc                                                                                                             >  2:1372851/1‑150 (MQ=255)
                                      cGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGc                                                                                                             <  2:65571/150‑1 (MQ=255)
                                      cGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGc                                                                                                             >  1:1523964/1‑150 (MQ=255)
                                       gagaTCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGcc                                                                                                            >  1:1442386/1‑150 (MQ=255)
                                       gagaTCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGcc                                                                                                            >  1:1442091/1‑150 (MQ=255)
                                         gaTCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACACCATAATGATCTTTTACGCCCa                                                                                                          >  1:1079992/1‑150 (MQ=255)
                                         gaTCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCa                                                                                                          >  2:226180/1‑150 (MQ=255)
                                         gaTCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCa                                                                                                          >  2:1461189/1‑150 (MQ=255)
                                         gaTCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCa                                                                                                          >  1:1079986/1‑150 (MQ=255)
                                           tcatcTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAgg                                                                                                        <  1:1080282/150‑1 (MQ=255)
                                           tcatcTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAgg                                                                                                        <  1:291377/150‑1 (MQ=255)
                                            catcTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAggg                                                                                                       >  2:677658/1‑150 (MQ=255)
                                             atcTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTTACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGc                                                                                                      <  2:230167/150‑1 (MQ=255)
                                             atcTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTTACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGc                                                                                                      <  2:619247/150‑1 (MQ=255)
                                                tGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATc                                                                                                   >  1:533748/1‑150 (MQ=255)
                                                tGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATc                                                                                                   >  1:588029/1‑150 (MQ=255)
                                                tGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATc                                                                                                   <  2:415360/150‑1 (MQ=255)
                                                       ggCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTcc                                                                                            <  1:810136/150‑1 (MQ=255)
                                                         cccACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCg                                                                                          <  2:758880/150‑1 (MQ=255)
                                                         cccACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCg                                                                                          <  2:758866/150‑1 (MQ=255)
                                                            acaGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGgcg                                                                                       >  1:1284673/1‑150 (MQ=255)
                                                                gTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCCACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCAtc                                                                                   <  2:1103667/150‑1 (MQ=255)
                                                                      acaACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCttcttc                                                                             <  1:1543339/150‑1 (MQ=255)
                                                                      acaACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCttcttc                                                                             <  1:1543364/150‑1 (MQ=255)
                                                                      acaACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCttcttc                                                                             <  2:864938/150‑1 (MQ=255)
                                                                      acaACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCttcttc                                                                             >  2:1103640/1‑150 (MQ=255)
                                                                      acaACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCttcttc                                                                             <  2:267116/150‑1 (MQ=255)
                                                                                  ccTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCg                                                                 <  2:588029/150‑1 (MQ=255)
                                                                                  ccTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCg                                                                 <  2:533748/150‑1 (MQ=255)
                                                                                        tGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTaa                                                           <  1:836161/150‑1 (MQ=255)
                                                                                              aGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTTACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCAtt                                                     >  2:2245/1‑150 (MQ=255)
                                                                                              aGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTTACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCAtt                                                     >  2:2282/1‑150 (MQ=255)
                                                                                              aGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCAtt                                                     <  2:757927/150‑1 (MQ=255)
                                                                                              aGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCAtt                                                     <  2:1415282/150‑1 (MQ=255)
                                                                                              aGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCAtt                                                     <  2:1471471/150‑1 (MQ=255)
                                                                                               gcgcATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCAttg                                                    >  1:50967/1‑150 (MQ=255)
                                                                                               gcgcATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCAttg                                                    >  1:852228/1‑150 (MQ=255)
                                                                                                 gcATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCAttgtt                                                  <  2:1195185/150‑1 (MQ=255)
                                                                                                 gcATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCAttgtt                                                  <  2:277002/150‑1 (MQ=255)
                                                                                                       aTAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCg                                            <  2:1222190/150‑1 (MQ=255)
                                                                                                          aCCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGgtg                                         <  2:1431246/150‑1 (MQ=255)
                                                                                                             ttagccaaTTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAAGGATCTTGTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAg                                      <  2:733425/144‑1 (MQ=255)
                                                                                                             ttAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAg                                      >  2:1467211/1‑150 (MQ=255)
                                                                                                               aGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGaa                                    <  2:1650720/150‑1 (MQ=255)
                                                                                                               aGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGaa                                    <  2:1718178/150‑1 (MQ=255)
                                                                                                                gAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGaaa                                   >  2:1555320/1‑150 (MQ=255)
                                                                                                                 aaaaTTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAAt                                  <  1:1352253/150‑1 (MQ=255)
                                                                                                                 aaaaTTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAAt                                  >  1:1146704/1‑150 (MQ=255)
                                                                                                                 aaaaTTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAAt                                  >  2:1424493/1‑150 (MQ=255)
                                                                                                                 aaaaTTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAAt                                  >  1:1663429/1‑150 (MQ=255)
                                                                                                                 aaaaTTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAAt                                  >  1:1663430/1‑150 (MQ=255)
                                                                                                                 aaaaTTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAAt                                  >  2:1231857/1‑150 (MQ=255)
                                                                                                                    aTTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCgg                               <  2:1591249/150‑1 (MQ=255)
                                                                                                                    aTTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCgg                               <  2:1591308/150‑1 (MQ=255)
                                                                                                                    aTTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCgg                               <  2:267187/150‑1 (MQ=255)
                                                                                                                    aTTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCgg                               <  2:790721/150‑1 (MQ=255)
                                                                                                                     tttATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGt                              <  1:660958/150‑1 (MQ=255)
                                                                                                                     tttATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGt                              <  1:1019235/150‑1 (MQ=255)
                                                                                                                      ttATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTa                             <  2:1674228/150‑1 (MQ=255)
                                                                                                                       tATCGTTATCCCAACACTGAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTaa                            <  2:762556/150‑1 (MQ=255)
                                                                                                                           gTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGGAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATAcc                        >  1:1169197/1‑150 (MQ=255)
                                                                                                                              aTCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTaa                     >  1:1644067/1‑150 (MQ=255)
                                                                                                                              aTCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTaa                     >  1:69208/1‑150 (MQ=255)
                                                                                                                                 ccAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATgg                  <  2:940938/150‑1 (MQ=255)
                                                                                                                                 ccAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATgg                  <  2:944479/150‑1 (MQ=255)
                                                                                                                                 ccAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATgg                  <  2:1016064/150‑1 (MQ=255)
                                                                                                                                 ccAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGCGTAGAAATCGGTAATACCTAATgg                  <  2:859019/150‑1 (MQ=255)
                                                                                                                                   aaCACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCt                >  1:390563/1‑150 (MQ=255)
                                                                                                                                    acacTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTc               <  1:443631/150‑1 (MQ=255)
                                                                                                                                    acacTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTc               <  1:443632/150‑1 (MQ=255)
                                                                                                                                       cTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTCATACCGAATGGCTCAAg            >  2:662397/1‑150 (MQ=255)
                                                                                                                                       cTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTCATACCAAATGGCTCAAg            >  2:662342/1‑150 (MQ=255)
                                                                                                                                        ttAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAgg           >  2:861266/1‑150 (MQ=255)
                                                                                                                                          aaaTATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAAt                                                                                               >  2:708819/1‑64 (MQ=255)
                                                                                                                                          aaaTATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAAt                                                                                               <  1:708819/64‑1 (MQ=255)
                                                                                                                                          aaaTATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGAc         <  1:536166/150‑1 (MQ=255)
                                                                                                                                          aaaTATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGAc         <  2:923067/150‑1 (MQ=255)
                                                                                                                                           aaTATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGAcc        <  2:1585349/150‑1 (MQ=255)
                                                                                                                                              atCTTACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGAt     <  1:2282/150‑1 (MQ=255)
                                                                                                                                              atCTTACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGAt     <  1:2245/150‑1 (MQ=255)
                                                                                                                                                 tGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATaaa  <  1:797501/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAA  >  NC_000913/3807409‑3807703

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: