Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,288,033 T→G 6.4% intergenic (+23/‑57) yraH → / → yraI putative fimbrial‑like adhesin protein/putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,288,0330TG 6.4% 751.7 / 33.2 229intergenic (+23/‑57)yraH/yraIputative fimbrial‑like adhesin protein/putative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (105/108);  new base G (7/8);  total (113/116)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 231 total aligned reads displayed.

CCCGGCGGAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACG  >  NC_000913/3287884‑3288182
                                                                                                                                                     |                                                                                                                                                      
cccGGCGGAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCt                                                                                                                                                        <  1:3601649/150‑1 (MQ=255)
  cggcggAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGt                                                                                                                                                      >  1:2401480/1‑150 (MQ=255)
     cggAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAg                                                                                                                                                   <  2:1386942/150‑1 (MQ=255)
       gAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGAt                                                                                                                                                 <  1:1000914/150‑1 (MQ=255)
         aTCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATaa                                                                                                                                               <  2:1790109/150‑1 (MQ=255)
         aTCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATaa                                                                                                                                               <  1:2808822/150‑1 (MQ=255)
         aTCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATaa                                                                                                                                               <  2:1326081/150‑1 (MQ=255)
         aTCCAATGGTATTGTCAATGCTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCCTCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATaa                                                                                                                                               <  2:3124626/150‑1 (MQ=255)
          tCCACTGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAAc                                                                                                                                              <  1:3109844/150‑1 (MQ=255)
          tCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAAc                                                                                                                                              <  1:475733/150‑1 (MQ=255)
          tCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAAc                                                                                                                                              <  2:1627874/150‑1 (MQ=255)
           ccAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACg                                                                                                                                             <  1:3320823/150‑1 (MQ=255)
            cAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACGc                                                                                                                                            <  2:2365951/150‑1 (MQ=255)
             aaTGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGcc                                                                                                                                           >  2:2557299/1‑150 (MQ=255)
             aaTGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGGTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGcc                                                                                                                                           <  2:3950152/150‑1 (MQ=255)
                 gTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCttt                                                                                                                                       >  1:860129/1‑150 (MQ=255)
                  tATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCtttt                                                                                                                                      >  1:2542523/1‑150 (MQ=255)
                  tATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCtttt                                                                                                                                      >  2:3648208/1‑150 (MQ=255)
                   aTTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCACTTttc                                                                                                                                     >  2:2109364/1‑150 (MQ=255)
                      gTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttcctt                                                                                                                                  >  1:1033154/1‑150 (MQ=255)
                      gTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttcctt                                                                                                                                  <  2:943628/150‑1 (MQ=255)
                      gTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttcctt                                                                                                                                  >  2:3485318/1‑150 (MQ=255)
                       tCAATGTTGACCTCATCGCTTACTCCGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttccttc                                                                                                                                 <  2:3144143/150‑1 (MQ=255)
                         aaTGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGGTAGGGCGTTATCTGCCCGGTCAGATACCGCCCTGTtccgtcct                                                                                                                               >  2:1072664/1‑150 (MQ=255)
                          aTGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCctc                                                                                                                              <  2:3482689/150‑1 (MQ=255)
                             ttGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCttt                                                                                                                           >  2:1864001/1‑150 (MQ=255)
                               gACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTtct                                                                                                                         >  1:2307140/1‑150 (MQ=255)
                                aCCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTtctc                                                                                                                        >  1:2214362/1‑150 (MQ=255)
                                   tCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGtt                                                                                                                     >  2:1790813/1‑150 (MQ=255)
                                     aTCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGt                                                                                                                   >  1:2900436/1‑150 (MQ=255)
                                       cGCTTACAACGAAACAACCGTAGTTCCAACAGAAAACACCCCACGGGAAGCACACGCAAACGTAAACGTAACGCTGCATTACCGTTAATACGTTAGGGCGTTAACTGACCTGTCAGATAACGCCCATTTCCTTCCTCTTTCACGTTGTAt                                                                                                                 <  2:3998025/150‑1 (MQ=255)
                                        gCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATc                                                                                                                >  2:2304745/1‑150 (MQ=255)
                                          ttACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAg                                                                                                              >  1:2754777/1‑150 (MQ=255)
                                           tactacGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAgg                                                                                                             >  2:976314/1‑150 (MQ=255)
                                            actacGTATCAACCGTAGTTCCATCAGAAATCACCCCCGGGGAAGCAGACGCTATCGTTAACGTAACGCTGTATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGt                                                                                                            <  1:681670/150‑1 (MQ=255)
                                               ccGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑a                                                                                                        <  1:3632226/149‑1 (MQ=255)
                                               acGTATCAACCGTAGTTCCATCAGAAATCCCCCCAGGGGACGCAGACGATATCGTTAACGTAACGCGGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑a                                                                                                        <  1:439043/150‑1 (MQ=255)
                                               acGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑a                                                                                                        >  1:2873325/1‑150 (MQ=255)
                                                    tCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCATATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAAt                                                                                                   <  1:2813438/150‑1 (MQ=255)
                                                    tCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAAt                                                                                                   >  1:2061775/1‑150 (MQ=255)
                                                      aaCCGTAGTTCCATCAGAAATCACCCCCGGGGAAGCAGACGCCATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATaa                                                                                                 <  1:483644/150‑1 (MQ=255)
                                                       aCCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAAc                                                                                                >  1:3979884/1‑150 (MQ=255)
                                                           tAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATg                                                                                            >  1:817513/1‑150 (MQ=255)
                                                              ttCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCa                                                                                         >  1:1777476/1‑150 (MQ=255)
                                                               tCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCaa                                                                                        >  2:527942/1‑150 (MQ=255)
                                                                 cATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCaaaa                                                                                      <  2:630609/150‑1 (MQ=255)
                                                                  aacAGAAATCACCCCTGGTTAAGCAGACGCTATCGTTAACGTAACGCTTGATTACCGTTAATACTTTAGGGCTTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAAc                                                                                     <  1:1791470/148‑1 (MQ=255)
                                                                   tCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACg                                                                                    >  1:2662043/1‑150 (MQ=255)
                                                                   tCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCATTTACATCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAAAGATGTCaaaaag                                                                                    >  1:1117327/1‑148 (MQ=255)
                                                                    cAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTACCGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGa                                                                                   <  2:3329814/150‑1 (MQ=255)
                                                                    cAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGa                                                                                   >  2:3088820/1‑150 (MQ=255)
                                                                      gAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAAc                                                                                 >  1:1291437/1‑150 (MQ=255)
                                                                        aaTCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACAt                                                                               >  1:3079223/1‑150 (MQ=255)
                                                                         aTCACCCCAGGGGAAGCAGACTCTATCTTTAACGTAACGCTGTTTTACCTTTAATACTTTAGGGCTTTATCTGACGTTTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACAtt                                                                              <  1:1072536/150‑1 (MQ=255)
                                                                         aTCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACAtt                                                                              >  2:1019403/1‑150 (MQ=255)
                                                                         aTCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACAtt                                                                              <  2:1165706/150‑1 (MQ=255)
                                                                              cccAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCATATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCgg                                                                         <  1:2155085/150‑1 (MQ=255)
                                                                               ccAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGTATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGt                                                                        <  2:2815982/150‑1 (MQ=255)
                                                                                   gggAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGata                                                                    >  2:747177/1‑150 (MQ=255)
                                                                                    ggAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGatat                                                                   <  1:3886311/150‑1 (MQ=255)
                                                                                     gAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATAtt                                                                  >  1:301087/1‑150 (MQ=255)
                                                                                        gCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAAc                                                               <  1:976314/150‑1 (MQ=255)
                                                                                        gCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAAc                                                               <  1:3168866/150‑1 (MQ=255)
                                                                                         cAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAAcc                                                              <  1:2162395/150‑1 (MQ=255)
                                                                                           gACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCtg                                                            >  2:1505644/1‑149 (MQ=255)
                                                                                            aCGCTATCGTTAACGTAACGCTGGATTACCGGTAATACTTTAGGGCGTTATCTGTCCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTGGTATCAGGTTG‑CCAAATAACTATGTCAAAACGAACCAGCGCGGTGATATTAACCtgg                                                           >  2:2421400/1‑148 (MQ=255)
                                                                                              gCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgtt                                                         <  2:2214362/150‑1 (MQ=255)
                                                                                                tATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTtgttgt                                                       <  2:3914877/150‑1 (MQ=255)
                                                                                                  tCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTgtgt                                                     <  1:1469144/150‑1 (MQ=255)
                                                                                                   cGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTa                                                    <  1:3485318/150‑1 (MQ=255)
                                                                                                   cGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTa                                                    <  1:2888155/150‑1 (MQ=255)
                                                                                                    gTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTAACTGACCTGTCAGATAACGCCCATTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCATGGTGTGTAg                                                   >  1:525614/1‑150 (MQ=255)
                                                                                                        acgtaacgCTGTATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTc                                               <  1:750858/150‑1 (MQ=255)
                                                                                                         cgtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCt                                              <  2:3079223/150‑1 (MQ=255)
                                                                                                         cgtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCt                                              <  2:3774293/150‑1 (MQ=255)
                                                                                                          gtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTg                                             <  1:2843703/150‑1 (MQ=255)
                                                                                                               gCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTg                                        >  2:851637/1‑150 (MQ=255)
                                                                                                                cTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTTTCAGATAAAGCCCATTTCCATCCACTTTATCGTTTTATAAGGTGG‑AAAAATAAAAATGTAAAAAAGAACAGTCGCGGTGCTATTAACAGTGTTGTGTGGCGTATGTTTTgg                                       >  1:193710/1‑150 (MQ=255)
                                                                                                                   gATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCa                                    <  1:3725441/150‑1 (MQ=255)
                                                                                                                   gATTACAGTTAATACGTTAGGGCGTTATCTGACCTGTAAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCa                                    <  2:3250390/150‑1 (MQ=255)
                                                                                                                      tACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGgc                                 <  1:3232963/150‑1 (MQ=255)
                                                                                                                              aTACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCag                         <  2:1757087/150‑1 (MQ=255)
                                                                                                                               tACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCagg                        <  1:527942/150‑1 (MQ=255)
                                                                                                                                aCGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCCGGCGCTTGCagg                        >  2:2941450/1‑150 (MQ=255)
                                                                                                                                aCGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCagga                       <  2:2061775/150‑1 (MQ=255)
                                                                                                                                     aGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAAt                  <  1:354801/150‑1 (MQ=255)
                                                                                                                                        gCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGt               <  2:1177301/150‑1 (MQ=255)
                                                                                                                                        gCGTTATCTGACAGGTCCGATAACGCCCTTTTCCTTCCTCTATCTCGTTGTATCAGGTTG‑ACAAAAAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTAGTGTAGCTTCTGTCTAGGCCAGGCGCTTGCAGGAGGAATCGt               >  2:3827834/1‑150 (MQ=255)
                                                                                                                                          gTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGttt             <  2:3171891/150‑1 (MQ=255)
                                                                                                                                          gTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGttt             <  1:3446716/150‑1 (MQ=255)
                                                                                                                                          gTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGGATTGGCCAGGCGCTTGCAGGAGGACTCGttt             >  1:152857/1‑150 (MQ=255)
                                                                                                                                           ttATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGtttt            >  2:3284771/1‑150 (MQ=255)
                                                                                                                                           ttATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGtttt            >  2:3893605/1‑150 (MQ=255)
                                                                                                                                            tATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTa           <  1:2207262/150‑1 (MQ=255)
                                                                                                                                            tATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTa           <  1:3897377/150‑1 (MQ=255)
                                                                                                                                                 gACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCg      >  2:2539364/1‑150 (MQ=255)
                                                                                                                                                 gACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCg      >  1:2597441/1‑150 (MQ=255)
                                                                                                                                                 gACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCGTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCg      >  2:1619224/1‑150 (MQ=255)
                                                                                                                                                  aCCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGa     >  1:3841652/1‑150 (MQ=255)
                                                                                                                                                     tGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAAcg  >  2:3687115/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCCGGCGGAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTG‑AAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAACG  >  NC_000913/3287884‑3288182

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