Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,495,434 A→G 5.8% Q49Q (CAA→CAG ydcJ → DUF1338 domain‑containing protein YdcJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,495,4340AG 5.8% 210.5 / 2.5 69Q49Q (CAA→CAGydcJDUF1338 domain‑containing protein YdcJ
Reads supporting (aligned to +/- strand):  ref base A (35/30);  new base G (2/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CTATGGCGAACAGCATCACGGCGGATGAGATTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATC  >  NC_000913/1495286‑1495582
                                                                                                                                                    |                                                                                                                                                    
cTATGGCGAACAGCATCACGGCGGATGAGATTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAAc                                                                                                                                                     >  1:1633303/1‑150 (MQ=255)
  aTGGCGAACAGCATCACGGCGGATGAGATTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTg                                                                                                                                                   >  1:438182/1‑150 (MQ=255)
            gCATCACGGCGGATGAGATTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGaaaaaa                                                                                                                                         >  2:1277357/1‑150 (MQ=255)
            gCATAACGGCGGATGAGATTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGaaaaaa                                                                                                                                         >  2:533525/1‑150 (MQ=255)
             cATCTCGGCGGATGAGATTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACg                                                                                                                                               <  2:1419153/143‑1 (MQ=255)
             cATCACGGCGGATGAGATTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACg                                                                                                                                               >  1:1419153/1‑143 (MQ=255)
               tCACGGCGGATGAGATTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATgg                                                                                                                                      <  2:11647/150‑1 (MQ=255)
                cACGGCGGATGAGATTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGt                                                                                                                                     <  1:71004/150‑1 (MQ=255)
                 aCGGCGGATGAGATTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTa                                                                                                                                    <  1:1762541/150‑1 (MQ=255)
                  cggcggATGAGATTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTaa                                                                                                                                   >  1:703134/1‑150 (MQ=255)
                      ggATGAGATTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGc                                                                                                                               <  1:1277357/150‑1 (MQ=255)
                        aTGAGATTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAg                                                                                                                             <  2:1549850/150‑1 (MQ=255)
                         tGAGATTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGa                                                                                                                            >  1:911264/1‑150 (MQ=255)
                            gaTTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGa                                                                                                                         >  2:674069/1‑150 (MQ=255)
                                cGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTg                                                                                                                     >  2:358161/1‑150 (MQ=255)
                                  ggAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGgc                                                                                                                   >  2:1757716/1‑150 (MQ=255)
                                       aGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGAc                                                                                                              <  1:1656978/150‑1 (MQ=255)
                                        gTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACt                                                                                                             <  2:758903/150‑1 (MQ=255)
                                            tCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAAt                                                                                                         <  1:706550/150‑1 (MQ=255)
                                              gCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGt                                                                                                       >  2:1665450/1‑150 (MQ=255)
                                              gCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGt                                                                                                       <  1:533525/150‑1 (MQ=255)
                                                     aTGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCa                                                                                                                              <  1:1574492/120‑1 (MQ=255)
                                                     aTGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCa                                                                                                                              >  2:1574492/1‑120 (MQ=255)
                                                                 tACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGAtt                                                                                    >  1:1196471/1‑150 (MQ=255)
                                                                   ccAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAACACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTcg                                                                                  >  2:1610841/1‑150 (MQ=255)
                                                                     aGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTcgcg                                                                                >  1:1840165/1‑150 (MQ=255)
                                                                       cAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGtt                                                                              <  1:358161/150‑1 (MQ=255)
                                                                         agaagTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTgg                                                                            >  1:322405/1‑150 (MQ=255)
                                                                          gaagTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCATCAACTGCACGAAAAAATGTTAAATTCAGACGAGATGGCGAGACTGAATGTTGAACTTCATGGGGCGATTCGCGTTggg                                                                           >  2:265292/1‑150 (MQ=255)
                                                                            agTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGAc                                                                         >  2:1191709/1‑150 (MQ=255)
                                                                             gTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATg                                                                                                        >  1:261413/1‑118 (MQ=255)
                                                                             gTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATg                                                                                                        <  2:261413/118‑1 (MQ=255)
                                                                             gTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACt                                                                        <  2:1391266/150‑1 (MQ=255)
                                                                              ttCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTg                                                                       <  1:1279648/150‑1 (MQ=255)
                                                                               tCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGc                                                                      >  1:141230/1‑150 (MQ=255)
                                                                                   cAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACaa                                                                  >  1:1473821/1‑150 (MQ=255)
                                                                                    aaTATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAag                                                                 >  1:275882/1‑150 (MQ=255)
                                                                                       atGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGc                                                              >  2:1532562/1‑150 (MQ=255)
                                                                                         ggCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCtt                                                            >  1:799285/1‑150 (MQ=255)
                                                                                           cACGCTGCTGGAACTGGTAGCTGATGGGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGc                                                          <  2:1698364/150‑1 (MQ=255)
                                                                                                  ctgGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAACACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCtt                                                   <  1:1610841/150‑1 (MQ=255)
                                                                                                   tgGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATCCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTc                                                  <  2:438182/150‑1 (MQ=255)
                                                                                                    ggAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTcg                                                 >  2:1004202/1‑150 (MQ=255)
                                                                                                    ggAACTGGTAGCTGATGTGAATCTGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTcgg                                                >  2:1669507/1‑150 (MQ=255)
                                                                                                        cTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTcggcgg                                             <  1:1858121/150‑1 (MQ=255)
                                                                                                        cTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTcggcgg                                             <  1:674069/150‑1 (MQ=255)
                                                                                                           gTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATg                                          >  1:119139/1‑150 (MQ=255)
                                                                                                           gTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGCTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATg                                          <  2:16333/150‑1 (MQ=255)
                                                                                                            tAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGt                                         >  2:1102244/1‑150 (MQ=255)
                                                                                                                    gtgAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCAtt                                 >  1:94845/1‑150 (MQ=255)
                                                                                                                     tgAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTa                                <  1:1004202/150‑1 (MQ=255)
                                                                                                                       aaTCTGGCTGTGCTGGAAAACAATCCTCAGCTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATg                              >  1:1438239/1‑150 (MQ=255)
                                                                                                                       aaTCTGGCTGTGCTGGAAAACAATCCTCAGCTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATg                              >  1:1445243/1‑150 (MQ=255)
                                                                                                                          cTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATgggg                           <  2:1633303/150‑1 (MQ=255)
                                                                                                                             gCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATg                        <  2:275882/150‑1 (MQ=255)
                                                                                                                              cTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGt                       >  2:1293283/1‑150 (MQ=255)
                                                                                                                                gtgCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGc                                                          <  2:859819/113‑1 (MQ=255)
                                                                                                                                gtgCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGc                                                          >  1:859819/1‑113 (MQ=255)
                                                                                                                                gtgCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGcc                                    >  2:119974/1‑135 (MQ=255)
                                                                                                                                gtgCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGcc                                    <  1:119974/135‑1 (MQ=255)
                                                                                                                                   cTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCg                  >  2:1778983/1‑150 (MQ=255)
                                                                                                                                     ggAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGt                >  2:222351/1‑150 (MQ=255)
                                                                                                                                       aaaaCAATCCTCAGCTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGa              <  2:1445243/150‑1 (MQ=255)
                                                                                                                                       aaaaCAATCCTCAGCTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGa              <  2:1438239/150‑1 (MQ=255)
                                                                                                                                        aaaCAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAg             <  1:1192255/150‑1 (MQ=255)
                                                                                                                                        aaaCAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAg             >  2:1649115/1‑150 (MQ=255)
                                                                                                                                              tCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTa       <  2:322405/150‑1 (MQ=255)
                                                                                                                                              tCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTa       <  1:1191709/150‑1 (MQ=255)
                                                                                                                                                  cAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAt   <  2:799285/150‑1 (MQ=255)
                                                                                                                                                   aaCTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtc  <  1:222351/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTATGGCGAACAGCATCACGGCGGATGAGATTCGGGAACAGTTTTCGCAGGCAATGTCAGCCATGTACCAGCAAGAAGTTCCGCAATATGGCACGCTGCTGGAACTGGTAGCTGATGTGAATCTGGCTGTGCTGGAAAACAATCCTCAACTGCACGAAAAAATGGTAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATC  >  NC_000913/1495286‑1495582

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: