Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,636,291 A→C 5.4% L356V (TTG→GTG)  der ← GTPase; multicopy suppressor of ftsJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,636,2910AC 5.4% 585.7 / 22.8 168L356V (TTG→GTG) derGTPase; multicopy suppressor of ftsJ
Reads supporting (aligned to +/- strand):  ref base A (78/81);  new base C (4/5);  total (82/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.03e-01

Only 100 of 170 total aligned reads displayed.

ACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAA  >  NC_000913/2636143‑2636435
                                                                                                                                                    |                                                                                                                                                
acacGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCaa                                                                                                                                                  >  1:803926/1‑149 (MQ=255)
 cacGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCaa                                                                                                                                                  >  3:859079/1‑148 (MQ=255)
    gacgacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGc                                                                                                                                               <  4:595961/148‑1 (MQ=255)
    gacgacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGGCCACACGCGACGGGTGGAGCTGTCAAACGCTTCACGTACTGATTCAAACAAGGTACCAACTCCACCGCCGTGCAAGGCa                                                                                                                                              <  2:902298/149‑1 (MQ=255)
     acgacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCa                                                                                                                                              >  3:275592/1‑148 (MQ=255)
      cgacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCag                                                                                                                                             <  3:209522/148‑1 (MQ=255)
      cgacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGAGCAAGGCag                                                                                                                                             <  3:67079/148‑1 (MQ=255)
       gacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCTTGCACGGCagtg                                                                                                                                           <  2:824597/149‑3 (MQ=255)
       gacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCagag                                                                                                                                           >  2:198783/1‑149 (MQ=255)
          cgcgTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCTTGCACGGCagtggt                                                                                                                                         <  3:339697/148‑5 (MQ=255)
           gcgTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATa                                                                                                                                        >  4:1044350/1‑148 (MQ=255)
            cgTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATaaa                                                                                                                                      <  1:975305/149‑1 (MQ=255)
            cgTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATaa                                                                                                                                       <  4:146101/148‑1 (MQ=255)
             gTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATaaa                                                                                                                                      <  3:868027/148‑1 (MQ=255)
              tACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGt                                                                                                                                    <  1:76208/149‑1 (MQ=255)
                        ggTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGc                                                                                                                          >  2:743938/1‑149 (MQ=255)
                         gTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCa                                                                                                                         >  1:276469/1‑149 (MQ=255)
                              tGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACAGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCag                                                                                                                                             >  4:759237/1‑124 (MQ=255)
                               gATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATc                                                                                                                    >  3:125602/1‑148 (MQ=255)
                                aTCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCTTGCACGGCagt                                                                                                                                            >  3:282730/1‑122 (MQ=255)
                                aTCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCTTGCACGGCagt                                                                                                                                            <  4:282730/123‑2 (MQ=255)
                                aTCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGa                                                                                                                  <  2:803926/149‑1 (MQ=255)
                                 tCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAg                                                                                                                                                 >  3:890290/1‑117 (MQ=255)
                                   ttCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATaa                                                                                                               >  1:39307/1‑149 (MQ=255)
                                     cAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAg                                                                                                             <  1:1038785/149‑1 (MQ=255)
                                      aaCAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGc                                                                                                            >  2:28198/1‑149 (MQ=255)
                                       aCAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGcc                                                                                                           >  2:952236/1‑149 (MQ=255)
                                              tCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACg                                                                                                     >  4:828849/1‑148 (MQ=255)
                                                gTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGaa                                                                                                  <  1:904349/149‑1 (MQ=255)
                                                gTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGaa                                                                                                  >  1:534631/1‑149 (MQ=255)
                                                gTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGa                                                                                                   <  4:701423/148‑1 (MQ=255)
                                                gTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGa                                                                                                   <  3:905872/148‑1 (MQ=255)
                                                 tCATGAAGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAg                                                                                                 >  2:209759/1‑149 (MQ=255)
                                                    tgatgCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTcc                                                                                              <  2:534631/149‑1 (MQ=255)
                                                    tgatgCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTcc                                                                                              <  2:886810/149‑1 (MQ=255)
                                                      atgCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAg                                                                                            >  1:622925/1‑149 (MQ=255)
                                                             tCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGtttct                                                                                     <  1:33289/149‑1 (MQ=255)
                                                             tCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGtttc                                                                                      <  3:91491/148‑1 (MQ=255)
                                                             tCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGtttc                                                                                      >  3:216169/1‑148 (MQ=255)
                                                                  aTAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCa                                                                                 >  4:750628/1‑148 (MQ=255)
                                                                     gagGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAg                                                                                            <  3:72662/134‑1 (MQ=255)
                                                                     gagGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAg                                                                                            >  4:72662/1‑134 (MQ=255)
                                                                     gagGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCt                                                                              >  3:980755/1‑148 (MQ=255)
                                                                     gagGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTg                                                                             >  2:1068151/1‑149 (MQ=255)
                                                                       ggTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGct                                                                           <  2:39048/149‑1 (MQ=255)
                                                                        gTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGctc                                                                          >  2:850590/1‑149 (MQ=255)
                                                                           cccACACGACGAGTGGAGCTGTCATACGCTTCCCGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCttt                                                                       >  1:560999/1‑149 (MQ=255)
                                                                            ccACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCttt                                                                       >  3:529542/1‑148 (MQ=255)
                                                                            ccACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTc                                                                      <  1:901687/149‑1 (MQ=255)
                                                                              acacGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCAc                                                                    >  1:637659/1‑149 (MQ=255)
                                                                              acacGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCAc                                                                    >  2:29980/1‑149 (MQ=255)
                                                                                acgacgGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACtt                                                                  <  2:626796/149‑1 (MQ=255)
                                                                                  gacgGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTcc                                                                <  2:733704/149‑1 (MQ=255)
                                                                                   acgGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTcc                                                                <  4:856805/148‑1 (MQ=255)
                                                                                     gggTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGa                                                             >  1:627835/1‑149 (MQ=255)
                                                                                      ggTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGAc                                                            <  1:480126/149‑1 (MQ=255)
                                                                                                 cATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATc                                                  >  4:573085/1‑148 (MQ=255)
                                                                                                  aTACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATccc                                                >  2:567694/1‑149 (MQ=255)
                                                                                                   tACGCTTCACGTACTGATTCAAACAAGTTACCAACGCCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATccc                                                >  3:510174/1‑148 (MQ=255)
                                                                                                   tACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACtt                                                                  >  1:155927/1‑130 (MQ=255)
                                                                                                   tACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACtt                                                                  <  2:155927/130‑1 (MQ=255)
                                                                                                   tACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATccc                                                >  4:427938/1‑148 (MQ=255)
                                                                                                    aCGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCaaa                                                                                                                       <  4:680374/76‑1 (MQ=255)
                                                                                                    aCGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCa                                               <  4:225623/148‑1 (MQ=255)
                                                                                                       cTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACtt                                            <  3:427938/148‑1 (MQ=255)
                                                                                                       cTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACtt                                            <  3:8025/148‑1 (MQ=255)
                                                                                                       cTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTa                                           >  1:219860/1‑149 (MQ=255)
                                                                                                        ttCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTAt                                          >  1:1058206/1‑149 (MQ=255)
                                                                                                         tCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTAt                                          <  3:735537/148‑1 (MQ=255)
                                                                                                         tCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTAt                                          <  4:275592/148‑1 (MQ=255)
                                                                                                          cACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTc                                        <  1:567694/149‑1 (MQ=255)
                                                                                                             gTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCa                                                    >  3:57651/1‑134 (MQ=255)
                                                                                                                 tGATTCAAACAAGTTACCAACACCACTGACGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACaa                                 <  2:286620/149‑1 (MQ=255)
                                                                                                                  gATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACaa                                 >  4:444298/1‑148 (MQ=255)
                                                                                                                        aaCAAGTTACCCACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAg                          >  1:800272/1‑149 (MQ=255)
                                                                                                                        aaCAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAg                          >  1:167656/1‑149 (MQ=255)
                                                                                                                         aCAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAg                          <  3:828849/148‑1 (MQ=255)
                                                                                                                         aCAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAg                          >  4:810830/1‑148 (MQ=255)
                                                                                                                            aGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTg                        >  4:52668/1‑147 (MQ=255)
                                                                                                                             gTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAgc                     <  2:1058206/149‑1 (MQ=255)
                                                                                                                             gTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAg                      <  3:573085/148‑1 (MQ=255)
                                                                                                                              ttACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAgcg                    >  2:461386/1‑149 (MQ=255)
                                                                                                                              ttACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAgcg                    >  1:476015/1‑149 (MQ=255)
                                                                                                                               tACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAgcg                    >  4:105380/1‑148 (MQ=255)
                                                                                                                               tACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAgcg                    >  3:24165/1‑148 (MQ=255)
                                                                                                                                     caccacTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTa             <  2:637659/149‑1 (MQ=255)
                                                                                                                                      accacTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTAt            >  1:1025584/1‑149 (MQ=255)
                                                                                                                                      accacTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTGCACCTGCTCTTACACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTAt            >  2:767364/1‑149 (MQ=255)
                                                                                                                                       tctcTGCCTTGCACGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTAt            <  4:904142/145‑1 (MQ=255)
                                                                                                                                        ctcTGCCTTGCACGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTg          >  2:315068/3‑149 (MQ=255)
                                                                                                                                        ctcTGCCTTGCACGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTg          >  2:18379/3‑149 (MQ=255)
                                                                                                                                        cacTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTg          >  1:627283/1‑149 (MQ=255)
                                                                                                                                         ccTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAATATCGATAACGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACACGTGAGCGCCCACTATTg          <  4:378452/147‑1 (MQ=255)
                                                                                                                                         acTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTga         <  1:853544/149‑1 (MQ=255)
                                                                                                                                         acTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTga         <  1:224604/149‑1 (MQ=255)
                                                                                                                                         acTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTg          <  4:529542/148‑1 (MQ=255)
                                                                                                                                            gCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTgaga       <  4:859097/148‑1 (MQ=255)
                                                                                                                                                tGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATaaa  <  2:627835/149‑1 (MQ=255)
                                                                                                                                                 gCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTGTTCACCACAATGACAAGTGAGCGCCCACTATTgaga       >  4:248019/1‑143 (MQ=255)
                                                                                                                                                 gCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTGTTCACCACAATGACAAGTGAGCGCCCACTATTgaga       <  3:248019/143‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAA  >  NC_000913/2636143‑2636435

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: