Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,327,471 C→T 5.2% P113S (CCC→TCC)  crfC → clamp‑binding sister replication fork colocalization protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,327,4710CT 5.2% 241.3 / 2.6 77P113S (CCC→TCC) crfCclamp‑binding sister replication fork colocalization protein
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base T (2/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATAGATAAAGATATGCGTGCGCTTATGCAG  >  NC_000913/4327322‑4327620
                                                                                                                                                     |                                                                                                                                                     
gcaccaTTTCCAGTCGGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCCCCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGCCTGCGCTGCCGACGCTTATTCGCCATACGc                                                                                                                                                       >  1:913101/6‑150 (MQ=255)
 tAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGcc                                                                                                                                                      <  2:31124/150‑1 (MQ=255)
    aaTTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCgg                                                                                                                                                   <  1:1392655/150‑1 (MQ=255)
     aTTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCggg                                                                                                                                                  <  2:1376658/150‑1 (MQ=255)
         cccGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCaaa                                                                                                                                              >  2:1010739/1‑150 (MQ=255)
           cGACTGGAAATGGTGCTTGCAATTGTCGGTACAATGAAAGCAGGGAAATCAACCAACATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGAAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAg                                                                                                                                            >  2:1562339/1‑150 (MQ=255)
               tGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAAc                                                                                                                                        >  1:1284010/1‑150 (MQ=255)
                ggAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAAcc                                                                                                                                       >  1:919685/1‑150 (MQ=255)
                ggAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACCCCAGGGCAAAAGGacac                                                                                                                                       >  1:688077/1‑147 (MQ=255)
                 gAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCg                                                                                                                                      >  2:714136/1‑150 (MQ=255)
                     tGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTAc                                                                                                                                  <  2:372922/150‑1 (MQ=255)
                       gTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTg                                                                                                                                <  1:1445466/150‑1 (MQ=255)
                        tGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCAGCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACTCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGc                                                                                                                               <  2:763268/150‑1 (MQ=255)
                        tGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGc                                                                                                                               <  2:919685/150‑1 (MQ=255)
                        tGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGc                                                                                                                               >  1:163827/1‑150 (MQ=255)
                              cAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCAttttt                                                                                                                         >  2:1414179/1‑150 (MQ=255)
                               aaTTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTAGGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTc                                                                                                                        <  1:912383/150‑1 (MQ=255)
                               aaTTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTc                                                                                                                        <  1:5831/150‑1 (MQ=255)
                                 ttGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTcac                                                                                                                      <  2:1357414/150‑1 (MQ=255)
                                  tGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTcaca                                                                                                                     >  2:787097/1‑150 (MQ=255)
                                    tCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATg                                                                                                                   <  1:322324/150‑1 (MQ=255)
                                     cGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGt                                                                                                                  <  2:683226/150‑1 (MQ=255)
                                      ggTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTc                                                                                                                 >  1:126305/1‑150 (MQ=255)
                                        tACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTcgc                                                                                                               >  1:1395671/1‑150 (MQ=255)
                                           cATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGcc                                                                                                            <  2:1319668/150‑1 (MQ=255)
                                           cATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGcc                                                                                                            >  2:742494/1‑150 (MQ=255)
                                            aTGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCa                                                                                                           >  1:1250452/1‑150 (MQ=255)
                                               aaaGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATc                                                                                                        >  2:1335253/1‑150 (MQ=255)
                                                    aGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTg                                                                                                   >  2:272548/1‑150 (MQ=255)
                                                       gAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGttt                                                                                                <  1:987805/150‑1 (MQ=255)
                                                        aaaTCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTa                                                                                               <  2:1482324/150‑1 (MQ=255)
                                                         aaTCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTaa                                                                                              <  1:695930/150‑1 (MQ=255)
                                                          aTCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGTCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAAt                                                                                             >  1:542061/1‑150 (MQ=255)
                                                          aTCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGTCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAAt                                                                                             >  1:536518/1‑150 (MQ=255)
                                                                caccaTTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACTCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTcaaca                                                                                       <  1:1485608/150‑1 (MQ=255)
                                                                  ccaTTTATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTGTAATTCAACAAt                                                                                     >  2:838306/1‑150 (MQ=255)
                                                                  ccaTTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAAt                                                                                     <  2:1221485/150‑1 (MQ=255)
                                                                  ccaTTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAAt                                                                                     <  2:126305/150‑1 (MQ=255)
                                                                      tAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTaca                                                                                 <  1:1509987/150‑1 (MQ=255)
                                                                         tGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTCCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTacaaca                                                                              <  1:1285774/150‑1 (MQ=255)
                                                                               tgttgGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCt                                                                        >  2:756210/1‑150 (MQ=255)
                                                                                gttgGTACGGAGGTTCTGCGTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTg                                                                       <  1:664262/150‑1 (MQ=255)
                                                                                   ggTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGt                                                                    >  1:982003/1‑150 (MQ=255)
                                                                                     tACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGa                                                                  >  2:120714/1‑150 (MQ=255)
                                                                                       cGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGAtt                                                                >  2:1542597/1‑150 (MQ=255)
                                                                                            gTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGat                                                           <  2:1284010/150‑1 (MQ=255)
                                                                                               cTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATAtt                                                        <  2:186716/150‑1 (MQ=255)
                                                                                                tGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTa                                                       <  1:1010739/150‑1 (MQ=255)
                                                                                                 gCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTaa                                                      >  2:197249/1‑150 (MQ=255)
                                                                                                  acTAATCGTAATCGCCCAATGAATGCGCTGCCGACGCTTATTCGCCATACGACCCGGCAAAACGAACCGGTACTCCAATTTTAACATGTCGCGCCAATCGATTGTTTAATTCAAAAATTAAAACAGAGCCTGCGTGATTGAGATATTaac                                                     >  2:1310441/2‑149 (MQ=255)
                                                                                                   cTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGc                                                    >  2:498624/1‑150 (MQ=255)
                                                                                                   cTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGc                                                    <  1:479643/150‑1 (MQ=255)
                                                                                                    taatcgtaatcgCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCa                                                   >  1:180306/1‑150 (MQ=255)
                                                                                                      atcgtaatcgCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATc                                                 >  1:87295/1‑150 (MQ=255)
                                                                                                      atcgtaatcgCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATc                                                 >  2:750044/1‑150 (MQ=255)
                                                                                                            atcgCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCg                                           <  1:1414179/150‑1 (MQ=255)
                                                                                                             tcgCCCAATGACTGAGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGAGTGAGTGCGATATTAAGGCTCTGACCgg                                          >  1:1061904/1‑149 (MQ=255)
                                                                                                               gCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGAtg                                        <  2:1395671/150‑1 (MQ=255)
                                                                                                               gCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGAtg                                        >  1:1308767/1‑150 (MQ=255)
                                                                                                                 ccAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGAtgtg                                      <  1:1335253/150‑1 (MQ=255)
                                                                                                                   aaTGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCt                                    <  2:127872/150‑1 (MQ=255)
                                                                                                                    aTGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTg                                   <  1:9118/150‑1 (MQ=255)
                                                                                                                        cTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGaaa                               >  2:565839/1‑150 (MQ=255)
                                                                                                                          gcgcTGCCGACGCTTATTCTCCATACGCCCGGGCAAAAGTAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATa                             <  1:1562339/150‑1 (MQ=255)
                                                                                                                          gcgcTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATa                             >  2:545118/1‑150 (MQ=255)
                                                                                                                              tGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATAGATa                         <  1:272548/150‑1 (MQ=255)
                                                                                                                               gCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATAGATaa                        >  1:1106740/1‑150 (MQ=255)
                                                                                                                               gCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGAGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATAGATaa                        >  2:760972/1‑150 (MQ=255)
                                                                                                                                  gACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATAGATAAAGa                     >  1:464549/1‑150 (MQ=255)
                                                                                                                                   aCGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATAGATAAAGat                    <  2:464549/150‑1 (MQ=255)
                                                                                                                                   aCGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATAGATAAAGat                    <  2:1106740/150‑1 (MQ=255)
                                                                                                                                    ggCTTATTCGCCATACGTCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAAGAGCGCCTGCGTGATTGCGATATTAAGCATCTGAGCGATGTGCTGGAAATAGATAAAGata                   <  2:542061/149‑1 (MQ=255)
                                                                                                                                    cGCTTATTCGCCATACGTCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATAGATAAAGata                   <  2:536518/150‑1 (MQ=255)
                                                                                                                                        tATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATAGATAAAGATAtgcg               >  2:318248/1‑150 (MQ=255)
                                                                                                                                            cGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTg                                               <  1:19767/114‑1 (MQ=255)
                                                                                                                                            cGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTg                                               >  2:19767/1‑114 (MQ=255)
                                                                                                                                                 tACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATAGATAAAGATATGCGTGCGCTTAt      >  1:435853/1‑150 (MQ=255)
                                                                                                                                                 tACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATAGATAAAGATATGCGTGCGCTTAt      <  1:1358806/150‑1 (MQ=255)
                                                                                                                                                   cGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATAGATAAAGATATGCGTGCGCTTATGc    >  2:1391937/1‑150 (MQ=255)
                                                                                                                                                     cccGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATAGATAAAGATATGCGTGCGCTTATGCAg  >  1:682888/1‑150 (MQ=255)
                                                                                                                                                     cccGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATAGATAAAGATATGCGTGCGCTTATGCAg  >  2:1036360/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAAAATTTCCCGACTGGAAATGGTGCTTGCAATTGTCGGTACCATGAAAGCAGGGAAATCAACCACCATTAATGCCATTGTTGGTACGGAGGTTCTGCCTAATCGTAATCGCCCAATGACTGCGCTGCCGACGCTTATTCGCCATACGCCCGGGCAAAAGGAACCGGTACTGCATTTTTCACATGTCGCGCCAATCGATTGTTTAATTCAACAATTACAACAGCGCCTGCGTGATTGCGATATTAAGCATCTGACCGATGTGCTGGAAATAGATAAAGATATGCGTGCGCTTATGCAG  >  NC_000913/4327322‑4327620

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