Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,485,806 Δ1 bp coding (146/1512 nt) pepA ← multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,485,8030C.98.5% 316.5 / ‑3.4 68coding (149/1512 nt)pepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (39/28);  minor base A (1/0);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTC  >  NC_000913/4485653‑4485948
                                                                                                                                                      |                                                                                                                                                 
aGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGc                                                                                                                                                   >  1:71022/1‑151 (MQ=255)
 gCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGcc                                                                                                                                                  <  2:368264/151‑1 (MQ=255)
  cGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGccc                                                                                                                                                 >  1:100610/1‑151 (MQ=255)
   gTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGcccg                                                                                                                                                >  1:338952/1‑150 (MQ=255)
   gTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGcccg                                                                                                                                                <  2:411078/151‑2 (MQ=255)
     aTTAAAGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCTCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGcccgac                                                                                                                                              <  2:733981/151‑4 (MQ=255)
      ttAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGCGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGcccgacg                                                                                                                                             >  1:300790/1‑147 (MQ=255)
          tGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGc                                                                                                                                        >  1:863749/1‑151 (MQ=255)
          tGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGc                                                                                                                                        >  1:634410/1‑151 (MQ=255)
              tttCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGAAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAgg                                                                                                                                     <  1:861151/150‑1 (MQ=255)
                 cTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGgcgc                                                                                                                                 >  1:104814/1‑151 (MQ=255)
                  tGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCt                                                                                                                                >  2:614606/1‑151 (MQ=255)
                   gAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCt                                                                                                                                <  2:645568/150‑1 (MQ=255)
                    aaTAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTg                                                                                                                               >  1:93854/1‑150 (MQ=255)
                        aCCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTTCCCCGGTTTTCCTTCCCGTTCG‑CCCGACGTAGCAGGGCGCTGAt                                                                                                                             <  2:437185/148‑1 (MQ=255)
                        aCCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGt                                                                                                                           >  2:229097/1‑150 (MQ=255)
                         ccTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTAc                                                                                                                         >  2:590780/1‑151 (MQ=255)
                         ccTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTC‑ACCCGACGTAGCAGGGCGCTGATGTAc                                                                                                                         >  1:322570/1‑151 (MQ=255)
                            gCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTAccc                                                                                                                       >  1:573283/1‑150 (MQ=255)
                            gCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTAcc                                                                                                                        <  1:411092/149‑1 (MQ=255)
                               tGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATc                                                                                                                    <  1:106121/150‑1 (MQ=255)
                                gTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCg                                                                                                                   <  1:13539/150‑1 (MQ=255)
                                  aCTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCt                                                                                                                 <  2:235629/150‑1 (MQ=255)
                                  aCTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTg                                                                                                                >  1:212738/1‑151 (MQ=255)
                                   cTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTg                                                                                                                >  1:383986/1‑150 (MQ=255)
                                         cTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTAt                                                                                                         >  2:187654/1‑151 (MQ=255)
                                          tCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATc                                                                                                        >  2:818929/1‑151 (MQ=255)
                                              ccAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAg                                                                                                     <  1:785551/150‑1 (MQ=255)
                                                   tCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTc                                                                                               >  1:66271/1‑151 (MQ=255)
                                                    cACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCt                                                                                              >  2:31900/1‑151 (MQ=255)
                                                     aCGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTg                                                                                             <  2:283173/151‑1 (MQ=255)
                                                      cGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGc                                                                                            <  1:88149/151‑1 (MQ=255)
                                                          cTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATc                                                                                        >  2:875262/1‑151 (MQ=255)
                                                            ttGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATc                                                                                        <  1:391879/149‑1 (MQ=255)
                                                              gCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGga                                                                                     <  1:11754/150‑1 (MQ=255)
                                                                cGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGgaga                                                                                   >  2:275667/1‑150 (MQ=255)
                                                                   aaCCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAgg                                                                               <  1:496218/151‑1 (MQ=255)
                                                                    aCCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGc                                                                              <  1:32046/151‑1 (MQ=255)
                                                                     ccAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAgg                                                                               <  2:552588/149‑1 (MQ=255)
                                                                     ccAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCg                                                                             <  1:495260/151‑1 (MQ=255)
                                                                      cAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCg                                                                             >  1:614988/1‑150 (MQ=255)
                                                                       aaTAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGAc                                                                           <  1:777354/151‑1 (MQ=255)
                                                                         tAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAGTCGGAGAAAGGCGACGt                                                                         >  2:447951/1‑151 (MQ=255)
                                                                             gagaGTTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAGCAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGtt                                                                     <  2:38289/151‑1 (MQ=255)
                                                                              agaATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTc                                                                    >  1:528458/1‑151 (MQ=255)
                                                                                  ttCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGa                                                                  >  1:115202/1‑149 (MQ=255)
                                                                                            aaGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacg                                                       >  1:824824/1‑150 (MQ=255)
                                                                                             aGTACATTCGGGACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacgac                                                     >  1:24250/1‑151 (MQ=255)
                                                                                             aGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacgac                                                     <  2:291340/151‑1 (MQ=255)
                                                                                             aGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacgac                                                     >  1:1279/1‑151 (MQ=255)
                                                                                                aCATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGAt                                                  >  2:343573/1‑151 (MQ=255)
                                                                                                 cATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATg                                                 <  2:909302/151‑1 (MQ=255)
                                                                                                  aTTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGc                                                >  2:587625/1‑151 (MQ=255)
                                                                                                      ggAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAg                                              <  2:881734/149‑1 (MQ=255)
                                                                                                      ggAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGc                                            >  1:352385/1‑151 (MQ=255)
                                                                                                      ggAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGc                                            >  1:871413/1‑151 (MQ=255)
                                                                                                       gAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGc                                            >  2:247906/1‑150 (MQ=255)
                                                                                                       gAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGc                                            >  2:760960/1‑150 (MQ=255)
                                                                                                           aTGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTcc                                       <  2:824757/151‑1 (MQ=255)
                                                                                                              gTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCt                                    <  2:93621/151‑1 (MQ=255)
                                                                                                                  agcaACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGttt                                >  2:332833/1‑151 (MQ=255)
                                                                                                                   gcaACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTtc                               >  1:734382/1‑151 (MQ=255)
                                                                                                                     aacaaTGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTtctc                             >  1:180317/1‑151 (MQ=255)
                                                                                                                     aacaaTGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTtctc                             <  2:71418/151‑1 (MQ=255)
                                                                                                                     aacaaTGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTtctc                             >  1:67939/1‑151 (MQ=255)
                                                                                                                      acaaTGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTcc                            <  1:289676/151‑1 (MQ=255)
                                                                                                                         aTGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCggg                         <  2:573233/151‑1 (MQ=255)
                                                                                                                           gTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGc                        <  2:644202/150‑1 (MQ=255)
                                                                                                                             cTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTAc                     >  2:156844/1‑151 (MQ=255)
                                                                                                                               gCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCg                   >  2:825351/1‑151 (MQ=255)
                                                                                                                                   cGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCttt               <  2:512624/151‑1 (MQ=255)
                                                                                                                                       tttCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTac           >  1:61735/1‑151 (MQ=255)
                                                                                                                                           cttccAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt        >  2:686076/1‑150 (MQ=255)
                                                                                                                                              ccAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc    >  2:464051/1‑151 (MQ=255)
                                                                                                                                                aGTTCG‑CCCGACGTAGCAGGGCGCTGGTGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc  >  1:159119/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
AGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTC  >  NC_000913/4485653‑4485948

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑