Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 665,094 G→C 8.6% 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 8.6% 272.1 / 13.0 81P33A (CCT→GCT) rlpAseptal ring protein, suppressor of prc, minor lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (39/35);  new base C (3/4);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

TCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGA  >  NC_000913/664948‑665238
                                                                                                                                                  |                                                                                                                                                
tCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCt                                                                                                                                                >  4:466068/1‑149 (MQ=255)
 cccGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCt                                                                                                                                                >  2:196716/1‑148 (MQ=255)
  ccGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGCCTGc                                                                                                                                              >  4:791656/1‑149 (MQ=255)
       tGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGCCTGCGGTat                                                                                                                                         <  3:791656/149‑2 (MQ=255)
       tGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGCCTGCGGTat                                                                                                                                         <  3:327515/149‑2 (MQ=255)
         gCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTacac                                                                                                                                       <  4:172459/149‑1 (MQ=255)
         gCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCATTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTacac                                                                                                                                       >  4:741582/1‑149 (MQ=255)
          cTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTacac                                                                                                                                       >  2:352216/1‑148 (MQ=255)
            aaaTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACtt                                                                                                                                     >  2:822422/1‑148 (MQ=255)
             aaTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTa                                                                                                                                    >  1:158826/1‑148 (MQ=255)
             aaTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTAc                                                                                                                                   >  4:44256/1‑149 (MQ=255)
             aaTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGCCTGCGGTatgtaacg                                                                                                                                   <  4:724057/149‑8 (MQ=255)
              aTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTAcc                                                                                                                                  >  4:257366/1‑149 (MQ=255)
              aTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTAc                                                                                                                                   <  2:158826/148‑1 (MQ=255)
                cGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCg                                                                                                                                 >  1:481140/1‑148 (MQ=255)
                cGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCg                                                                                                                                 >  1:300017/1‑148 (MQ=255)
                   gaCGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGCCTGCGGTa                                                                                                                                          <  1:244315/136‑1 (MQ=255)
                   gaCGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGCCTGCGGTa                                                                                                                                          >  2:244315/1‑136 (MQ=255)
                        aTCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCt                                                                                                                        >  3:793662/1‑149 (MQ=255)
                          ccTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTg                                                                                                                       <  2:300017/148‑1 (MQ=255)
                          ccTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGa                                                                                                                      <  4:296036/149‑1 (MQ=255)
                               aCGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcat                                                                                                                  >  2:339331/1‑148 (MQ=255)
                                  aTTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcatcat                                                                                                               >  1:37117/1‑148 (MQ=255)
                                   ttttGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCg                                                                                                             >  4:890911/1‑149 (MQ=255)
                                    tttGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGc                                                                                                            <  3:44256/149‑1 (MQ=255)
                                    tttGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGc                                                                                                            >  4:163334/1‑149 (MQ=255)
                                      tGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCtt                                                                                                          <  4:793662/149‑1 (MQ=255)
                                        tAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGt                                                                                                        <  3:257366/149‑1 (MQ=255)
                                         aGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGt                                                                                                        >  1:400160/1‑148 (MQ=255)
                                         aGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGt                                                                                                        <  2:481140/148‑1 (MQ=255)
                                             tttACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACAt                                                                                                    >  2:108620/1‑148 (MQ=255)
                                                 ccGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCcgc                                                                                               <  3:466068/149‑1 (MQ=255)
                                                   gTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGa                                                                                             >  3:477315/1‑149 (MQ=255)
                                                       agcTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCa                                                                                          <  1:719485/147‑1 (MQ=255)
                                                             tAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTg                                                                                    >  1:395255/1‑148 (MQ=255)
                                                                  cTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgcg                                                                               <  1:339331/148‑1 (MQ=255)
                                                                  cTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgcg                                                                               <  2:422836/148‑1 (MQ=255)
                                                                    gATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGAt                                                                             >  2:530330/1‑148 (MQ=255)
                                                                      tttGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGc                                                                           <  2:400160/148‑1 (MQ=255)
                                                                        tGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAg                                                                         <  1:315014/148‑1 (MQ=255)
                                                                         gCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGa                                                                        <  2:822450/148‑1 (MQ=255)
                                                                           cGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATcc                                                                     <  3:163334/149‑1 (MQ=255)
                                                                            gTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATccc                                                                    <  3:488524/149‑1 (MQ=255)
                                                                               gcgTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAg                                                                 >  3:358428/1‑149 (MQ=255)
                                                                               gcgTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACAGACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCTTCGCTTGAACATGCCGCGAGCATTCCTGCCGCGATGCAGATCTCGAg                                                                 >  4:563263/1‑149 (MQ=255)
                                                                                   tCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCAc                                                             <  4:656673/149‑1 (MQ=255)
                                                                                       tGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGc                                                          >  1:812268/1‑148 (MQ=255)
                                                                                         gTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTa                                                       >  4:695627/1‑149 (MQ=255)
                                                                                           tCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTAc                                                      >  2:457097/1‑148 (MQ=255)
                                                                                              aaaCGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTCCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCa                                                   >  2:265278/1‑148 (MQ=255)
                                                                                                aCGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTg                                                <  4:837725/149‑1 (MQ=255)
                                                                                                aCGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTg                                                >  3:720055/1‑149 (MQ=255)
                                                                                                     ggTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCa                                            >  1:638875/1‑148 (MQ=255)
                                                                                                        ccGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACct                                         <  1:352216/148‑1 (MQ=255)
                                                                                                         cGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACctc                                        <  2:841140/148‑1 (MQ=255)
                                                                                                            cccGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTcc                                                                                      >  4:566802/1‑99 (MQ=255)
                                                                                                            cccGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTcc                                                                                      <  3:566802/99‑1 (MQ=255)
                                                                                                            cccGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTa                                     >  2:319429/1‑148 (MQ=255)
                                                                                                              cGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGc                                                          >  1:127495/1‑125 (MQ=255)
                                                                                                              cGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGc                                                          <  2:127495/125‑1 (MQ=255)
                                                                                                                     ttCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCtttt                            <  1:129690/148‑1 (MQ=255)
                                                                                                                       cAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCg                          >  1:792237/1‑148 (MQ=255)
                                                                                                                         aCTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGac                        >  1:105814/1‑148 (MQ=255)
                                                                                                                         aCTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGac                        >  2:401099/1‑148 (MQ=255)
                                                                                                                           tataGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGacaac                     >  4:164740/1‑149 (MQ=255)
                                                                                                                              aGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACAtt                  >  4:716057/1‑149 (MQ=255)
                                                                                                                                ggCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACAttt                 <  2:303533/148‑1 (MQ=255)
                                                                                                                                 gCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACAtttt                <  2:812268/148‑1 (MQ=255)
                                                                                                                                 gCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTc               <  4:358428/149‑1 (MQ=255)
                                                                                                                                     ttACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtgt           <  3:695627/149‑1 (MQ=255)
                                                                                                                                     ttACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtg            <  1:401099/148‑1 (MQ=255)
                                                                                                                                       aCATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtgtg          <  1:457097/148‑1 (MQ=255)
                                                                                                                                       aCATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTgg         >  4:755548/1‑149 (MQ=255)
                                                                                                                                       aCATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGACGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTgg         <  3:741582/149‑1 (MQ=255)
                                                                                                                                           aCCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGgtgt     >  3:161352/1‑149 (MQ=255)
                                                                                                                                           aCCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGgtg      <  2:380545/148‑1 (MQ=255)
                                                                                                                                           aCCGCAGCCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGgtg      >  2:869399/1‑148 (MQ=255)
                                                                                                                                              gCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTgg   >  2:662210/1‑148 (MQ=255)
                                                                                                                                              gCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTgg   <  1:319429/148‑1 (MQ=255)
                                                                                                                                              gCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTgg   <  1:270973/148‑1 (MQ=255)
                                                                                                                                               cAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCTTTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGa  <  2:269845/148‑1 (MQ=255)
                                                                                                                                                  gCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtgt           >  4:209644/1‑136 (MQ=255)
                                                                                                                                                  gCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtgt           <  3:209644/136‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGA  >  NC_000913/664948‑665238

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