Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,875,413 T→C T3A (ACC→GCC)  cysN ← sulfate adenylyltransferase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,875,4130TC100.0% 1034.3 / NA 240T3A (ACC→GCC) cysNsulfate adenylyltransferase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (129/111);  total (129/111)

Only 100 of 249 total aligned reads displayed.

GCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAAC  >  NC_000913/2875315‑2875510
                                                                                                  |                                                                                                 
gCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGgcgt                                                                                                 <  2:4564071/101‑1 (MQ=255)
gCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGgcgt                                                                                                 >  1:4564392/1‑101 (MQ=255)
 cTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGtt                                                                                                >  1:4564395/1‑101 (MQ=255)
  tACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTc                                                                                               <  2:4563816/101‑1 (MQ=255)
   aCCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCa                                                                                              >  2:4564399/1‑101 (MQ=255)
   aCCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCa                                                                                              >  1:4564397/1‑101 (MQ=255)
    ccACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCAt                                                                                             <  1:4563766/101‑1 (MQ=255)
    ccACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCAt                                                                                             <  1:4563955/101‑1 (MQ=255)
    ccACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCAt                                                                                             >  2:4564403/1‑101 (MQ=255)
     cacaGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATc                                                                                            <  2:4564026/101‑1 (MQ=255)
        aGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTa                                                                                         >  1:4564407/1‑101 (MQ=255)
          gTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTaaa                                                                                       >  2:4564409/1‑101 (MQ=255)
           tCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTaaaa                                                                                      <  1:4564015/101‑1 (MQ=255)
           tCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTaaaa                                                                                      <  2:4563988/101‑1 (MQ=255)
            cAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTaaaaa                                                                                     >  1:4564411/1‑101 (MQ=255)
             aGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAAt                                                                                    >  2:4564414/1‑101 (MQ=255)
             aGAAAACGCAGCAGGCTATTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTaa                                                                                        <  1:4564177/97‑1 (MQ=255)
              gAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATa                                                                                   <  1:4563821/101‑1 (MQ=255)
                aaaCGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATAcc                                                                                 <  2:4563948/101‑1 (MQ=255)
                 aaCGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATAccc                                                                                >  2:4564423/1‑101 (MQ=255)
                 aaCGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATAccc                                                                                >  1:4564419/1‑101 (MQ=255)
                 aaCGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATAccc                                                                                >  1:4564421/1‑101 (MQ=255)
                  aCGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATAcccc                                                                               >  2:4564426/1‑101 (MQ=255)
                   cGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCt                                                                              >  1:4564427/1‑101 (MQ=255)
                     cagcagGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATCCCCCTGa                                                                            <  1:4563870/101‑1 (MQ=255)
                      agcagGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGAc                                                                           <  1:4563911/101‑1 (MQ=255)
                      agcagGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGAc                                                                           >  1:4564430/1‑101 (MQ=255)
                        cagGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGt                                                                         >  1:4564434/1‑101 (MQ=255)
                         agGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGtt                                                                        <  2:4564050/101‑1 (MQ=255)
                          ggCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGttt                                                                       >  1:4564436/1‑101 (MQ=255)
                           gCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGtttt                                                                      >  2:4564438/1‑101 (MQ=255)
                           gCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGtttt                                                                      <  2:4564063/101‑1 (MQ=255)
                            cTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGttttt                                                                     <  1:4564227/101‑1 (MQ=255)
                            cTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGttttt                                                                     <  2:4564300/101‑1 (MQ=255)
                             ttttATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGtttttt                                                                    >  2:4564440/1‑101 (MQ=255)
                             ttttATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGtttttt                                                                    <  2:4564091/101‑1 (MQ=255)
                              tttATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTc                                                                   >  1:4564444/1‑101 (MQ=255)
                               ttATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCa                                                                  <  2:4564180/101‑1 (MQ=255)
                               ttATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCa                                                                  >  1:4564446/1‑101 (MQ=255)
                                tATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAg                                                                 <  1:4564074/101‑1 (MQ=255)
                                   gttgttGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTc                                                              >  2:4564451/1‑101 (MQ=255)
                                   gttgttGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTc                                                              >  2:4564449/1‑101 (MQ=255)
                                   gttgttGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTc                                                              <  2:4563973/101‑1 (MQ=255)
                                   gttgttGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTc                                                              <  1:4564100/101‑1 (MQ=255)
                                   gttgttGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTc                                                              <  1:4564004/101‑1 (MQ=255)
                                    ttgttgCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTcc                                                             <  2:4564120/101‑1 (MQ=255)
                                     tgttgCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATCCCCCTGACGTTTTTTCAGCTCCa                                                            <  1:4563993/101‑1 (MQ=255)
                                      gttgCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCAt                                                           >  1:4564454/1‑101 (MQ=255)
                                       ttgCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATa                                                          <  1:4564161/101‑1 (MQ=255)
                                       ttgCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATa                                                          <  1:4564023/101‑1 (MQ=255)
                                        tgCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAg                                                         >  1:4564456/1‑101 (MQ=255)
                                         gcgcAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGa                                                        <  2:4564080/101‑1 (MQ=255)
                                          cgcAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGAc                                                       >  2:4564458/1‑101 (MQ=255)
                                           gcAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGAcc                                                      <  1:4564090/101‑1 (MQ=255)
                                           gcAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGAcc                                                      >  2:4564460/1‑101 (MQ=255)
                                              atcatcCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGAc                                                       >  1:4564464/1‑97 (MQ=255)
                                              atcatcCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCg                                                   <  2:4563761/101‑1 (MQ=255)
                                                catcCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGcc                                                 >  1:4564470/1‑101 (MQ=255)
                                                catcCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGcc                                                 >  2:4564468/1‑101 (MQ=255)
                                                catcCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGcc                                                 >  2:4564466/1‑101 (MQ=255)
                                                 atcCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCt                                                >  2:4564473/1‑101 (MQ=255)
                                                  tcCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATTCCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTg                                               <  2:4564093/101‑1 (MQ=255)
                                                  tcCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGAcc                                                      <  2:4564109/94‑1 (MQ=255)
                                                     aGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTc                                            >  1:4564477/1‑101 (MQ=255)
                                                     aGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTc                                            <  1:4564048/101‑1 (MQ=255)
                                                      ggCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTcg                                           <  2:4564218/101‑1 (MQ=255)
                                                        cTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTcgcg                                         >  2:4564480/1‑101 (MQ=255)
                                                         ttCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCgg                                        >  1:4564482/1‑101 (MQ=255)
                                                          tCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCt                                                               >  2:4564486/1‑77 (MQ=255)
                                                            gACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCa                                     <  2:4564037/101‑1 (MQ=255)
                                                            gACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCa                                     >  1:4564489/1‑101 (MQ=255)
                                                            gACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCa                                     >  2:4564490/1‑101 (MQ=255)
                                                              cgccgcCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAAt                                   <  2:4564106/101‑1 (MQ=255)
                                                              cgccgcCTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAAt                                   >  2:4564494/1‑101 (MQ=255)
                                                                ccgccTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCa                                 >  1:4564495/1‑101 (MQ=255)
                                                                   ccTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCAcgc                              <  1:4563934/101‑1 (MQ=255)
                                                                    cTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCAcgcg                             <  1:4564279/101‑1 (MQ=255)
                                                                    cTTCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCAcgcg                             >  2:4564499/1‑101 (MQ=255)
                                                                     ttCATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCgg                            >  2:4564502/1‑101 (MQ=255)
                                                                       cATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGcc                          >  1:4564504/1‑101 (MQ=255)
                                                                       cATTGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGcc                          >  2:4564505/1‑101 (MQ=255)
                                                                          tGGCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTg                       <  1:4564209/101‑1 (MQ=255)
                                                                          tGGCGATTTGTTGAGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTg                       >  1:4564509/1‑101 (MQ=255)
                                                                           ggCGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGa                      >  2:4564510/1‑101 (MQ=255)
                                                                             cGATTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATCCCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACg                    <  2:4564057/101‑1 (MQ=255)
                                                                               aTTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGtt                  <  2:4564164/101‑1 (MQ=255)
                                                                               aTTTGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGtt                  <  2:4564066/101‑1 (MQ=255)
                                                                                tttGTTGTGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTc                 >  2:4564512/1‑101 (MQ=255)
                                                                                   gttgtGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACt              <  1:4563946/101‑1 (MQ=255)
                                                                                   gttgtGCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACt              >  2:4564516/1‑101 (MQ=255)
                                                                                       tgCAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACtggtg          <  1:4564152/101‑1 (MQ=255)
                                                                                         cAAGTGCGGCGTTCATCTTAAAAATCCCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACtggtggt        <  2:4564149/101‑1 (MQ=255)
                                                                                         cAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACtggtggt        <  1:4564190/101‑1 (MQ=255)
                                                                                         cAAGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACtggtggt        >  2:4564519/1‑101 (MQ=255)
                                                                                           aGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACtggtggtgg      >  1:4564522/1‑101 (MQ=255)
                                                                                           aGTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACtggtggtgg      >  2:4564523/1‑101 (MQ=255)
                                                                                            gTGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGa     >  2:4564525/1‑101 (MQ=255)
                                                                                             tGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGaa    <  2:4564135/101‑1 (MQ=255)
                                                                                             tGCGGCGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGaa    >  1:4564527/1‑101 (MQ=255)
                                                                                              gcggcgTTCATCTTAAAAATCCCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGaaa   <  2:4564198/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
GCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAAC  >  NC_000913/2875315‑2875510

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