Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,116,857 G→T 6.6% M17I (ATG→ATT)  yceA → UPF0176 protein YceA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,116,8570GT 6.6% 234.5 / 2.2 70M17I (ATG→ATT) yceAUPF0176 protein YceA
Reads supporting (aligned to +/- strand):  ref base G (29/36);  new base T (3/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACG  >  NC_000913/1116710‑1117005
                                                                                                                                                   |                                                                                                                                                    
gcgcAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCACTTCCAACGACGCGCTAAAAGCCAAAATGtt                                                                                                                                                    >  2:942836/1‑150 (MQ=255)
 cgcAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTg                                                                                                                                                   <  1:1211419/150‑1 (MQ=255)
   cAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGc                                                                                                                                                 <  1:1045451/150‑1 (MQ=255)
        gTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTgaga                                                                                                                                            >  1:447680/1‑150 (MQ=255)
         tAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAACATTTTGGCttttag                                                                                                                                           <  1:506108/150‑6 (MQ=255)
          aGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgc                                                                                                                                          >  2:89859/1‑150 (MQ=255)
          aGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgc                                                                                                                                          >  2:1215138/1‑150 (MQ=255)
           gCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTCTCAGCCAACATTTTGGCttttagcg                                                                                                                                         <  1:1262901/150‑8 (MQ=255)
            cGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCGAAAATGTTGGCTGAGAgcga                                                                                                                                        <  2:1163097/150‑1 (MQ=255)
            cGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgcga                                                                                                                                        >  1:955709/1‑150 (MQ=255)
            cGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATATCGCAGGAACCGTACAACATGACAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgcga                                                                                                                                        >  1:650154/1‑150 (MQ=255)
                     ggtgGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCTCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgcgaacagctaac                                                                                                                               <  2:981439/149‑1 (MQ=255)
                        ggATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAaccac                                                                                                                            >  2:1503660/1‑150 (MQ=255)
                              tAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTc                                                                                                                      <  2:1539429/150‑1 (MQ=255)
                                     cGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCTAACCACCATTTCGTTTTAc                                                                                                               >  1:908430/1‑150 (MQ=255)
                                      gATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgcgaaccgcga                                                                                                                                 >  1:586985/1‑131 (MQ=255)
                                      gATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgcgaaccgcga                                                                                                                                 <  2:586985/131‑1 (MQ=255)
                                      gATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACa                                                                                                              >  1:307623/1‑150 (MQ=255)
                                      gATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACa                                                                                                              <  2:955709/150‑1 (MQ=255)
                                      gATAAATGCTATATCATGCAGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACa                                                                                                              >  2:1547545/1‑150 (MQ=255)
                                        tAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAg                                                                                                            <  1:412330/150‑1 (MQ=255)
                                        tAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAg                                                                                                            <  1:406246/150‑1 (MQ=255)
                                               tataTCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTcc                                                                                                     <  1:1178654/150‑1 (MQ=255)
                                                atatCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCGAGTGTTACACAACCGCATATCCAACTACTGTCTAAAAGTCAAAATGTTGTCTTAGAGGGTACCGGCAACCAGCATTTGTTTTTACAAGTATTTAca                                                                                                    <  2:684880/150‑1 (MQ=255)
                                                atatCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCca                                                                                                    >  2:1662436/1‑150 (MQ=255)
                                                atatCATGCCGCGGAGTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAATCCAAAATGTTGGCTGAGAGCGAACCGCGAACCTCCATTTCGTTTTACAAGTATTTCca                                                                                                    <  2:919708/150‑1 (MQ=255)
                                                  atCATGCCGCGAATTTTTACTTTCCCATCTCGCAGGACCCGTACACCATGCCAGTGTTACACAACCGCATTTACAACGCCGCGCTAAAAGCCAAAATGTTGGCTGAGAGCTAACCGCGAACCACCAATTCGTTTTACAAGTATTTCcaca                                                                                                  >  2:149363/1‑150 (MQ=255)
                                                      tGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCCCCATTTCGTTTTACAAGTATTTCCCCATCGc                                                                                              >  2:1278606/1‑150 (MQ=255)
                                                       gCCGCTGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGcc                                                                                             <  1:874651/150‑1 (MQ=255)
                                                          gcgGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGAt                                                                                          <  1:75494/150‑1 (MQ=255)
                                                            ggATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATcc                                                                                        <  2:447680/150‑1 (MQ=255)
                                                              aTTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTa                                                                                      <  1:1503660/150‑1 (MQ=255)
                                                               tttttACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATTATCGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTAAAAGAATTTCCACATCGCCGATCCaaa                                                                                     >  2:330878/1‑150 (MQ=255)
                                                                        tCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCg                                                                            >  2:702413/1‑150 (MQ=255)
                                                                        tCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCg                                                                            <  1:801664/150‑1 (MQ=255)
                                                                          ccATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTg                                                                          >  1:125821/1‑150 (MQ=255)
                                                                            aTCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGAc                                                                        <  1:1215138/150‑1 (MQ=255)
                                                                              ctcGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGc                                                                      <  2:650154/150‑1 (MQ=255)
                                                                               tcGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAACTATTTCCAAATAGCCGATCCTAAGGCGACCCGTGACGCt                                                                     >  1:500663/1‑150 (MQ=255)
                                                                                    ggAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGAGGCTCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATATCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACGCGTGACGCTTtata                                                                <  2:745116/150‑1 (MQ=255)
                                                                                      aaCCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATc                                                              >  2:1345408/1‑150 (MQ=255)
                                                                                      aaCCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATc                                                              <  1:130802/150‑1 (MQ=255)
                                                                                       aCCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCa                                                             <  2:154884/150‑1 (MQ=255)
                                                                                        ccGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAg                                                            >  2:747222/1‑150 (MQ=255)
                                                                                          gTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCt                                                          >  2:713/1‑150 (MQ=255)
                                                                                           tACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACGATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATTTCAGCTg                                                         <  2:318422/150‑1 (MQ=255)
                                                                                               ccATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTa                                                     >  1:335780/1‑150 (MQ=255)
                                                                                               ccAAGCCAGTGTTACACAACCGCATTTCCAACGACTCGCTAAAATCAAAAATTTTGGCTGAGAGCGAACCGCGACCCACAATTTCGTTTAACAAGTATTTACACATCGCCGATCCTAAGGCTACACGTGACGCTTTATAAAACCTGTTTa                                                     >  2:1259107/1‑150 (MQ=255)
                                                                                                  tGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACcg                                                  <  1:420843/150‑1 (MQ=255)
                                                                                                  tGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGCTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACcg                                                  <  2:1150622/150‑1 (MQ=255)
                                                                                                   gCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCTAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACcgc                                                 <  1:455331/150‑1 (MQ=255)
                                                                                                     cAGTGTTACACAACCGCATTTCCAAAGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACTAGTTTTTCCACATCGCAGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCgcgc                                               <  1:1281320/150‑1 (MQ=255)
                                                                                                          ttACACAACCGCACTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAAt                                          <  1:942836/150‑1 (MQ=255)
                                                                                                             cacaACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGtt                                       >  2:1516333/1‑150 (MQ=255)
                                                                                                               caACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGtttt                                     >  2:44807/1‑150 (MQ=255)
                                                                                                                        ttCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAgt                            <  1:702413/150‑1 (MQ=255)
                                                                                                                         tCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAgtg                           <  1:402489/150‑1 (MQ=255)
                                                                                                                         tCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCGTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAgtg                           <  2:899090/150‑1 (MQ=255)
                                                                                                                           cAACGACGCGCTAAAAGCCAAACTGGTGGCGGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTCAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTa                         <  2:1364346/150‑1 (MQ=255)
                                                                                                                             acgacgCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATc                       >  1:764623/1‑150 (MQ=255)
                                                                                                                                acgCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTgg                    <  2:307317/150‑1 (MQ=255)
                                                                                                                                  gcgcTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGgcg                  <  1:1313854/150‑1 (MQ=255)
                                                                                                                                   cgcTAAAAGCCAAAATGTTCGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCAACATCGCCGAACCTAAGGCTACCCGACACGCTTTATATCAGATGTTTACAGAGCTGAATGTTTTTGTGCGAGTGTATATGgagc                 >  2:824561/1‑150 (MQ=255)
                                                                                                                                    gcTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCa                >  1:841203/1‑150 (MQ=255)
                                                                                                                                        aaaGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTACCACATCGCCGATCATAAGGCGACCCGTGACGCTTTATATCAGATGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATATGGCGCATGAg            >  1:641565/1‑150 (MQ=255)
                                                                                                                                        aaaGCCAAAAAGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAg            <  2:653501/150‑1 (MQ=255)
                                                                                                                                         aaGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAgg           >  1:586121/1‑150 (MQ=255)
                                                                                                                                             cAACATTTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCAt       >  1:74889/1‑150 (MQ=255)
                                                                                                                                                aaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGAGCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTaa    <  1:1547545/150‑1 (MQ=255)
                                                                                                                                                  tGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAAcg  >  2:417082/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGCAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACG  >  NC_000913/1116710‑1117005

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: ‑