Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,041,221 G→A 6.3% A71A (GCC→GCT mobB ← molybdopterin‑guanine dinucleotide biosynthesis protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,041,2210GA 6.3% 341.9 / 7.2 111A71A (GCC→GCTmobBmolybdopterin‑guanine dinucleotide biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base G (57/46);  new base A (5/2);  total (63/48)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

Only 100 of 113 total aligned reads displayed.

ATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAA  >  NC_000913/4041081‑4041366
                                                                                                                                            |                                                                                                                                                 
aTCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGAt                                                                                                                                           >  2:1222842/1‑149 (MQ=255)
 tCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATg                                                                                                                                          <  1:712050/149‑1 (MQ=255)
  cTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATgg                                                                                                                                         <  2:508840/149‑1 (MQ=255)
       cGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTtgt                                                                                                                                    <  1:1341764/149‑1 (MQ=255)
        gagctctTCATGCTTAAACCCTTCGCCCAGAATCAAAACCAGCTTTGAGGTATCCATTCGCCTTGCGAGAAACTGTAGATGCAGCAGTTCTTGGTCTGGTGTGTCCGGCATCAAGGCCCATCGTTGCTGGCTGGCAACGAGGGGTtggg                                                                                                                                   >  2:54331/4‑147 (MQ=255)
        gATCTCTTCATGCTTAAACCCTTCGCCCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGATACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGGCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTtgtg                                                                                                                                   >  2:270993/1‑149 (MQ=255)
           ctctTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGGCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCcg                                                                                                                                >  1:1576998/1‑149 (MQ=255)
              ttCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCgcgc                                                                                                                             >  2:705299/1‑149 (MQ=255)
              ttCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCgcgc                                                                                                                             >  1:1397434/1‑149 (MQ=255)
              ttCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCgcgc                                                                                                                             >  1:976008/1‑149 (MQ=255)
              ttCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCgcgc                                                                                                                             >  1:807849/1‑149 (MQ=255)
               tCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGcc                                                                                                                            >  2:1502943/1‑149 (MQ=255)
                  tGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGc                                                                                                                         <  1:526091/149‑1 (MQ=255)
                  tGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACATGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCTTTGCTGGCTGCCAACGATGGTTATTGCCGCGCCAGc                                                                                                                         >  1:819555/1‑149 (MQ=255)
                   gCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGcc                                                                                                                        >  2:1125412/1‑149 (MQ=255)
                    cTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCt                                                                                                                       >  1:1335525/1‑149 (MQ=255)
                    cTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCt                                                                                                                       >  1:598490/1‑149 (MQ=255)
                       aaaCCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGACTGGTGTTTCCGTCATCAAGGCCCATCGTTTCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCttcc                                                                                                                    >  2:456879/1‑147 (MQ=255)
                           ccTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAg                                                                                                                >  1:439480/1‑149 (MQ=255)
                            cTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGc                                                                                                               <  1:1125412/149‑1 (MQ=255)
                            cTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGc                                                                                                               <  1:1310917/149‑1 (MQ=255)
                            attCGACCAGAATCAAATCCAGCTTTGAGGTATAAATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTACGTCATCAAGGCCCATCGTAGCTGGCTGGCAACGATGGTTTTTGCAGCGCCAGCCTTTCGCAGc                                                                                                               >  2:346273/2‑149 (MQ=255)
                                 aCCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATa                                                                                                          >  2:1322151/1‑149 (MQ=255)
                                  ccAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAg                                                                                                         <  1:51572/149‑1 (MQ=255)
                                  ccAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCGCTTCTTCGTCTGGGGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAg                                                                                                         <  2:1048481/149‑1 (MQ=255)
                                   cAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACTATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGc                                                                                                        >  2:291462/1‑149 (MQ=255)
                                   cAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGc                                                                                                        <  2:1472026/149‑1 (MQ=255)
                                    aGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCt                                                                                                       <  2:1369887/149‑1 (MQ=255)
                                     gAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTa                                                                                                      >  2:931412/1‑149 (MQ=255)
                                      aaTCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTATTCTTCGTCTGGTGTTTACGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGTTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAACTAt                                                                                                     >  2:308807/1‑149 (MQ=255)
                                        tCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCt                                                                                                   >  1:1062868/1‑149 (MQ=255)
                                        tCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCACATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCt                                                                                                   >  2:490482/1‑149 (MQ=255)
                                            aTCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGc                                                                                               <  2:55615/149‑1 (MQ=255)
                                              ccAGCGTGGAGGGATCCATGCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCGGGTGTTGCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCt                                                                                             <  1:273591/149‑1 (MQ=255)
                                               cAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTg                                                                                            >  1:193250/1‑149 (MQ=255)
                                                aGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTgg                                                                                           >  1:205178/1‑149 (MQ=255)
                                                 gCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGc                                                                                          >  2:615628/1‑149 (MQ=255)
                                                   tttGAGGTATCCATTAGACTTGCGAGAAACTGTAGATACAGCTCTTCTTCGTCTGGTGTTTCCGTCATAAAGGCCCATAGTTGCTGGCTTGCAAAGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGATATCTTTTCCTGGCtt                                                                                        >  2:464856/1‑149 (MQ=255)
                                                     tGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTAt                                                                                      <  2:246941/149‑1 (MQ=255)
                                                      gAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATc                                                                                     <  2:1208191/149‑1 (MQ=255)
                                                      gAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCGGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCGGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATc                                                                                     <  1:1427558/149‑1 (MQ=255)
                                                         gTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCCTCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCacc                                                                                  >  1:219363/1‑147 (MQ=255)
                                                         gTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAAc                                                                                  >  2:1296663/1‑149 (MQ=255)
                                                            tCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATc                                                                               >  2:88439/1‑149 (MQ=255)
                                                             ccATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCATGGCTTATCAACATcc                                                                              >  1:1145175/1‑149 (MQ=255)
                                                               aTTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCat                                                                            <  1:543287/149‑1 (MQ=255)
                                                               aTTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCat                                                                            >  2:454018/1‑149 (MQ=255)
                                                                ttCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGATATCTTTGCCTGGCTTATCAACATCCata                                                                           >  2:671097/1‑149 (MQ=255)
                                                                  cGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATc                                                                         <  1:1502943/149‑1 (MQ=255)
                                                                   gACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCa                                                                        >  1:1318321/1‑149 (MQ=255)
                                                                     cTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATg                                                                      <  2:1062868/149‑1 (MQ=255)
                                                                     cTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCATGGCTTATAAACATCCATATCATg                                                                      >  2:968792/1‑149 (MQ=255)
                                                                     cTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCAGAGCTATCTTTGCCTGGGTTATCAACATCCATATCATg                                                                      >  2:265905/1‑149 (MQ=255)
                                                                       tGCGAGAAACTGTAGATCAAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATATCTATCTTTGCCTGGCTTATCAACATCCATATCATGGt                                                                    <  1:346273/149‑1 (MQ=255)
                                                                       tGCGAGAAACGGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCAGCAAGGCCCATCGTTGCTGGCGGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGt                                                                    <  1:944840/149‑1 (MQ=255)
                                                                          gagaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCg                                                                 >  1:399235/1‑149 (MQ=255)
                                                                           agaAACTGTAGATCCAGCTCTTCTTCGGCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCACCATCCATATCATGGTGCGt                                                                >  2:1241721/1‑149 (MQ=255)
                                                                            taaaCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTa                                                               <  2:1251454/148‑1 (MQ=255)
                                                                            gaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTa                                                               <  1:805387/149‑1 (MQ=255)
                                                                             aaaCTGTAGATCCAGCTCTTCATCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTAt                                                              <  2:1160336/149‑1 (MQ=255)
                                                                              aaCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATAGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGTTGCGTATg                                                             >  1:559701/1‑149 (MQ=255)
                                                                               aCTGTAGATCCAGCGCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCAGCGCCAGCCTTTCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTAttc                                                            <  2:48429/149‑3 (MQ=255)
                                                                                 tGTAGATCCAGATCTTCTTCGTCTGGTTTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATTGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCtt                                                          <  1:456879/149‑1 (MQ=255)
                                                                                  ttAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTa                                                         <  1:884909/148‑1 (MQ=255)
                                                                                  gTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTa                                                         >  2:575911/1‑149 (MQ=255)
                                                                                  gTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTa                                                         >  2:879502/1‑149 (MQ=255)
                                                                                  gTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTa                                                         <  1:270993/149‑1 (MQ=255)
                                                                                      aTCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCa                                                     >  2:49103/1‑149 (MQ=255)
                                                                                       tCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAg                                                    <  2:976008/149‑1 (MQ=255)
                                                                                           gCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCt                                                >  1:624698/1‑149 (MQ=255)
                                                                                               ttcttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGAc                                            >  1:1593750/1‑149 (MQ=255)
                                                                                                 cttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                            >  1:1134008/1‑131 (MQ=255)
                                                                                                  ttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGa                                         >  2:380775/1‑149 (MQ=255)
                                                                                                  ttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGACAGCCTTTCGCAGCTAAAACATATCTTTGCCTGGCTTATCAACATCCATATAAAGGTGCGTATGCTTAATCAGCCCTGTACGGa                                         >  1:972114/1‑149 (MQ=255)
                                                                                                    cGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATc                                       <  2:60954/149‑1 (MQ=255)
                                                                                                    cGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATc                                       <  2:1335525/149‑1 (MQ=255)
                                                                                                       cTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCct                                    <  1:705299/149‑1 (MQ=255)
                                                                                                          gtgtTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTg                                 <  2:1515490/149‑1 (MQ=255)
                                                                                                            gtTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGgc                               <  1:401607/149‑1 (MQ=255)
                                                                                                             tttCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGgcg                              >  2:1483023/1‑149 (MQ=255)
                                                                                                                   tcatcaAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCGCATAGCTATCTTTGCCTGGCTTTTCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAAt                        <  1:971301/149‑1 (MQ=255)
                                                                                                                      tcaAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTGGCGCAGCTCATAGCTATCTTTGCCTGGCTTTTCAACATCCATATCATGGTGCGTATGCTTATTCAGCCCTGGTCGGATCCCTCTGGCGCATAATggc                     >  1:252828/1‑147 (MQ=255)
                                                                                                                       caAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCg                    <  1:264514/149‑1 (MQ=255)
                                                                                                                       caAGGCCCATCGTTGCTGGCTGGAAACGATGGTTTGTGCCGCGCGAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATTGTGCGTACGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCg                    <  2:819555/149‑1 (MQ=255)
                                                                                                                        aaGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCgg                   <  2:652391/149‑1 (MQ=255)
                                                                                                                           gCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATcc                                      <  1:1322443/127‑1 (MQ=255)
                                                                                                                            cccATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATAAATATCATGGTGCGTATGCTTAATCAGACCTGGACGGATCCCTCTGGCGCAGAATGCCGGGATc               >  1:1353929/1‑149 (MQ=255)
                                                                                                                             ccATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCa              <  1:437549/149‑1 (MQ=255)
                                                                                                                             ccATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCa              >  2:1370880/1‑149 (MQ=255)
                                                                                                                              cATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAg             >  1:1564327/1‑149 (MQ=255)
                                                                                                                              cATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCGCTGGACGGATCCCTCTGGCGCATAATGCAGGGATCat             >  1:415322/1‑148 (MQ=255)
                                                                                                                              cATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAg             >  2:891381/1‑149 (MQ=255)
                                                                                                                              cATCGTTGCCAGCCAGCAACGATGGTATGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTAAGCTTTATCAGCCCTGGACGGATCCCTCTGTCGCATAATGCCGGGCTCAg             >  2:151189/1‑149 (MQ=255)
                                                                                                                               aTCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGt            <  1:668774/149‑1 (MQ=255)
                                                                                                                                  gTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGtttt         <  2:205178/149‑1 (MQ=255)
                                                                                                                                  gTTGCTGGCTGGCAACGATGGTTTGGTCCGCGCCAGCCTTGCGCAGCTCATATCTATCTTTGCCTGTCTTATCCACATCCATATCATGGTGCGTATGCTTAATCAGAACTGGACGGATCACTCTGGCGCATAATGCCGGGATCAGtttt         <  1:869335/149‑1 (MQ=255)
                                                                                                                                    tGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGtttttt       <  1:583091/149‑1 (MQ=255)
                                                                                                                                    tGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGtttttt       <  1:746110/149‑1 (MQ=255)
                                                                                                                                    tGCTGGCTGGCAACGATGGGTTGTGGCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTAAGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGtttttt       <  2:127492/149‑1 (MQ=255)
                                                                                                                                         gctggcAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCaaaa  <  1:1222842/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
ATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAA  >  NC_000913/4041081‑4041366

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