Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,693,755 A→T 5.2% R494R (CGT→CGA bcsA ← cellulose synthase, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,693,7550AT 5.2% 336.5 / 7.4 97R494R (CGT→CGAbcsAcellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (47/44);  new base T (2/3);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01

CAGTCAGGAAGATCAGCCGTGGAATGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTCACGGTTTCGACAGCAATGCCGC  >  NC_000913/3693610‑3693897
                                                                                                                                                 |                                                                                                                                              
cAGTCAGGAAGATCAGCCGTGGAATGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGc                                                                                                                                             <  2:416497/149‑1 (MQ=255)
 aGTCAGGAAGATCAGCCGTGGAATGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGTCCAGCGAATACGCt                                                                                                                                            <  1:806805/149‑1 (MQ=255)
    cAGGAAGATCAGCCGTGGAATGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGa                                                                                                                                          <  3:1023129/148‑1 (MQ=255)
    cAGGAAGATCAGCCGTGGAATGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGAc                                                                                                                                         <  2:751134/149‑1 (MQ=255)
     aGGAAGATCAGCCGTGGAATGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGTATTCGCTggg                                                                                                                                         <  3:77517/148‑3 (MQ=255)
      ggAAGATCAGCCGTGGAATGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGAcc                                                                                                                                        >  4:495424/1‑148 (MQ=255)
        aaGATCAGCCGTGGAATGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGa                                                                                                                                      >  3:714370/1‑148 (MQ=255)
                ccGTGGAATGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACg                                                                                                                                              <  3:357380/132‑1 (MQ=255)
                ccGTGGAATGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACg                                                                                                                                              >  4:357380/1‑132 (MQ=255)
                 cGTGGAATGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGTATTCGCTg                                                                                                                                           >  3:1054311/1‑134 (MQ=255)
                 cGTGGAATGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGTATTCGCTg                                                                                                                                           <  4:1054311/134‑1 (MQ=255)
                 cGTGGAATGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCg                                                                                                                             >  4:988203/1‑148 (MQ=255)
                 cGTGGAATGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCAGCGCGCCCAGCGAATACGCTGACCGATATGCGCCg                                                                                                                             >  4:690045/1‑148 (MQ=255)
                      aaTGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGa                                                                                                                        <  3:721621/148‑1 (MQ=255)
                          cccGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACttt                                                                                                                    <  3:495424/148‑1 (MQ=255)
                          cccGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTc                                                                                                                   <  2:713805/149‑1 (MQ=255)
                          cccGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTc                                                                                                                   <  2:73633/149‑1 (MQ=255)
                           ccGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTc                                                                                                                   <  3:488289/148‑1 (MQ=255)
                           ccGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCg                                                                                                                  >  1:1647/1‑149 (MQ=255)
                           ccGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCg                                                                                                                  >  1:41168/1‑149 (MQ=255)
                           ccGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCg                                                                                                                  <  1:1076884/149‑1 (MQ=255)
                                   aagTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGTATTCGCTgg                                                                                                                                          >  2:172248/1‑116 (MQ=255)
                                   aagTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGTATTCGCTgg                                                                                                                                          <  1:172248/117‑2 (MQ=255)
                                   aagTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCa                                                                                                           >  4:734319/1‑148 (MQ=255)
                                   aagTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAg                                                                                                          >  2:253092/1‑149 (MQ=255)
                                         aaCATGGCGTCGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCccg                                                                                                     >  3:684047/1‑148 (MQ=255)
                                             tGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCcgccgc                                                                                                 >  4:606464/1‑148 (MQ=255)
                                              ggCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCgccgcc                                                                                                <  4:36734/148‑1 (MQ=255)
                                                 gTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCg                                                                                            <  1:352395/149‑1 (MQ=255)
                                                 gTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCg                                                                                            >  2:450399/1‑149 (MQ=255)
                                                  ttGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCg                                                                                            >  4:737741/1‑148 (MQ=255)
                                                  ttGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCg                                                                                            >  3:826898/1‑148 (MQ=255)
                                                  ttGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCg                                                                                            >  4:117586/1‑148 (MQ=255)
                                                    gatgTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGa                                                                                          >  4:165630/4‑148 (MQ=255)
                                                     aCGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAAt                                                                                        >  1:257024/1‑149 (MQ=255)
                                                       gTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATgcgc                                                                                                                               >  2:890209/1‑108 (MQ=255)
                                                       gTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATgcgc                                                                                                                               <  1:890209/108‑1 (MQ=255)
                                                       gTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATAc                                                                                      <  1:449727/149‑1 (MQ=255)
                                                         aaCATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGc                                                                                    <  1:940528/149‑1 (MQ=255)
                                                         aaCATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGc                                                                                    >  1:70129/1‑149 (MQ=255)
                                                            aTAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTc                                                                                                               <  1:135608/119‑1 (MQ=255)
                                                            aTAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTc                                                                                                               >  2:135608/1‑119 (MQ=255)
                                                            aTAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCata                                                                                 >  1:658425/1‑149 (MQ=255)
                                                             tAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCatat                                                                                >  1:997336/1‑149 (MQ=255)
                                                              aGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTGCCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCAtata                                                                               >  2:631337/1‑149 (MQ=255)
                                                              aGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCatat                                                                                <  3:734319/148‑1 (MQ=255)
                                                                 cGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATAcgc                                                                            <  2:84829/149‑1 (MQ=255)
                                                                    tGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCgg                                                                          <  3:165630/148‑1 (MQ=255)
                                                                    tGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGa                                                                         <  2:907546/149‑1 (MQ=255)
                                                                    tGAGCAAACTTCAGCCCTTTACCGGGGAGAGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCTGCCTGCGGAATACGCATATACGCGGa                                                                         <  2:1647/149‑1 (MQ=255)
                                                                       gCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAgg                                                                       <  4:714370/148‑1 (MQ=255)
                                                                       gCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGt                                                                      >  2:722849/1‑149 (MQ=255)
                                                                          accTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATgcg                                                                                                                                <  3:688022/86‑1 (MQ=255)
                                                                          aaCTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATgcg                                                                                                                                >  4:688022/1‑88 (MQ=255)
                                                                          aaCTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAAt                                                                                        <  4:870559/128‑1 (MQ=255)
                                                                          aaCTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAAt                                                                                        >  3:870559/1‑128 (MQ=255)
                                                                              tCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGcc                                                                >  3:656803/1‑148 (MQ=255)
                                                                                    cTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACgg                                                         >  1:897402/1‑149 (MQ=255)
                                                                                     tttACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGt                                                        <  2:61291/149‑1 (MQ=255)
                                                                                                gCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCag                                              <  4:166038/148‑1 (MQ=255)
                                                                                                 cGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCaga                                             >  4:188076/1‑148 (MQ=255)
                                                                                                   ggTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGaa                                          <  1:762543/149‑1 (MQ=255)
                                                                                                    gTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAg                                         >  1:283331/1‑149 (MQ=255)
                                                                                                     ttATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATAc                                                                                      <  1:779211/103‑1 (MQ=255)
                                                                                                     ttATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATAc                                                                                      >  2:779211/1‑103 (MQ=255)
                                                                                                     ttATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAg                                         >  3:826652/1‑148 (MQ=255)
                                                                                                     ttATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCAATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGta                                                                     <  4:174257/120‑1 (MQ=255)
                                                                                                     ttATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCAATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGta                                                                     >  3:174257/1‑120 (MQ=255)
                                                                                                         cGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATgc                                    <  2:257024/149‑1 (MQ=255)
                                                                                                         cGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATg                                     >  4:36904/1‑148 (MQ=255)
                                                                                                              cGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCAt                                <  3:606464/148‑1 (MQ=255)
                                                                                                               ggAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCt                              >  1:625959/1‑149 (MQ=255)
                                                                                                                gAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCa                                                                                   >  2:318612/1‑95 (MQ=255)
                                                                                                                gAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCa                                                                                   <  1:318612/95‑1 (MQ=255)
                                                                                                                gAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCtt                             >  1:804693/1‑149 (MQ=255)
                                                                                                                gAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCtt                             >  2:401067/1‑149 (MQ=255)
                                                                                                                 aaGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTc                            >  2:179026/1‑149 (MQ=255)
                                                                                                                 aaGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTc                            >  1:498829/1‑149 (MQ=255)
                                                                                                                  aGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCa                           <  1:722849/149‑1 (MQ=255)
                                                                                                                   gATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCa                           <  3:36904/148‑1 (MQ=255)
                                                                                                                      ttGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTc                        >  4:543243/1‑148 (MQ=255)
                                                                                                                           ccATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCaa                                                    >  4:26598/1‑115 (MQ=255)
                                                                                                                           ccATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCaa                                                    <  3:26598/115‑1 (MQ=255)
                                                                                                                            cATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTCACGGttt                 <  2:737280/149‑1 (MQ=255)
                                                                                                                             aTCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTCACGGttt                 <  3:988203/148‑1 (MQ=255)
                                                                                                                              tCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTCACGGTTTc                >  4:894654/1‑148 (MQ=255)
                                                                                                                               cccGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTCa                       <  3:215017/140‑1 (MQ=255)
                                                                                                                               cccGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTCa                       >  4:215017/1‑140 (MQ=255)
                                                                                                                               cccGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTCACGGTTTCg               <  3:655374/148‑1 (MQ=255)
                                                                                                                                  gcgcgcCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTCACGGTTTCGACAg           >  1:854272/1‑149 (MQ=255)
                                                                                                                                      gcCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTCACGGTTTCGACAGCAAt       <  1:179026/149‑1 (MQ=255)
                                                                                                                                        ccAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTCACGGTTTCGACAGCAATgc     <  2:679440/149‑1 (MQ=255)
                                                                                                                                        ccAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTCACGGTTTCGACAGCAATgc     >  1:777008/1‑149 (MQ=255)
                                                                                                                                          aGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTCACGGTTTCGACAGCAATgccg   <  2:498829/149‑1 (MQ=255)
                                                                                                                                           gCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTCACGGTTTCGACAGCAATgccgc  <  2:658425/149‑1 (MQ=255)
                                                                                                                                            cGAATCCGCTGACCGATATGCGTCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCAACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCCGTCACGGTTTCGACAGCAATgccga  <  4:570642/148‑2 (MQ=255)
                                                                                                                                            cGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTCACGGTTTCGACAGCAATgccgc  >  3:467162/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CAGTCAGGAAGATCAGCCGTGGAATGCCCGACAAGAAGTGGAACATGGCGTTGACGTAACATAGCCGCTGAGCAAACTTCAGCCCTTTACCGGTGAGCGGGTTATCGAGACGGAAGATTTGTACCATCCCGCGCGCCCAGCGAATACGCTGACCGATATGCGCCGACAGACTTTCGGTCGCCAGCCCCGCCGCCTGCGGAATACGCATATACGCGGAGGTATAGCCACGACGGTGCAACCGCAGAGAAGTATGCGCATCTTCAGTCACGGTTTCGACAGCAATGCCGC  >  NC_000913/3693610‑3693897

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