Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,785,780 T→C 6.0% K445E (AAA→GAA)  ppsA ← phosphoenolpyruvate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,785,7800TC 6.0% 203.1 / 2.8 66K445E (AAA→GAA) ppsAphosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base T (26/36);  new base C (2/2);  total (28/38)
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.92e-01

CGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACGATGGCAGATGCTTTCTTCATGATCGGTTCCCAGTCCGGGTCGGTCA  >  NC_000913/1785634‑1785926
                                                                                                                                                  |                                                                                                                                                  
cGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCa                                                                                                                                                  <  7:50775/149‑1 (MQ=255)
        cGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCt                                                                                                                                          <  1:170878/149‑1 (MQ=255)
        cGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCt                                                                                                                                          <  7:99608/149‑1 (MQ=255)
         gACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGTAGCTTTTCACGCTAAATTCCAGGATGTCCTCATTGAGGTATCGGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCTTCCGTTCTg                                                                                                                                         <  6:369766/149‑1 (MQ=255)
           cGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGtt                                                                                                                                       <  4:308622/149‑1 (MQ=255)
           cGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGtt                                                                                                                                       <  1:42877/149‑1 (MQ=255)
             ttCATCATCACTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTg                                                                                                                                      >  7:36315/1‑148 (MQ=255)
               catcatcaCTTTCACCGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACTGTCTTTCATCCGTTCTGTTGCAt                                                                                                                                   >  4:156026/1‑149 (MQ=255)
                atcatcaCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCAtc                                                                                                                                  >  1:212058/1‑149 (MQ=255)
                 tcatcaCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCAtct                                                                                                                                 >  6:315525/1‑149 (MQ=255)
                  catcaCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCAtctc                                                                                                                                >  7:187040/1‑149 (MQ=255)
                    tcaCTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCca                                                                                                                              >  1:381925/1‑149 (MQ=255)
                       cTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAg                                                                                                                           <  5:187843/149‑1 (MQ=255)
                       cTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAg                                                                                                                           >  7:228075/1‑149 (MQ=255)
                          tCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGc                                                                                                                                     <  5:86493/136‑1 (MQ=255)
                          tCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGc                                                                                                                                     >  6:86493/1‑136 (MQ=255)
                           cAACGGCAGATCCGGCATCGTTTCTACTCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGccc                                                                                                                        <  2:235394/148‑1 (MQ=255)
                            aaCGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCa                                                                                                                       <  8:228075/148‑1 (MQ=255)
                            aaCGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCAc                                                                                                                      >  2:72542/1‑149 (MQ=255)
                            aaCGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTCCATCCGTTCTGTTGCATCTCCACATCCCAc                                                                                                                      >  4:216135/1‑149 (MQ=255)
                            aaCGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTCCATCCGTTCTGTTGCATCTCCACAGCCCa                                                                                                                       >  4:211276/1‑148 (MQ=255)
                             aCGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCAc                                                                                                                      >  3:171478/1‑148 (MQ=255)
                              cGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTa                                                                                                                    <  5:212883/149‑1 (MQ=255)
                                   gATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTAccgcc                                                                                                               >  7:86778/1‑149 (MQ=255)
                                     tCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCgg                                                                                                             <  5:220001/149‑1 (MQ=255)
                                           aTCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGcc                                                                                                       >  5:239630/1‑149 (MQ=255)
                                           aTCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGcc                                                                                                       <  1:332056/149‑1 (MQ=255)
                                              gTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCAtct                                                                                                                                 <  7:369772/120‑1 (MQ=255)
                                              gTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCAtct                                                                                                                                 >  8:369772/1‑120 (MQ=255)
                                                 tCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGtt                                                                                                  <  4:171478/148‑1 (MQ=255)
                                                     cGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCAc                                                                                               <  8:36315/147‑1 (MQ=255)
                                                                   aCGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATcgccgc                                                                               <  8:106367/149‑1 (MQ=255)
                                                                     gCTAAATTCCAGCAACTCCGCATAGACGTACCCGGTAGCACCTTCGGCACAAGAAACAGTGACGTGCTCACCGTCTTCCATCCGTTCTGTGGCATCTCCACAGCCCACTACCGCCGGAATGCACAGTTCACGAGCGATGATCGCCGCGt                                                                             <  3:216135/149‑1 (MQ=255)
                                                                     gCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCCCTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGt                                                                             >  1:256252/1‑149 (MQ=255)
                                                                     gCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGt                                                                             <  2:212058/149‑1 (MQ=255)
                                                                     gCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTCCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGt                                                                             <  3:211276/149‑1 (MQ=255)
                                                                      cTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGt                                                                             >  3:126699/1‑148 (MQ=255)
                                                                         aaTTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGAc                                                                          <  3:175259/148‑1 (MQ=255)
                                                                           ttCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAgg                                                                       >  4:194717/1‑149 (MQ=255)
                                                                              cAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTAc                                                                    <  4:81548/149‑1 (MQ=255)
                                                                              cAGCAACTCCGCAGAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGGTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTAc                                                                    <  4:341169/149‑1 (MQ=255)
                                                                               aGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACg                                                                   <  8:187040/149‑1 (MQ=255)
                                                                                      ccGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATcg                                                                                   >  4:219398/1‑126 (MQ=255)
                                                                                      ccGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATcg                                                                                   <  3:219398/126‑1 (MQ=255)
                                                                                       cGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGcc                                                             <  6:239630/147‑1 (MQ=255)
                                                                                        gCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACg                                                          <  8:86778/149‑1 (MQ=255)
                                                                                                gTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGa                                                  <  5:315525/149‑1 (MQ=255)
                                                                                                gTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGa                                                  <  3:288896/149‑1 (MQ=255)
                                                                                                  aaCCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGAc                                                 >  6:123273/1‑148 (MQ=255)
                                                                                                       gTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACGATgg                                            >  8:210214/1‑148 (MQ=255)
                                                                                                       gTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACGATGGc                                           <  1:72542/149‑1 (MQ=255)
                                                                                                                     acaAGAAACAGTGACGTTCTCACCGTCTTTCTTCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACGATGGCAGATGCTTTCTTCa                             <  4:184849/149‑1 (MQ=255)
                                                                                                                      caAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCgatgat                                                                                     <  4:365368/92‑1 (MQ=255)
                                                                                                                      caAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCgatgat                                                                                     >  3:365368/1‑92 (MQ=255)
                                                                                                                          aaaCAGTGACGTTCTCACTGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACGATGGCAGATGCTTTCTTCATGATc                        <  3:156026/149‑1 (MQ=255)
                                                                                                                          aaaCAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACGATGGCAGATGCTTTCTTCATGATc                        >  8:314618/1‑149 (MQ=255)
                                                                                                                            aCAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACGATGGCAGATGCTTTCTTCATGATCgg                      >  5:256355/1‑149 (MQ=255)
                                                                                                                                    gTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGa                                                  >  3:161255/1‑113 (MQ=255)
                                                                                                                                    gTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGa                                                  <  4:161255/113‑1 (MQ=255)
                                                                                                                                    gTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACGATGGCAGATGCTTTCTTCATGATCGGTTCCCAGt              <  6:256355/149‑1 (MQ=255)
                                                                                                                                        tcACCGTCTTTCATCCGTTCTGTTGCATCTCCACCGCCCACTACAGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACCATGGCAGATGCATTCTTCATGATCGGTTCACAgaacgg          >  5:255200/1‑144 (MQ=255)
                                                                                                                                          aCCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCCCTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACGATGGCAGATGCTTTCTTCATGATCGGTTCCCAGTCCgg          <  2:256252/147‑1 (MQ=255)
                                                                                                                                          aCCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACGATGGCAGATGCTTTCTTCATGATCGGTTCCCAGTCCgg          <  4:126699/147‑1 (MQ=255)
                                                                                                                                          aCCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACGATGGCAGATGCTTTCTTCATGATCGGTTCCCAGTCCgg          <  1:160229/147‑1 (MQ=255)
                                                                                                                                          aCCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACGATGGCAGATGCTTTCTTCATGATCGGTTCCCAGTCCGggt        >  7:170764/1‑149 (MQ=255)
                                                                                                                                            cGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACGATGGCAGATGCTTTCTTCATGATCGGTTCCCAGTCCGggtcg      >  5:162642/1‑149 (MQ=255)
                                                                                                                                                tttCATCCGTTCTGTTGCATCTCCACCGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACGATGGCAGATGCTTTCTTCATGATCGGTTCCCAGTCCGGGTCGGTCa  >  3:306207/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGGGTTACCGACGTTCATCATCACTTTCAACGGCAGATCCGGCATCGTTTCTACGCTGGAGCTTTTCACGCTAAATTCCAGCAACTCCGCATAGACGTAACCGGTATCACCTTCGGCACAAGAAACAGTGACGTTCTCACCGTCTTTCATCCGTTCTGTTGCATCTCCACAGCCCACTACCGCCGGAATGCCCAGTTCACGAGCGATGATCGCCGCGTGACAGGTACGACCGCCACGGTTGGTGACGATGGCAGATGCTTTCTTCATGATCGGTTCCCAGTCCGGGTCGGTCA  >  NC_000913/1785634‑1785926

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