Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 665,094 G→C 5.9% P33A (CCT→GCT)  rlpA ← septal ring protein, suppressor of prc, minor lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913665,0940GC 5.9% 267.2 / 6.9 85P33A (CCT→GCT) rlpAseptal ring protein, suppressor of prc, minor lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (30/50);  new base C (2/3);  total (32/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

CAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGA  >  NC_000913/664945‑665242
                                                                                                                                                     |                                                                                                                                                    
cAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCTACTATAGGGCCGTTACATACCGCAgg                                                                                                                                                      <  1:938459/150‑1 (MQ=255)
 aGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGc                                                                                                                                                     >  2:986330/1‑150 (MQ=255)
 aGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGAATTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGc                                                                                                                                                     <  1:671393/150‑1 (MQ=255)
    cccGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGTCCGTTACATACCGCAGGCTGc                                                                                                                                                  <  1:1558357/150‑1 (MQ=255)
       accTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCACCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGt                                                                                                                                               <  2:1657528/149‑1 (MQ=255)
        ccTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTa                                                                                                                                              <  1:915446/150‑1 (MQ=255)
        ccTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTa                                                                                                                                              >  1:365757/1‑150 (MQ=255)
        ccTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTATTACTGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTa                                                                                                                                              >  2:1633761/1‑150 (MQ=255)
          tGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTTGCGCTGGTAATAATTATTTGCCGTCGAGTTCAGTGGTTCGAAACGCGTGTCCGCCCTGCTAATTTGAACTCTATGTCCGTTACATACCGCAGCCTGCGGTatg                                                                                                                                            <  2:1106077/150‑3 (MQ=255)
          tGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGCCTGCGGTatg                                                                                                                                            <  1:576252/150‑3 (MQ=255)
          ctacaaaaTCGTGCAGGATCCGTCTCGATTTTGTAGCTTTTACCGTCGCGCTGGTAATACTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAAATTCAACTATAGGGCCGTTACATACCGCCGGCTGCGGGaca                                                                                                                                            >  1:358484/6‑150 (MQ=255)
           ggCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGCCTGCGGTatgt                                                                                                                                           <  2:1338495/150‑4 (MQ=255)
            gCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACt                                                                                                                                          >  1:1357096/1‑150 (MQ=255)
             cTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACtt                                                                                                                                         >  1:27331/1‑150 (MQ=255)
                caacGAGACGGATCCTGCAGGATTTTGTAGCGCTGACCGTCGCGCGGGTAATCCTGAATTGCCGTCGCGCTCAGTGGTTCGAAACGCGGGTCCGCCCCGCAAATTTCAACTATAGGGCCGTTAAATACCGCAGGCTGCGGTACACTTAcc                                                                                                                                      <  1:269747/147‑1 (MQ=255)
                aaTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTAcc                                                                                                                                      <  1:938215/150‑1 (MQ=255)
                aaTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTCCCGTCGCGTTCAGTGGGTCGAACCGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTAcc                                                                                                                                      <  1:1043550/150‑1 (MQ=255)
                 aTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCg                                                                                                                                     <  1:426641/150‑1 (MQ=255)
                 aTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCg                                                                                                                                     >  2:1589457/1‑150 (MQ=255)
                 aTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCg                                                                                                                                     <  1:1185433/150‑1 (MQ=255)
                    gagaCGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCt                                                                                                                                  >  1:507934/1‑150 (MQ=255)
                    gagaCGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCt                                                                                                                                  <  2:1396458/150‑1 (MQ=255)
                     agaCGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCGGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTg                                                                                                                                 <  1:351623/150‑1 (MQ=255)
                           aTCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTg                                                                                                                           <  2:176732/150‑1 (MQ=255)
                               tGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACca                                                                                                                       <  1:138526/150‑1 (MQ=255)
                                gCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCGGCGGTACACTTACCGTCTGTTGCTGACcat                                                                                                                      >  2:1429126/1‑150 (MQ=255)
                                gCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGGTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACAAACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcat                                                                                                                      <  1:1159312/150‑1 (MQ=255)
                                 cACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcata                                                                                                                     >  1:1633615/1‑149 (MQ=255)
                                  aCGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcatca                                                                                                                    <  2:555712/150‑1 (MQ=255)
                                    gATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTAGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCatcatc                                                                                                                  <  2:316941/150‑1 (MQ=255)
                                     aTTTTGTAGCTTTTACGGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCg                                                                                                                 <  1:1029851/150‑1 (MQ=255)
                                         tGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTg                                                                                                             <  2:964631/150‑1 (MQ=255)
                                           tAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTa                                                                                                           <  1:1595484/150‑1 (MQ=255)
                                               ttttACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATg                                                                                                       <  1:1354355/150‑1 (MQ=255)
                                                   aCCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACAGTCTGTTGCTGACCATCATCGCTAGTACATGCcgg                                                                                                   >  2:705987/1‑149 (MQ=255)
                                                      gTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAg                                                                                                <  2:494023/150‑1 (MQ=255)
                                                       tCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGc                                                                                               <  2:844482/150‑1 (MQ=255)
                                                       tCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGc                                                                                               <  2:831565/150‑1 (MQ=255)
                                                          cgcTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCAtt                                                                                            <  1:986330/150‑1 (MQ=255)
                                                                   tCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgcg                                                                                   <  2:415296/150‑1 (MQ=255)
                                                                   tCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgcg                                                                                   <  1:1419571/150‑1 (MQ=255)
                                                                   accTGATTAGCCGTCGCGTTAAGTTGTTCGCAACGCTGGTAAGCCCCGCTAATTTCAACTTTAGGGCCGTTACATTCCGCAGGCTGCTGTACACTTAACGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCTAGAATTCCTGCcgcg                                                                                   <  2:107464/149‑1 (MQ=255)
                                                                        agttGCCGGCGCGTTCAGTGGTGCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCa                                                                              <  2:1561893/148‑1 (MQ=255)
                                                                        aTTTGCCGTCTCGTTCAGTGGTTCGAAAGGCTGGTGCGCCCCGCTAATTTCATCTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTTACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCa                                                                              <  1:269881/150‑1 (MQ=255)
                                                                             ccGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATcc                                                                         >  2:758993/1‑150 (MQ=255)
                                                                                 cgcgTTCAGGGCTTCGAAACGCGGGGCCGCCCCGCTAAGTTCACCTAGAGGGCCGTTACATACCGCCGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCAACGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAg                                                                     <  1:121016/150‑1 (MQ=255)
                                                                                   cgTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGcc                                                                   >  1:693352/1‑150 (MQ=255)
                                                                                      gcAGTGGGTCGAAACGCGGGGCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACt                                                                <  2:291406/149‑1 (MQ=255)
                                                                                             gtCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTCCATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACg                                                         <  2:277583/149‑1 (MQ=255)
                                                                                               cGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCa                                                       >  2:1291938/1‑150 (MQ=255)
                                                                                               cGAAACGCGGGTCCGCCCCGCTAATTTCAACTATACGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCa                                                       <  2:528792/150‑1 (MQ=255)
                                                                                                gAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCAt                                                      >  2:778517/1‑150 (MQ=255)
                                                                                                gAAACGATGGTCCGCCCCGCTAATTTCAACTATATGGCCGTTACAAACCGCAGGCTGAGTTACACTTTCCGTCTTTTTCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCAt                                                      <  2:1197101/150‑1 (MQ=255)
                                                                                                   aCGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTgc                                                   <  2:365757/150‑1 (MQ=255)
                                                                                                   aCGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTgc                                                   <  2:1559619/150‑1 (MQ=255)
                                                                                                   aCGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTgc                                                   <  2:507934/150‑1 (MQ=255)
                                                                                                      cgGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCa                                                <  1:155273/150‑1 (MQ=255)
                                                                                                       gggTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCAc                                               >  2:1348277/1‑150 (MQ=255)
                                                                                                       gggTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCAc                                               >  2:1499018/1‑150 (MQ=255)
                                                                                                       gggTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCAc                                               <  1:88070/150‑1 (MQ=255)
                                                                                                        ggTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCAcc                                              <  2:220101/150‑1 (MQ=255)
                                                                                                             gCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCGTTGCGCACCTCTTa                                         >  1:1050100/1‑150 (MQ=255)
                                                                                                               cccGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTaca                                       >  1:1466232/1‑150 (MQ=255)
                                                                                                                ccGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTacac                                      <  2:600405/150‑1 (MQ=255)
                                                                                                                       tttCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTc                               <  1:692335/150‑1 (MQ=255)
                                                                                                                        ttCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCg                              >  1:764490/1‑150 (MQ=255)
                                                                                                                        ttCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCg                              <  2:1590305/150‑1 (MQ=255)
                                                                                                                          cAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTTACCATCATAGCTTGTACATGCCGCGATCATTCCTGCCGCGATGCAGATCCAGAGCCACTTCTTACGCATTGCGCACCTCTTCCACGATTTTCGac                            >  2:129560/1‑150 (MQ=255)
                                                                                                                            aCTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGacaa                          >  1:1054619/1‑150 (MQ=255)
                                                                                                                             cTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGCGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGacaac                         >  1:1560667/1‑150 (MQ=255)
                                                                                                                              tataGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTACTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGacaaca                        >  2:1554063/1‑150 (MQ=255)
                                                                                                                               ataGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACAt                       >  1:921158/1‑150 (MQ=255)
                                                                                                                                 aGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCTCGATCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCATTTCGACAACAttt                     >  2:689739/1‑150 (MQ=255)
                                                                                                                                   ggCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTc                   <  2:600964/150‑1 (MQ=255)
                                                                                                                                    gCCGCTACATACCGCAGGCTGCCGTACACTTACCGTCAGTTGCTGACCATCCTCGCGTGTACAAGCCGCGAGCATTCCGGCCGCGATGCAGATCCCGAGCCACTGCATACTCATTGCCCACCTCTAACACGCTTTTCTACACCATTTTcc                  <  2:655308/150‑1 (MQ=255)
                                                                                                                                     ccGTTACAGACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCt                 <  1:866269/150‑1 (MQ=255)
                                                                                                                                      gtaagtgtACCGCAGCCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtg                >  2:749383/8‑150 (MQ=255)
                                                                                                                                      cGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtg                >  2:1601713/1‑150 (MQ=255)
                                                                                                                                       gTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtgt               <  2:693352/150‑1 (MQ=255)
                                                                                                                                        ttACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtgtg              <  2:1338063/150‑1 (MQ=255)
                                                                                                                                         agtgtACCGCAGCCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGCTGCAGATCCCGAGCCACTGCTTCCGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTgg             >  2:1505422/5‑150 (MQ=255)
                                                                                                                                          aCATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTggg            >  1:599943/1‑150 (MQ=255)
                                                                                                                                          aCATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTggg            <  2:559792/150‑1 (MQ=255)
                                                                                                                                              aCCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGtgtg        <  1:481338/150‑1 (MQ=255)
                                                                                                                                               ccGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTgg       >  2:1160863/1‑150 (MQ=255)
                                                                                                                                                 gCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGAt     >  1:419291/1‑150 (MQ=255)
                                                                                                                                                    ggCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGa  <  1:1499018/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGA  >  NC_000913/664945‑665242

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: