Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,796,013 G→A 5.9% G12D (GGC→GAC)  waaC → ADP‑heptose:LPS heptosyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,796,0130GA 5.9% 758.7 / 21.8 236G12D (GGC→GAC) waaCADP‑heptose:LPS heptosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (104/118);  new base A (5/9);  total (109/127)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 240 total aligned reads displayed.

CAAAGTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCT  >  NC_000913/3795864‑3796162
                                                                                                                                                     |                                                                                                                                                     
caAAGTGCGTAAAGGTGACGCTGCGGAGGGTTATCATCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATggg                                                                                                                                                       <  2:1173546/150‑1 (MQ=255)
caAAGTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATggg                                                                                                                                                       >  2:1783115/1‑150 (MQ=255)
 aaaGTTCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCATGACGGATGCGTGTTTTGATCGTTAAAACATCGTCGATGGGc                                                                                                                                                      <  2:812142/150‑1 (MQ=255)
    gTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGTAATAACTCAACGCGCTATTGTTAAAAGAGGAAGCCTGACGGATGCTTGTTTTGATGGTTAAAATATCGTCCATCGACGAt                                                                                                                                                   <  2:883628/150‑1 (MQ=255)
       cGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGtt                                                                                                                                                >  2:1567644/1‑150 (MQ=255)
       cGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGtt                                                                                                                                                >  1:1492287/1‑150 (MQ=255)
         tAAAGGTGACGCTGCGGATGGTTATCACCAGAGCTTAATCGACATTACTGCCCAGCGCTTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGCCCATCGACGATGtttt                                                                                                                                              <  2:1043656/150‑3 (MQ=255)
         tAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGAt                                                                                                                                                   >  1:1461764/1‑145 (MQ=255)
         tAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTtct                                                                                                                                              >  1:301162/1‑150 (MQ=255)
         tAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTtct                                                                                                                                              >  1:338826/1‑150 (MQ=255)
          aaaGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTtctc                                                                                                                                             <  1:1821949/150‑1 (MQ=255)
           aaGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGTGTTTTTATCGTTAAAACTTCTTCGATGGGCGATGTTCTcc                                                                                                                                            <  2:1514450/150‑1 (MQ=255)
           aaGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTcc                                                                                                                                            <  2:1258722/150‑1 (MQ=255)
           aaGGGGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTcc                                                                                                                                            <  2:1677888/150‑1 (MQ=255)
            aGGTGACGCTGTGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGCCCATCGACGATGttttaac                                                                                                                                           <  1:1099090/150‑6 (MQ=255)
            aGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGCCCATCGACGATGttttaac                                                                                                                                           <  1:650334/150‑6 (MQ=255)
             ggTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCAt                                                                                                                                          <  1:205168/150‑1 (MQ=255)
             ggTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGCCCATCGACGATGttttaacg                                                                                                                                          <  2:1353173/150‑7 (MQ=255)
              gTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATa                                                                                                                                         <  2:1709231/150‑1 (MQ=255)
                 aCGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGt                                                                                                                                      <  1:152020/150‑1 (MQ=255)
                  cGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGtt                                                                                                                                     >  1:1657697/1‑150 (MQ=255)
                   gCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTg                                                                                                                                    <  1:855760/150‑1 (MQ=255)
                    cTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGc                                                                                                                                   >  1:1416013/1‑150 (MQ=255)
                     tGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGcc                                                                                                                                  >  2:114739/1‑150 (MQ=255)
                     tGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGcc                                                                                                                                  <  1:1200512/150‑1 (MQ=255)
                      gCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGccc                                                                                                                                 <  1:1303266/150‑1 (MQ=255)
                            ggTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTGTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGcactc                                                                                                                           <  1:1476439/150‑1 (MQ=255)
                              ttATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGcactcac                                                                                                                         <  2:1080666/150‑1 (MQ=255)
                                   aCCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATg                                                                                                                    >  2:49285/1‑150 (MQ=255)
                                   aCCAGAGCTGAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGACGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATg                                                                                                                    >  2:932265/1‑150 (MQ=255)
                                    ccAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGc                                                                                                                   <  1:544455/150‑1 (MQ=255)
                                           ttAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCcagcag                                                                                                            <  1:1533752/150‑1 (MQ=255)
                                             aaTCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGc                                                                                                          <  2:1416013/150‑1 (MQ=255)
                                             aaTCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGc                                                                                                          <  1:329843/150‑1 (MQ=255)
                                                cGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCg                                                                                                                                                     <  1:1747347/104‑1 (MQ=255)
                                                cGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCg                                                                                                                                                     >  2:1747347/1‑104 (MQ=255)
                                                      tACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAgg                                                                                                 <  1:44168/150‑1 (MQ=255)
                                                        cTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGa                                                                                               >  2:181881/1‑150 (MQ=255)
                                                          ccccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGAtt                                                                                             >  2:1996066/1‑150 (MQ=255)
                                                          ccccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGAtt                                                                                             >  1:1773026/1‑150 (MQ=255)
                                                                        ggAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggt                                                                               >  2:1447900/1‑150 (MQ=255)
                                                                        ggAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggt                                                                               >  2:1903084/1‑150 (MQ=255)
                                                                          aagaaCTCAACGCGCTATTGTTACACGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTTTCTATAGGTTTCCGGCACTCACTGATTCTAATCAGGCATTCCCTGTGATTAATTTTGACTGtgtgg                                                                             <  2:1775365/150‑1 (MQ=255)
                                                                          aagaaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggtgg                                                                             <  1:114739/150‑1 (MQ=255)
                                                                              aCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTGACTGATGcc                                                                                                                  >  1:1894872/1‑109 (MQ=255)
                                                                               cTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaa                                                                        <  2:1724815/150‑1 (MQ=255)
                                                                               cTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTCGTGgaa                                                                        >  1:1009925/1‑150 (MQ=255)
                                                                                 cAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGACGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaaga                                                                      >  2:1627831/1‑150 (MQ=255)
                                                                                    cgcgCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAgg                                                                   >  1:1296886/1‑150 (MQ=255)
                                                                                      cgcTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGCGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCAAGGGATAAAGTTTGACTGGGTGGTGGAAGAAGGGt                                                                 >  2:243283/1‑150 (MQ=255)
                                                                                        cTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAgg                                                                   <  2:216213/146‑1 (MQ=255)
                                                                                        cTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAgg                                                                   >  1:216213/1‑146 (MQ=255)
                                                                                         tATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATCAAGTTTGACTGGGTGGTGGAAGAAGGGTTCg                                                              <  1:572195/150‑1 (MQ=255)
                                                                                          aTTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGc                                                             <  1:403842/150‑1 (MQ=255)
                                                                                          aTTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAAACCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGc                                                             >  2:598143/1‑150 (MQ=255)
                                                                                           ttGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGca                                                            <  2:1053156/150‑1 (MQ=255)
                                                                                            tGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaa                                                                        <  2:424498/137‑1 (MQ=255)
                                                                                             gTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGcaca                                                          <  2:840302/150‑1 (MQ=255)
                                                                                               tACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGa                                                        >  1:1585665/1‑150 (MQ=255)
                                                                                                aCAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGc                                                                                                          <  1:1719399/99‑1 (MQ=255)
                                                                                                aCAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATccc                                                                                                    <  1:832686/105‑1 (MQ=255)
                                                                                                aCAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAt                                                       <  2:787221/150‑1 (MQ=255)
                                                                                                 cAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAtt                                                      <  1:1947945/150‑1 (MQ=255)
                                                                                                   agagGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttcc                                                    >  2:6186/1‑150 (MQ=255)
                                                                                                       gAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttccttcc                                                <  1:879055/150‑1 (MQ=255)
                                                                                                        aaGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATtccttcct                                               <  1:105971/150‑1 (MQ=255)
                                                                                                            cTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCa                                           >  1:1969042/1‑150 (MQ=255)
                                                                                                             tGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGcaca                                                          >  2:1887468/1‑134 (MQ=255)
                                                                                                              gACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACg                                         <  2:338826/150‑1 (MQ=255)
                                                                                                                   aTGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGcc                                    <  1:607836/150‑1 (MQ=255)
                                                                                                                     gCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGt                                  >  1:162604/1‑150 (MQ=255)
                                                                                                                      cGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGtt                                 <  2:309525/150‑1 (MQ=255)
                                                                                                                       gggTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGAAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGTCGTTg                                >  2:1324753/1‑150 (MQ=255)
                                                                                                                        ggTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGa                               <  1:618640/150‑1 (MQ=255)
                                                                                                                         gTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCTGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCAGGCTGCCGTTGAg                              <  1:1994227/150‑1 (MQ=255)
                                                                                                                         gTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGgtggtg                                                                           >  2:265204/1‑105 (MQ=255)
                                                                                                                         gTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGgtggtg                                                                           <  1:265204/105‑1 (MQ=255)
                                                                                                                         gTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAg                              >  1:761448/1‑150 (MQ=255)
                                                                                                                             tGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTg                                              <  2:203308/130‑1 (MQ=255)
                                                                                                                             tGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAg                          >  1:1045140/1‑150 (MQ=255)
                                                                                                                              gATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGt                         >  1:695741/1‑150 (MQ=255)
                                                                                                                              gATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGt                         >  1:1848487/1‑150 (MQ=255)
                                                                                                                              gATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGt                         >  1:1253025/1‑150 (MQ=255)
                                                                                                                               aTCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGtt                        <  1:49285/150‑1 (MQ=255)
                                                                                                                               aTCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGtt                        <  2:1561760/150‑1 (MQ=255)
                                                                                                                                 cGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTAt                      >  2:62446/1‑150 (MQ=255)
                                                                                                                                  gTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTAtt                     >  1:1365482/1‑150 (MQ=255)
                                                                                                                                   ggagaaCATCGCCCATCGACGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCATGCTGCCGTTGAGCGAGTTATTc                    >  2:738470/5‑150 (MQ=255)
                                                                                                                                   ggagaaCATCGCCCATCGACGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTc                    >  2:1278889/5‑150 (MQ=255)
                                                                                                                                    gagaaCATCGCCCATCGACGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTcc                   >  2:1135360/4‑150 (MQ=255)
                                                                                                                                      gaaCATCGCCCATCGACGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCtg                 >  2:1388424/2‑150 (MQ=255)
                                                                                                                                       aaCATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGGGGAAGAAGGGGTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCtgt                <  2:1939658/150‑1 (MQ=255)
                                                                                                                                        aCATCGCCCATCGACGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCtgtg               <  1:1278889/150‑1 (MQ=255)
                                                                                                                                          aTCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGc             <  2:549837/150‑1 (MQ=255)
                                                                                                                                          aTCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGc             >  1:1085090/1‑150 (MQ=255)
                                                                                                                                           tcgtcgACGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCa            <  1:932265/150‑1 (MQ=255)
                                                                                                                                             gtcgATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAAt          >  2:1605206/1‑150 (MQ=255)
                                                                                                                                               cgATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATAc        >  2:758687/1‑150 (MQ=255)
                                                                                                                                                 aTGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGt      <  1:1697792/150‑1 (MQ=255)
                                                                                                                                                    ggCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGc   <  1:1996066/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAAGTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCT  >  NC_000913/3795864‑3796162

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