Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,688,074 G→A 5.2% P367L (CCT→CTT)  bcsC ← cellulose synthase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,688,0740GA 5.2% 251.3 / 3.3 77P367L (CCT→CTT) bcsCcellulose synthase subunit
Reads supporting (aligned to +/- strand):  ref base G (33/40);  new base A (2/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

CATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑GGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGAT  >  NC_000913/3687932‑3688221
                                                                                                                                              |                                                                                                                                                    
cATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGattat                                                                                                                                                >  4:145273/1‑149 (MQ=255)
 aTACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGattatt                                                                                                                                               <  2:565640/149‑1 (MQ=255)
 aTACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGattatt                                                                                                                                               <  3:202090/149‑1 (MQ=255)
  tACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTg                                                                                                                                              <  1:316769/149‑1 (MQ=255)
  tACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTg                                                                                                                                              >  4:27355/1‑149 (MQ=255)
   aCGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTgg                                                                                                                                             <  4:66723/149‑1 (MQ=255)
    cGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGc                                                                                                                                            <  2:868292/149‑1 (MQ=255)
     gCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCt                                                                                                                                           >  1:111750/1‑149 (MQ=255)
     gCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCt                                                                                                                                           <  4:331672/149‑1 (MQ=255)
         ggTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCa                                                                                                                                       <  1:323821/149‑1 (MQ=255)
            ctgctgATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAgcg                                                                                                                                    <  3:357104/149‑1 (MQ=255)
              gctgATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCgca                                                                                                                                  <  4:277401/149‑1 (MQ=255)
              gctgATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCgca                                                                                                                                  >  1:679210/1‑149 (MQ=255)
                tgATAATAAAATTCGTCGGCGGGATAATCTTTTCGCGCCATCGCCTCTTCGCCCAGCCCTAGCAGTTCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACACGCGTTCTGCCCGGTCAGGCTCATTGGCTTTCATCccagc                                                                                                                                >  4:463885/1‑149 (MQ=255)
                    aataaCGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCg                                                                                                                            >  1:500792/1‑149 (MQ=255)
                     ataaCGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGc                                                                                                                           <  1:685893/149‑1 (MQ=255)
                          gTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCctgc                                                                                                                      >  4:554754/1‑149 (MQ=255)
                          gTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCctgc                                                                                                                      >  4:20030/1‑149 (MQ=255)
                          gTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCctgc                                                                                                                      >  3:664914/1‑149 (MQ=255)
                          gTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCctgc                                                                                                                      >  1:796476/1‑149 (MQ=255)
                               gcggcggGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGAt                                                                                                                 <  1:71494/149‑1 (MQ=255)
                               gcggcggGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGAt                                                                                                                 <  4:272295/149‑1 (MQ=255)
                                  gcggGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCgc                                                                                                              <  4:214264/149‑1 (MQ=255)
                                     ggATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCgccag                                                                                                           <  2:261231/149‑1 (MQ=255)
                                     ggATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCgccag                                                                                                           >  2:177338/1‑149 (MQ=255)
                                      gATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATc                                                                                                                >  3:471620/1‑143 (MQ=255)
                                      gATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATc                                                                                                                <  4:471620/143‑1 (MQ=255)
                                       aTAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGattat                                                                                                                                                >  4:198622/1‑110 (MQ=255)
                                       aTAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGattat                                                                                                                                                <  3:198622/110‑1 (MQ=255)
                                       aTAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCgccagcc                                                                                                         >  1:302012/1‑149 (MQ=255)
                                          aTCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGt                                                                                                      >  2:811698/1‑149 (MQ=255)
                                           tCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTa                                                                                                     <  3:61632/149‑1 (MQ=255)
                                           tCTTTTCGCGCCATCGCCACATCACCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTAGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTa                                                                                                     >  2:597012/1‑149 (MQ=255)
                                             ttttCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAgg                                                                                                   <  1:313173/149‑2 (MQ=255)
                                              tttCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCg                                                                                                  <  2:796476/149‑1 (MQ=255)
                                                tCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGt                                                                                                <  1:890069/149‑1 (MQ=255)
                                                   cgcCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTa                                                                                             <  1:117138/149‑1 (MQ=255)
                                                        tCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTc                                                                                        <  2:488326/149‑1 (MQ=255)
                                                        tCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTc                                                                                        <  2:641594/149‑1 (MQ=255)
                                                        tCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTc                                                                                        <  1:224137/149‑1 (MQ=255)
                                                         cGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTca                                                                                       <  2:111468/149‑1 (MQ=255)
                                                         cGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTca                                                                                       <  3:479975/149‑1 (MQ=255)
                                                         cGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTCCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTca                                                                                       <  1:41938/149‑1 (MQ=255)
                                                             acaTCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTcagcag                                                                                   >  2:515692/1‑149 (MQ=255)
                                                                     cagccccagcACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTcc                                                                           >  1:739548/1‑149 (MQ=255)
                                                                                 tGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCgttgt                                                               <  3:554754/149‑1 (MQ=255)
                                                                                   cATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGttgttg                                                             <  2:159982/149‑1 (MQ=255)
                                                                                    aTAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTtgc                                                            >  1:397425/1‑149 (MQ=255)
                                                                                       aCTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTtgctgc                                                         <  4:791794/149‑1 (MQ=255)
                                                                                           tCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCtgtg                                                     >  1:122784/1‑149 (MQ=255)
                                                                                            cGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGc                                                    <  1:223948/149‑1 (MQ=255)
                                                                                             ggTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCg                                                   >  4:407393/1‑149 (MQ=255)
                                                                                              gtgtTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCgg                                                  >  1:715830/1‑149 (MQ=255)
                                                                                               tgtTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCgt                                                                  <  4:28186/132‑1 (MQ=255)
                                                                                               tgtTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCgt                                                                  >  3:28186/1‑132 (MQ=255)
                                                                                                    tCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAAGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAg                                            >  3:800048/1‑149 (MQ=255)
                                                                                                    tCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAAGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAg                                            >  1:102240/1‑149 (MQ=255)
                                                                                                     cGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGt                                           <  3:20030/149‑1 (MQ=255)
                                                                                                      gACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTg                                          >  4:622213/1‑149 (MQ=255)
                                                                                                       aCATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGc                                         >  2:689140/1‑149 (MQ=255)
                                                                                                          ttACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑g                                     <  4:424436/149‑1 (MQ=255)
                                                                                                          ttACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑g                                     <  3:134236/149‑1 (MQ=255)
                                                                                                              gcgcCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑GGGCt                                 >  1:629480/1‑149 (MQ=255)
                                                                                                              gcgcCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑GGGCt                                 <  3:239904/149‑1 (MQ=255)
                                                                                                                  ctgctgGAACAGGCGTTCTGCCCGGTCAAGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑GGGCttttt                             <  4:800048/149‑1 (MQ=255)
                                                                                                                  ctgctgGAACAGGCGTTCTGCCCGGTCAAGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑GGGCttttt                             <  2:102240/149‑1 (MQ=255)
                                                                                                                   tgctgGAACAGGCGTTCTGCCCGGTCAGGACTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCCTTTGTCGTTGTTGCTGCTGTGCGGGTCCAGAGCGACGGGCttttt                             >  4:346558/1‑149 (MQ=255)
                                                                                                                     ctgGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑GGGCTTTTTCCa                          >  2:222131/1‑149 (MQ=255)
                                                                                                                     ctgGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑GGGCTTTTTCCa                          <  1:831758/149‑1 (MQ=255)
                                                                                                                       ggAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑GGGCTTTTTCCAGa                        >  2:663095/1‑149 (MQ=255)
                                                                                                                            aGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGTGTCCAGTGCGA‑GGGCTTTTTCCAGATTCGc                   >  1:271158/1‑149 (MQ=255)
                                                                                                                            aGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑GGGCTTTTTCCAGATTCGc                   <  3:646160/149‑1 (MQ=255)
                                                                                                                            aGGCGTTCTGCACGCTCAGGATTATTGGCTTTCAGCGCAGCAATGCCCTGCTGGATCGCCAGCCAGTAGCGGGTTACTTTCAGCATACTGTTCCAGTTGGCGTTGTTGCAGCTGTGCAGGGCCAGTGCGA‑TGGCTGTTTGCAGAttggt                   >  4:108146/1‑146 (MQ=255)
                                                                                                                                   cTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑GGGCTTTTTCCAGATTCGCCACTGCa            >  3:205380/1‑149 (MQ=255)
                                                                                                                                            cAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑GGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGa   <  1:222959/149‑1 (MQ=255)
                                                                                                                                             aGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑GGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGAt  >  4:29171/1‑149 (MQ=255)
                                                                                                                                             aGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑GGGCTTTCTCCAGATTCGCCACTGCATTGGCGCGAt  <  4:704379/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CATACGCAAGGTCTGCTGATAATAACGTTCGGCGGCGGGATAATCTTTTCGCGCCATCGCCACATCGCCCAGCCCCAGCACTGCATAACTGTCGGTGTTATCGACATTACGCGCCTGCTGGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGA‑GGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGAT  >  NC_000913/3687932‑3688221

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