Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,041,217 G→A 7.2% Q73* (CAG→TAG) ‡ mobB ← molybdopterin‑guanine dinucleotide biosynthesis adaptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,041,2170GA 7.2% 430.6 / 14.6 138Q73* (CAG→TAG) ‡mobBmolybdopterin‑guanine dinucleotide biosynthesis adaptor protein
Reads supporting (aligned to +/- strand):  ref base G (64/64);  new base A (4/6);  total (68/70)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 152 total aligned reads displayed.

CGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAA  >  NC_000913/4041069‑4041366
                                                                                                                                                    |                                                                                                                                                     
cgAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTgctgg                                                                                                                                                      >  1:887351/1‑150 (MQ=255)
cgAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTgctgg                                                                                                                                                      <  1:1357019/150‑1 (MQ=255)
  aaaCAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCTTCATGGCCTATCGTTgctggct                                                                                                                                                    <  2:1429777/150‑1 (MQ=255)
    aCAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGctggctgg                                                                                                                                                  <  1:303818/150‑1 (MQ=255)
      aGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCGTTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCa                                                                                                                                                >  1:1101314/1‑150 (MQ=255)
        cacaATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAAc                                                                                                                                              <  2:301895/150‑1 (MQ=255)
          caATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAACTCCAGCTTTGAGGGATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGa                                                                                                                                            <  1:503846/150‑1 (MQ=255)
           aaTCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGAt                                                                                                                                           <  2:181724/150‑1 (MQ=255)
            aTCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATg                                                                                                                                          <  1:478066/150‑1 (MQ=255)
              cTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATggg                                                                                                                                        <  1:1221176/150‑2 (MQ=255)
               tttGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGtt                                                                                                                                       >  1:1093186/1‑150 (MQ=255)
               tttGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATgggc                                                                                                                                       <  2:578825/150‑3 (MQ=255)
               tttGCGATCTATTCATGCTTAAACCCTTAGACCAGAATCAAATGCAGCTTTGAGGTATACATTCGACTTTCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCGATAGTTGCCATCCAGCAACGATgggc                                                                                                                                       <  2:65584/150‑3 (MQ=255)
                    gATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGGTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGc                                                                                                                                  <  2:321539/150‑1 (MQ=255)
                    gATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCtggctggc                                                                                                                                                 >  2:1234860/1‑135 (MQ=255)
                        tctTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCcgcg                                                                                                                              >  1:524610/1‑150 (MQ=255)
                         ctTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCgcgc                                                                                                                             >  1:1049555/1‑150 (MQ=255)
                              tGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGcc                                                                                                                        >  2:31870/1‑150 (MQ=255)
                                  tAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgc                                                                                                                    <  2:41926/150‑1 (MQ=255)
                                    aaCCGTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc                                                                                                                  <  1:638363/150‑1 (MQ=255)
                                      cccTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAg                                                                                                                >  1:805183/1‑150 (MQ=255)
                                        cTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCt                                                                                                              >  2:932362/1‑150 (MQ=255)
                                         ttCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTc                                                                                                             <  1:692227/150‑1 (MQ=255)
                                           cGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCAt                                                                                                           <  2:1146517/150‑1 (MQ=255)
                                              ccAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGc                                                                                                        >  1:288229/1‑150 (MQ=255)
                                                aGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTa                                                                                                      <  1:879547/150‑1 (MQ=255)
                                                 gAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAg                                                                                                         <  2:1349979/146‑1 (MQ=255)
                                                 gAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAg                                                                                                         >  1:1349979/1‑146 (MQ=255)
                                                 gAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGCCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATACCTAt                                                                                                     >  2:994980/1‑150 (MQ=255)
                                                   aTCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCt                                                                                                   <  2:102162/150‑1 (MQ=255)
                                                      aaaTCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTg                                                                                                >  1:65874/1‑150 (MQ=255)
                                                            aGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTgctggctg                                                                                                                                                   >  2:67748/1‑93 (MQ=255)
                                                            aGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGc                                                                                          >  1:237717/1‑150 (MQ=255)
                                                                  gAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCa                                                                                    >  1:687674/1‑150 (MQ=255)
                                                                    ggTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAAc                                                                                  >  1:387346/1‑150 (MQ=255)
                                                                     gTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACa                                                                                 >  2:1030923/1‑150 (MQ=255)
                                                                     gTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACa                                                                                 <  2:1424459/150‑1 (MQ=255)
                                                                      tATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACAt                                                                                >  1:999942/1‑150 (MQ=255)
                                                                         ccATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCa                                                                             <  2:900841/150‑1 (MQ=255)
                                                                          cATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCat                                                                            >  1:1152504/1‑150 (MQ=255)
                                                                          cATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCat                                                                            >  2:67807/1‑150 (MQ=255)
                                                                              cGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCa                                                                        >  1:761068/1‑150 (MQ=255)
                                                                                 cTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATgg                                                                     >  2:1157482/1‑150 (MQ=255)
                                                                                 cTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATgg                                                                     <  2:1360057/150‑1 (MQ=255)
                                                                                     cGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGCGGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATCTCATGGTGCg                                                                 >  2:970637/1‑150 (MQ=255)
                                                                                     cGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCg                                                                 >  1:302125/1‑150 (MQ=255)
                                                                                     cGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCg                                                                 <  1:618964/150‑1 (MQ=255)
                                                                                     cGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCg                                                                 <  2:887351/150‑1 (MQ=255)
                                                                                      gagaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGt                                                                >  1:422315/1‑150 (MQ=255)
                                                                                             tGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                            <  1:911743/147‑1 (MQ=255)
                                                                                             tGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                            >  2:911743/1‑147 (MQ=255)
                                                                                             tGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTa                                                         <  2:1101314/150‑1 (MQ=255)
                                                                                              gTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTaa                                                        >  2:437859/1‑150 (MQ=255)
                                                                                              ctAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTaa                                                        >  1:349812/2‑150 (MQ=255)
                                                                                                aGATCGAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGGCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCTTATGCTTAATc                                                      <  2:1184905/150‑1 (MQ=255)
                                                                                                aGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATc                                                      <  1:999816/150‑1 (MQ=255)
                                                                                                 gATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCa                                                     >  2:398741/1‑150 (MQ=255)
                                                                                                 gATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCa                                                     >  2:1023536/1‑150 (MQ=255)
                                                                                                   tCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGc                                                   <  2:792051/150‑1 (MQ=255)
                                                                                                   tCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGc                                                   <  1:899796/150‑1 (MQ=255)
                                                                                                       gCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTg                                               <  2:1049555/150‑1 (MQ=255)
                                                                                                          cttcttCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGAc                                            >  1:153952/1‑150 (MQ=255)
                                                                                                              ttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGAt                                        <  1:169839/150‑1 (MQ=255)
                                                                                                              ttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGAt                                        <  2:805183/150‑1 (MQ=255)
                                                                                                               tcGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTTGCAACGATGGTTTTTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTTGACGGATc                                       <  1:491202/150‑1 (MQ=255)
                                                                                                               tcGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATc                                       >  2:744543/1‑150 (MQ=255)
                                                                                                                 gTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCtt                                                          >  1:241342/1‑129 (MQ=255)
                                                                                                                 gTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATccc                                     >  1:1121733/1‑150 (MQ=255)
                                                                                                                 gTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATccc                                     >  2:516096/1‑150 (MQ=255)
                                                                                                                 gTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATccc                                     >  1:1385218/1‑150 (MQ=255)
                                                                                                                    tGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATc                                       >  2:229096/1‑145 (MQ=255)
                                                                                                                    tGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCctct                                  <  2:372047/150‑1 (MQ=255)
                                                                                                                      gtgtTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTgg                                <  1:1023634/150‑1 (MQ=255)
                                                                                                                      gtgtTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTgg                                >  1:652820/1‑150 (MQ=255)
                                                                                                                        gtTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGgcg                              >  1:648837/1‑150 (MQ=255)
                                                                                                                          ttCCGTCATCAAGGGCGATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCTGGCTTGCGCAGCTGATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATGAGCCCTGGACGGATCCCTCTGGCGCa                            <  2:433221/150‑1 (MQ=255)
                                                                                                                          ttCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCa                            >  1:937267/1‑150 (MQ=255)
                                                                                                                           tCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCAt                           <  1:676615/150‑1 (MQ=255)
                                                                                                                            ccGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATa                          >  1:1425086/1‑150 (MQ=255)
                                                                                                                                 atcaAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCCGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGcc                     >  2:649412/1‑150 (MQ=255)
                                                                                                                                 atcaAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGcc                     >  1:711282/1‑150 (MQ=255)
                                                                                                                                  tcaAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCg                    >  1:1272070/1‑150 (MQ=255)
                                                                                                                                  tcaAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCg                    >  1:1297531/1‑150 (MQ=255)
                                                                                                                                   caAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCgg                   <  2:292810/150‑1 (MQ=255)
                                                                                                                                   caAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCgg                   >  2:276389/1‑150 (MQ=255)
                                                                                                                                    aaGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCggg                  <  1:725075/150‑1 (MQ=255)
                                                                                                                                       gCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATc               >  1:991949/1‑150 (MQ=255)
                                                                                                                                         ccATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAg             >  1:823747/1‑150 (MQ=255)
                                                                                                                                         ccATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAg             >  1:373434/1‑150 (MQ=255)
                                                                                                                                         ccATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAg             >  1:1105222/1‑150 (MQ=255)
                                                                                                                                             cGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGtttt         >  1:323853/1‑150 (MQ=255)
                                                                                                                                             cGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGtttt         >  1:459972/1‑150 (MQ=255)
                                                                                                                                              gTTGCTGGCTGGCAACGATGGTTTGTGCCGGGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGttttt        <  1:994980/150‑1 (MQ=255)
                                                                                                                                              gTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGttttt        >  1:769623/1‑150 (MQ=255)
                                                                                                                                                tGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTc      <  2:1121733/150‑1 (MQ=255)
                                                                                                                                                tGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTc      <  1:1093839/150‑1 (MQ=255)
                                                                                                                                                 gctggctgGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTct     <  2:387766/150‑2 (MQ=255)
                                                                                                                                                 gctggctgGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCa     >  1:944665/1‑150 (MQ=255)
                                                                                                                                                   tggctggcAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCaaa   >  1:869100/1‑150 (MQ=255)
                                                                                                                                                    ggctggcAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCaaaa  >  1:1177890/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAA  >  NC_000913/4041069‑4041366

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: