Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,636,286 G→T 5.6% H357Q (CAC→CAA 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,2860GT 5.6% 612.7 / 17.5 171H357Q (CAC→CAAderGTPase; multicopy suppressor of ftsJ
Reads supporting (aligned to +/- strand):  ref base G (81/80);  new base T (5/5);  total (86/85)
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.83e-01

Only 100 of 171 total aligned reads displayed.

CTTCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGA  >  NC_000913/2636139‑2636430
                                                                                                                                                   |                                                                                                                                                
cTTCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGt                                                                                                                                                 >  1:39048/1‑149 (MQ=255)
 ttCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTg                                                                                                                                                >  2:853544/1‑149 (MQ=255)
    acacGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCaa                                                                                                                                             >  1:803926/1‑149 (MQ=255)
     cacGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAg                                                                                                                                            <  2:846472/149‑1 (MQ=255)
      acgatgACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAGCACCACTGCCGTGCAAg                                                                                                                                            <  4:63492/148‑1 (MQ=255)
       cgacgaCCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAgg                                                                                                                                           >  3:225623/1‑148 (MQ=255)
        gacgacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGc                                                                                                                                          <  4:595961/148‑1 (MQ=255)
        gacgacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGc                                                                                                                                          <  3:881078/148‑1 (MQ=255)
        gacgacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGGCCACACGCGACGGGTGGAGCTGTCAAACGCTTCACGTACTGATTCAAACAAGGTACCAACTCCACCGCCGTGCAAGGCa                                                                                                                                         <  2:902298/149‑1 (MQ=255)
          cgacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGAGCAAGGCag                                                                                                                                        <  3:67079/148‑1 (MQ=255)
            acCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCagag                                                                                                                                      >  3:13322/1‑148 (MQ=255)
              cgcgTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCTTGCACGGCagtggt                                                                                                                                    <  3:339697/148‑5 (MQ=255)
               gcgTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCTTGCACGGCagtggtgt                                                                                                                                  <  1:965019/149‑7 (MQ=255)
               gcgTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATa                                                                                                                                   >  4:1044350/1‑148 (MQ=255)
                cgTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATaaa                                                                                                                                 <  2:494142/149‑1 (MQ=255)
                  tACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGt                                                                                                                               <  1:76208/149‑1 (MQ=255)
                     cAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGt                                                                                                                                                 <  1:900960/128‑1 (MQ=255)
                     cAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGt                                                                                                                                                 >  2:900960/1‑128 (MQ=255)
                          gcggTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACg                                                                                                                        <  4:13322/148‑1 (MQ=255)
                          gcggTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGa                                                                                                                       <  2:284090/149‑1 (MQ=255)
                             gTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCa                                                                                                                    >  1:276469/1‑149 (MQ=255)
                                  tGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCag                                                                                                                                        <  3:759237/124‑1 (MQ=255)
                                  tGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATc                                                                                                               >  1:629752/1‑149 (MQ=255)
                                  tGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACAGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCag                                                                                                                                        >  4:759237/1‑124 (MQ=255)
                                    aTCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCTTGCACGGCagt                                                                                                                                       >  3:282730/1‑122 (MQ=255)
                                    aTCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCTTGCACGGCagt                                                                                                                                       <  4:282730/123‑2 (MQ=255)
                                    aTCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCg                                                                                                              <  4:339913/148‑1 (MQ=255)
                                    aTCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGa                                                                                                             <  2:803926/149‑1 (MQ=255)
                                     tCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAg                                                                                                                                            >  3:890290/1‑117 (MQ=255)
                                     tCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGa                                                                                                             >  3:859097/1‑148 (MQ=255)
                                       ttCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATaa                                                                                                          >  2:758748/1‑149 (MQ=255)
                                         cAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAg                                                                                                        <  1:1038785/149‑1 (MQ=255)
                                           aCAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGcc                                                                                                      >  2:952236/1‑149 (MQ=255)
                                                  tCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACg                                                                                                >  4:828849/1‑148 (MQ=255)
                                                    gTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAgg                                                                                                                                           >  2:22119/1‑103 (MQ=255)
                                                    gTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGaa                                                                                             >  1:534631/1‑149 (MQ=255)
                                                    gTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGaa                                                                                             <  1:904349/149‑1 (MQ=255)
                                                    gTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGa                                                                                              <  3:905872/148‑1 (MQ=255)
                                                    gTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGa                                                                                              <  4:701423/148‑1 (MQ=255)
                                                     tCATGAAGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAg                                                                                            >  2:209759/1‑149 (MQ=255)
                                                         gatgCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCa                                                                                        <  2:190156/149‑1 (MQ=255)
                                                          atgCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAg                                                                                       >  1:622925/1‑149 (MQ=255)
                                                                 tCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGtttct                                                                                <  2:629752/149‑1 (MQ=255)
                                                                 tCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGtttct                                                                                <  1:33289/149‑1 (MQ=255)
                                                                 tCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGtttc                                                                                 >  3:216169/1‑148 (MQ=255)
                                                                     cATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCa                                                                            <  2:314804/149‑1 (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)
                                                                               cccACACGACGAGTGGAGCTGTCATACGCTTCCCGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCttt                                                                  >  1:560999/1‑149 (MQ=255)
                                                                                ccACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCttt                                                                  >  3:529542/1‑148 (MQ=255)
                                                                                  acacGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCAc                                                               >  2:29980/1‑149 (MQ=255)
                                                                                  acacGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCAc                                                               >  1:637659/1‑149 (MQ=255)
                                                                                      gacgGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTcc                                                           <  2:733704/149‑1 (MQ=255)
                                                                                       acgGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTcc                                                           <  4:856805/148‑1 (MQ=255)
                                                                                          ggTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGAc                                                       <  1:480126/149‑1 (MQ=255)
                                                                                                     cATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATc                                             >  4:573085/1‑148 (MQ=255)
                                                                                                       tACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACtt                                                             >  1:155927/1‑130 (MQ=255)
                                                                                                       tACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATccc                                           >  4:427938/1‑148 (MQ=255)
                                                                                                        aCGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCa                                          >  3:279557/1‑148 (MQ=255)
                                                                                                        aCGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCa                                          <  4:225623/148‑1 (MQ=255)
                                                                                                         cGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCAc                                         >  3:587340/1‑148 (MQ=255)
                                                                                                           cTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGTACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACtt                                       <  3:935190/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                                     <  4:275592/148‑1 (MQ=255)
                                                                                                              cACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTc                                   <  1:567694/149‑1 (MQ=255)
                                                                                                               aCGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTca                                  <  2:622925/149‑1 (MQ=255)
                                                                                                               aCGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTca                                  <  2:1053478/149‑1 (MQ=255)
                                                                                                                      gATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACaa                            >  4:444298/1‑148 (MQ=255)
                                                                                                                           aaaCAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACa                       <  3:1005230/148‑1 (MQ=255)
                                                                                                                            aaCAAGTTACCCACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAg                     >  1:800272/1‑149 (MQ=255)
                                                                                                                            aaCAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAg                     >  2:540765/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)
                                                                                                                              cAAGTTACCAACAACACTGCCTTGCTAGGCAGACATAACGTCCACACGAGCAAAATCGTTAAAGCTGTGACGGAATTACAGCGGCACTTTCAACCGCTCTTTCACCTCCTAACTCAGGCCATCCCACTTATTCACCACAATGACAAGt                    >  4:16980/1‑148 (MQ=255)
                                                                                                                               aaGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGa                  <  2:333048/149‑1 (MQ=255)
                                                                                                                                aGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTg                   <  3:52668/147‑1 (MQ=255)
                                                                                                                                aGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTg                   >  4:52668/1‑147 (MQ=255)
                                                                                                                                 gTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAg                 <  3:573085/148‑1 (MQ=255)
                                                                                                                                  ttACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAgcg               >  2:461386/1‑149 (MQ=255)
                                                                                                                                  ttACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAgcg               >  1:476015/1‑149 (MQ=255)
                                                                                                                                      cAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCa           >  2:142634/1‑149 (MQ=255)
                                                                                                                                         caccacTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTa        <  2:167656/149‑1 (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)
                                                                                                                                           ccacTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTAt       <  3:242050/148‑1 (MQ=255)
                                                                                                                                            cacTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTg     >  1:627283/1‑149 (MQ=255)
                                                                                                                                             acTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTga    <  2:851296/149‑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)
                                                                                                                                                   |                                                                                                                                                
CTTCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGA  >  NC_000913/2636139‑2636430

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