Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,698,520 A→C T103P (ACC→CCC)  bcsG → DUF3260 family cellulose production inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,698,5200AC96.6% 492.6 / 2.4 151T103P (ACC→CCC) bcsGDUF3260 family cellulose production inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base C (76/68);  total (79/72)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 152 total aligned reads displayed.

TGCCCTGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTACCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTATTCCTGTCACAAT  >  NC_000913/3698423‑3698616
                                                                                                 |                                                                                                
tGCCCTGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTccc                                                                                                <  1:739535/100‑1 (MQ=255)
 gCCCTGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCg                                                                                               >  1:4034320/1‑100 (MQ=255)
 gCCCTGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGCAAACCATAATGAGCCAGGGTTCGCAGGTGACGGGGTTCAGTCCCg                                                                                               >  2:3811619/1‑100 (MQ=255)
  cccTGCCGATCGGCTTTGCTTCGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGa                                                                                              >  1:3483414/1‑100 (MQ=255)
   ccTGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGat                                                                                             <  2:1808938/100‑1 (MQ=255)
    cTGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGatt                                                                                            >  2:1390066/1‑100 (MQ=255)
    cTGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGatt                                                                                            <  1:3527589/100‑1 (MQ=255)
    cTGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGatt                                                                                            >  1:1511837/1‑100 (MQ=255)
     tGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTACCGatta                                                                                           >  1:455122/1‑100 (MQ=255)
      gCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGattat                                                                                          >  2:1500902/1‑100 (MQ=255)
      gCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGattat                                                                                          <  1:120751/100‑1 (MQ=255)
          aTCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGTATAATGAGCCAGGGTTCGCAGGTGGCGGGCTTCAGTCCCGATTATTTaa                                                                                      >  2:1191046/1‑100 (MQ=255)
          aTCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTaa                                                                                      <  1:2607730/100‑1 (MQ=255)
            cGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATc                                                                                    <  2:2854090/100‑1 (MQ=255)
            cGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTACCGATTATTTAATc                                                                                    >  1:1708487/1‑100 (MQ=255)
             ggCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCg                                                                                   <  2:207916/100‑1 (MQ=255)
             ggCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCg                                                                                   >  1:257402/1‑100 (MQ=255)
             ggCTTTGCTTTGTTCTGCCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCACAGGGGGCGGGGTGCAGTCCCGATTATTTaagct                                                                                   >  2:233801/1‑97 (MQ=255)
             ggCTTTGCGTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTTGCGGGGTTCAGTCCCAAATATTTAATCg                                                                                   <  2:651424/100‑1 (MQ=255)
                tttgctttgTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGAcc                                                                                <  1:1677149/100‑1 (MQ=255)
                   gctttgTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTGCCGATTATTTAATCGACCTTg                                                                             >  2:791714/1‑100 (MQ=255)
                   gctttgTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTg                                                                             <  1:875997/100‑1 (MQ=255)
                    ctttgTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGt                                                                            >  1:3637764/1‑100 (MQ=255)
                    ctttgTTCTGGCATGACACCTGGTTGCCTGGCCCAGAAAGCATAATGAGCCAGAGTTCGCAGGTGGCGGGGTTCGGTCCCGATTATTTAATCGACCTTGt                                                                            >  1:3792979/1‑100 (MQ=255)
                     tttgTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTc                                                                           <  1:2482849/100‑1 (MQ=255)
                     tttgTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTc                                                                           <  2:1710869/100‑1 (MQ=255)
                     tttgTTCTGGCATGACACCTGGTTGCCTGCCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTCATCGACCTTGTc                                                                           <  2:1728831/100‑1 (MQ=255)
                      ttgTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTca                                                                          <  1:1921063/100‑1 (MQ=255)
                       tgTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCGAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAGTCGACCTTGTcac                                                                         <  2:3601771/100‑1 (MQ=255)
                       tgTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTcac                                                                         <  1:2730723/100‑1 (MQ=255)
                       tgTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTcac                                                                         <  2:1038353/100‑1 (MQ=255)
                       tgTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTcac                                                                         <  1:183736/100‑1 (MQ=255)
                        gTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTcaca                                                                        >  1:2324498/1‑100 (MQ=255)
                          tCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACg                                                                      <  1:228193/100‑1 (MQ=255)
                            tGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCt                                                                    <  1:732836/100‑1 (MQ=255)
                             gacATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGATCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCtt                                                                   >  2:3706622/3‑100 (MQ=255)
                              gCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCttt                                                                  >  2:1080938/1‑100 (MQ=255)
                                atctcacCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTttat                                                                <  1:1080034/96‑1 (MQ=255)
                                aTGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTttat                                                                >  1:3181048/1‑100 (MQ=255)
                                aTGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTttat                                                                <  1:3528742/100‑1 (MQ=255)
                                aTGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTttat                                                                >  1:2860982/1‑100 (MQ=255)
                                aTGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTttat                                                                <  1:1475842/100‑1 (MQ=255)
                                 tGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTttatt                                                               <  1:3745909/100‑1 (MQ=255)
                                 tGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTttatt                                                               <  2:3143389/100‑1 (MQ=255)
                                 tGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTACCGATTATTTAATCGACCTTGTCACACGCTttatt                                                               <  2:446480/100‑1 (MQ=255)
                                    cacCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTttatagat                                                            >  2:1058078/1‑96 (MQ=255)
                                    cacCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAAc                                                            <  2:1158823/100‑1 (MQ=255)
                                     acCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACt                                                           >  1:2869988/1‑100 (MQ=255)
                                      ccTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTACCGATTATATAATCGACCTTGTCACACGCTTTATTAACTg                                                          <  1:3269080/100‑1 (MQ=255)
                                        tGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGc                                                        <  2:2826788/100‑1 (MQ=255)
                                        tGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGc                                                        <  1:1041803/100‑1 (MQ=255)
                                          gTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAg                                                      >  1:264296/1‑100 (MQ=255)
                                             gCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAgatg                                                   >  1:875445/1‑100 (MQ=255)
                                             gCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAgatg                                                   >  1:1007894/1‑100 (MQ=255)
                                              ccTGGCCCGGAAAGCATAATGCGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAgatga                                                  <  2:1237102/100‑1 (MQ=255)
                                              ccTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAgatga                                                  >  2:2813187/1‑100 (MQ=255)
                                                tGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGAtt                                                <  1:2499240/100‑1 (MQ=255)
                                                tGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGAtt                                                <  1:185946/100‑1 (MQ=255)
                                                   cccGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGCCACACGCTTTATTAACTGGCAGATGATTggg                                             <  2:1649976/100‑1 (MQ=255)
                                                    ccGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTgggg                                            >  1:1588438/1‑100 (MQ=255)
                                                    ccGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATGgggg                                            >  2:687111/1‑100 (MQ=255)
                                                    ccGGAAAGCAGAATGGGCCGGGGTTCGCAGGTGGCGGGGTTCAGTCGCGATTATTTAATCGACCTTGTCACACGCTTTATTCGCTGGCAGATGATTgggg                                            <  2:2466110/100‑1 (MQ=255)
                                                      ggAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGcc                                          >  1:3316153/1‑100 (MQ=255)
                                                         aaGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCAtt                                       >  1:45096/1‑100 (MQ=255)
                                                           gCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCAtttt                                     >  1:3917241/1‑100 (MQ=255)
                                                              tAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGt                                  >  2:1653992/1‑100 (MQ=255)
                                                                aTGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGttt                                <  1:886834/100‑1 (MQ=255)
                                                                aTGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATGGGGGCCATTTTTGttt                                >  2:603819/1‑100 (MQ=255)
                                                                 tGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGtttt                               >  2:2936431/1‑100 (MQ=255)
                                                                 tGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGtttt                               >  1:2605245/1‑100 (MQ=255)
                                                                  gAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTtta                              <  2:3703638/100‑1 (MQ=255)
                                                                   aGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGGTGAGTGGGGCCATTTTTGTTttat                             >  2:390742/1‑100 (MQ=255)
                                                                   aGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTttat                             <  1:575440/100‑1 (MQ=255)
                                                                   aGCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTttat                             <  1:3071859/100‑1 (MQ=255)
                                                                    gCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTttatt                            >  1:2023232/1‑100 (MQ=255)
                                                                    gCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTttatt                            <  2:2665135/100‑1 (MQ=255)
                                                                    gCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTttatt                            >  1:3914984/1‑100 (MQ=255)
                                                                    gCCAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTttatt                            >  1:3757115/1‑100 (MQ=255)
                                                                     ccAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGCGGCCATttttttttttttt                           >  2:3572936/1‑99 (MQ=255)
                                                                      cAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAg                          >  2:2701419/1‑100 (MQ=255)
                                                                      cAGGGTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAg                          <  2:1800281/100‑1 (MQ=255)
                                                                          gTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGc                      <  1:3945145/100‑1 (MQ=255)
                                                                          gTTCGCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGc                      >  2:3239006/1‑100 (MQ=255)
                                                                              gCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTgg                  <  2:222373/100‑1 (MQ=255)
                                                                              gCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTgg                  >  1:2930578/1‑100 (MQ=255)
                                                                              gCAGGTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTgg                  >  2:3859193/1‑100 (MQ=255)
                                                                               cAGGTGGCGGGGTGCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGt                 <  2:3132633/100‑1 (MQ=255)
                                                                                 ggTGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTa               >  2:1530301/1‑100 (MQ=255)
                                                                                   tGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTAtt             >  1:2895851/1‑100 (MQ=255)
                                                                                   tGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTAtt             <  1:2170255/100‑1 (MQ=255)
                                                                                   tGGCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTAtt             >  1:1787270/1‑100 (MQ=255)
                                                                                    ggCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTATTc            >  2:2332034/1‑100 (MQ=255)
                                                                                     gCGGGGTTGAGCCACGATTATTCGATCGACGGTGTCGCTCGTTTTTTTAGCTGGCTTATGAGTGGGGCCATTTTTGTTTTATTAGTGGCGTGGTTATTcc           <  2:3218842/100‑1 (MQ=255)
                                                                                     gCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTATTcc           <  1:3526611/100‑1 (MQ=255)
                                                                                     gCGGGGTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTATTcc           <  1:3828873/100‑1 (MQ=255)
                                                                                        gggTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTATTCCTGt        >  2:3808380/1‑100 (MQ=255)
                                                                                          gTTCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTGATTAGTGGCCTGGTTATTCCTGTca      >  1:3410012/1‑100 (MQ=255)
                                                                                           ttCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTTTTTTTTTTTTCGCCTGGATATCCCTGTcac     >  2:905791/1‑100 (MQ=255)
                                                                                           ttCAGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTATTCCTGTcac     <  2:3951035/100‑1 (MQ=255)
                                                                                              aGTCCCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTGATTAGTGGCCTGGTTATTCCTGTCACAAt  >  1:2210003/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TGCCCTGCCGATCGGCTTTGCTTTGTTCTGGCATGACACCTGGTTGCCTGGCCCGGAAAGCATAATGAGCCAGGGTTCGCAGGTGGCGGGGTTCAGTACCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTATTCCTGTCACAAT  >  NC_000913/3698423‑3698616

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