Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,085,098 C→T 5.4% P36L (CCT→CTT)  phoH → ATP‑binding protein; putative PhoH family P‑loop ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,085,0980CT 5.4% 234.6 / 2.5 74P36L (CCT→CTT) phoHATP‑binding protein; putative PhoH family P‑loop ATPase
Reads supporting (aligned to +/- strand):  ref base C (34/36);  new base T (2/2);  total (36/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 = 1.00e+00

GGCAGTGATAGTCAGTCACTATGGAGATCGCGGATGGTAACGTCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGCCAACGGCGCGTGTGATGAATA  >  NC_000913/1084959‑1085245
                                                                                                                                           |                                                                                                                                                   
ggCAGTGATAGTCAGTCACTATGGAGATCGCGGATGGTAACGTCGTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTc                                                                                                                                            <  7:138363/149‑1 (MQ=255)
  cAGTGATAGGCAGTCACTATGGAGATCGCGGATGGTAACGTCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCAt                                                                                                                                          <  2:29864/149‑1 (MQ=255)
    gTGATAGTCAGTCACTATGGAGATCGCGGATGGTAACGTCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCAtgt                                                                                                                                        <  1:285423/149‑1 (MQ=255)
        tAGTCAGTCACTATGGAGATCGCGGATGGTAACGTCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTa                                                                                                                                    <  4:247004/149‑1 (MQ=255)
          gtcaggcaCTATGGAGATCGCGGATGGTAACGTCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTagc                                                                                                                                  <  5:49351/149‑3 (MQ=255)
              gtcaCTATGGAGATCGCGGATGGTAACGTCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACg                                                                                                                              <  3:224686/149‑1 (MQ=255)
                 aCTATGGAGATCGCGGATGGTAACGTCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCg                                                                                                                           >  8:240516/1‑149 (MQ=255)
                  cTATGGAGATCGCGGATGGTAACGTCCTGTATTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGt                                                                                                                          <  4:160529/149‑1 (MQ=255)
                  cTATGGAGATCGCGGATGGTAACGTCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCg                                                                                                                           >  1:281068/1‑148 (MQ=255)
                     tGGAGATCGCGGATGGTAACGTCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCAtgt                                                                                                                                        <  2:21823/132‑1 (MQ=255)
                     tGGAGATCGCGGATGGTAACGTCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCAtgt                                                                                                                                        >  1:21823/1‑132 (MQ=255)
                     tGGAGATCGCGGATGGTAACGTCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTaaa                                                                                                                       >  2:188904/1‑149 (MQ=255)
                     tGGAGATCGCGGATGGTAACGTCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTaaa                                                                                                                       <  7:240516/149‑1 (MQ=255)
                          aTCGCGGATGGTAACGTCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTg                                                                                                                  >  2:51582/1‑149 (MQ=255)
                                       aCGTCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGttt                                                                                                     <  7:329721/149‑1 (MQ=255)
                                          tCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCCTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTc                                                                                                  >  6:322867/1‑149 (MQ=255)
                                            cTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCgg                                                                                                <  3:121893/149‑1 (MQ=255)
                                               tACTGGACTTGTTTTAGATAATCAACGCGCCACAACTCGCGGAGTTTTCTCCTCCGGTAGGCATCTCGTAACTTGGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGACGTAAACGGTGCGCGTAGCCGTTTTTCGGctg                                                                                             >  6:94673/1‑147 (MQ=255)
                                                 cTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCAttt                                                                                           >  6:287424/1‑149 (MQ=255)
                                                 cTGGACATGTTTTAGATAATCAACGCGCAACAACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCAtt                                                                                            <  5:147978/148‑1 (MQ=255)
                                                    gACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTAt                                                                                        <  1:51582/149‑1 (MQ=255)
                                                     aCATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTATTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATc                                                                                       <  4:275037/149‑1 (MQ=255)
                                                       aTGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGTTGCGCGTAGCCGTTTTTCGGCATTTTATCCa                                                                                     <  1:43093/149‑1 (MQ=255)
                                                          ttttAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCaaaa                                                                                  <  3:124629/149‑1 (MQ=255)
                                                          ttttAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCaa                                                                                    >  3:286253/1‑147 (MQ=255)
                                                            ttAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAg                                                                                 >  7:106990/1‑148 (MQ=255)
                                                               gATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCTTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCaaa                                                                             >  1:230017/1‑149 (MQ=255)
                                                               gATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCTTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCaaa                                                                             >  1:130810/1‑149 (MQ=255)
                                                                aTAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAAc                                                                            >  5:110752/1‑149 (MQ=255)
                                                                        cgcgCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCTTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCtt                                                                    <  2:130810/149‑1 (MQ=255)
                                                                        cgcgCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCTTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCtt                                                                    <  2:230017/149‑1 (MQ=255)
                                                                             caccacTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAg                                                               <  4:84444/149‑1 (MQ=255)
                                                                              accacTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAgg                                                              <  1:188904/149‑1 (MQ=255)
                                                                              accacTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAg                                                               >  4:146587/1‑148 (MQ=255)
                                                                                cacTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGaa                                                            <  4:38002/149‑1 (MQ=255)
                                                                                cacTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGaa                                                            <  5:122182/149‑1 (MQ=255)
                                                                                 acTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAg                                                           >  6:261100/1‑149 (MQ=255)
                                                                                  cTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCCTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGc                                                          <  5:322867/149‑1 (MQ=255)
                                                                                     gcgGAGTTTTCTCCTCCGGTAGTCATCTCTTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTTCAAGGAAGCCa                                                        <  2:281068/148‑1 (MQ=255)
                                                                                      cgGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAAc                                                      <  5:287424/149‑1 (MQ=255)
                                                                                       ggAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAAcc                                                     >  5:325167/1‑149 (MQ=255)
                                                                                        gAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAAccc                                                    <  6:43388/149‑1 (MQ=255)
                                                                                           ttttCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCctct                                                 >  8:10302/1‑149 (MQ=255)
                                                                                           ttttCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCctct                                                 >  8:5584/1‑149 (MQ=255)
                                                                                                  ctccGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTc                                          >  5:281484/1‑149 (MQ=255)
                                                                                                   tccGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCg                                         >  5:188296/1‑149 (MQ=255)
                                                                                                      ggTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTgc                                      <  4:286253/149‑1 (MQ=255)
                                                                                                       gTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTgcg                                     >  1:296696/1‑149 (MQ=255)
                                                                                                       gTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTgcg                                     >  6:159028/1‑149 (MQ=255)
                                                                                                       gTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTgcg                                     >  5:268736/1‑149 (MQ=255)
                                                                                                        tAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCTTTTTTCGGCATTTTATCCAAAAGCAATCTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTgcg                                     >  2:65806/1‑148 (MQ=255)
                                                                                                             aTCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCAttt                                                                                           >  7:143784/1‑89 (MQ=255)
                                                                                                             aTCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCAttt                                                                                           <  8:143784/89‑1 (MQ=255)
                                                                                                             aTCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAAt                               >  6:118259/1‑149 (MQ=255)
                                                                                                                     aCTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTg                        <  7:5584/148‑1 (MQ=255)
                                                                                                                     aCTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTg                        <  7:10302/148‑1 (MQ=255)
                                                                                                                     aCTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGc                       <  6:325167/149‑1 (MQ=255)
                                                                                                                      cTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGcc                      <  5:118259/149‑1 (MQ=255)
                                                                                                                       tttGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGcc                      >  1:126388/1‑148 (MQ=255)
                                                                                                                       tttGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGcc                      >  3:172683/1‑148 (MQ=255)
                                                                                                                         tGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTAACGACGCCGTAAACGGTGCGCGTAGCCGTTTTTGGGCTTTTTTTCCAAAAGCAAACTATGGCTTACAAGGCAGGCAACACTCTGATGTTCGTGAGCTAAATCGCGCTGccga                    <  4:32989/148‑3 (MQ=255)
                                                                                                                          gCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGCCAACg                  >  6:204575/1‑149 (MQ=255)
                                                                                                                           cATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCgcgc                          <  7:97751/140‑1 (MQ=255)
                                                                                                                           cATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCgcgc                          >  8:97751/1‑140 (MQ=255)
                                                                                                                                  aaCCTCATCCTTTCTTTTGATGTGTTACAGACGCCGTAAAAGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCTAACTATGGCTTACAAGGATGGCAACCCTCTGATGTTCGTGCGCAAAATCGCGCTGCCAACGGCGCgtgt          <  5:319706/149‑1 (MQ=255)
                                                                                                                                  aaCCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGACTCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGCCAACGCCGCgtgt          >  3:29782/1‑149 (MQ=255)
                                                                                                                                  aaCCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGCCAACGGCGCgtgt          >  3:299581/1‑149 (MQ=255)
                                                                                                                                   aCCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGCCAACGGCGCGtgtg         <  6:268736/149‑1 (MQ=255)
                                                                                                                                     cTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTGCGTGCACATAATCGCGCTGCCAACGGCGCGTGtgat       <  1:303152/149‑1 (MQ=255)
                                                                                                                                        aTCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGCCAACGGCGCGTGtgatga     >  2:238580/1‑148 (MQ=255)
                                                                                                                                        aTCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGCCAACGGCGCGTGTGATGaa    >  8:258589/1‑149 (MQ=255)
                                                                                                                                        aTCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGCCAACGGCGCGTGTGATGaa    >  4:125104/1‑149 (MQ=255)
                                                                                                                                         tCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGCCAACGGCGCGTGTGATGaa    <  7:296857/148‑1 (MQ=255)
                                                                                                                                         tCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGCCAACGGCGCGTGTGATGaa    >  8:296857/1‑148 (MQ=255)
                                                                                                                                           ctttctttTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGCCAACGGCGCGTGTGATGAATa  >  4:125433/1‑148 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GGCAGTGATAGTCAGTCACTATGGAGATCGCGGATGGTAACGTCCTGTACTGGACATGTTTTAGATAATCAACGCGCCACCACTCGCGGAGTTTTCTCCTCCGGTAGTCATCTCGTCACTTTGCATTTTCAACCTCATCCTTTCTTTTCATGTGTTACCGACGCCGTAAACGGTGCGCGTAGCCGTTTTTCGGCATTTTATCCAAAAGCAAACTATGGCTTACAAGGAAGCCAACCCTCTGATGTTCGTGCGCATAATCGCGCTGCCAACGGCGCGTGTGATGAATA  >  NC_000913/1084959‑1085245

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