Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,351,228 G→C 7.4% intergenic (‑189/+179) fabI ← / ← ycjD enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,351,2280GC 7.4% 350.0 / 13.2 106intergenic (‑189/+179)fabI/ycjDenoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein
Reads supporting (aligned to +/- strand):  ref base G (46/52);  new base C (3/5);  total (49/57)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 111 total aligned reads displayed.

AGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAG  >  NC_000913/1351082‑1351376
                                                                                                                                                  |                                                                                                                                                    
aGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCagag                                                                                                                                                    <  2:658665/149‑1 (MQ=255)
  gCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAggg                                                                                                                                                  >  4:656733/1‑149 (MQ=255)
   cGATTCGTCGTTTCAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgggg                                                                                                                                                 <  3:67856/149‑1 (MQ=255)
     aTTCGTCGTTTTAGTAAACATTACGAACAGATAAACAGTTAGTATAATCCACCTGGCTGTGCGTAGCTATAGTTGTCAGGTCCGGCCGGAGCATGCTGCGGCAGGGCGGGCGCTTTTCCACTCACCCTGACGCTCTTCCCACTgggggt                                                                                                                                               >  2:428260/1‑147 (MQ=255)
        cgtcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAgg                                                                                                                                            <  2:293119/149‑1 (MQ=255)
        cgtcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAgg                                                                                                                                            <  1:867501/149‑1 (MQ=255)
         gtcgtTTTAGTAATCAGTACGAACAGATAAACGGTTATTCTAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGTGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAggg                                                                                                                                           <  1:621259/149‑1 (MQ=255)
         gtcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAggg                                                                                                                                           >  4:380544/1‑149 (MQ=255)
         gtcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAggg                                                                                                                                           <  3:128229/149‑1 (MQ=255)
         gtcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAggg                                                                                                                                           <  2:666944/149‑1 (MQ=255)
            gtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGAc                                                                                                                                        >  1:696465/1‑149 (MQ=255)
                  gTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAtt                                                                                                                                  >  1:209225/1‑149 (MQ=255)
                  gTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAtt                                                                                                                                  >  3:90822/1‑149 (MQ=255)
                  gTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAtt                                                                                                                                  <  2:474322/149‑1 (MQ=255)
                  gTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAtt                                                                                                                                  <  3:380544/149‑1 (MQ=255)
                     aaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATtgtg                                                                                                                               >  1:594744/1‑149 (MQ=255)
                     aaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATtgtg                                                                                                                               <  2:335941/149‑1 (MQ=255)
                     aaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATtgtg                                                                                                                               <  3:122332/149‑1 (MQ=255)
                      aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGc                                                                                                                              >  3:244035/1‑149 (MQ=255)
                             gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTGTTGTGCAAATAtt                                                                                                                       >  3:108897/1‑149 (MQ=255)
                             gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATAtt                                                                                                                       <  2:331104/149‑1 (MQ=255)
                             gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATAtt                                                                                                                       <  2:140506/149‑1 (MQ=255)
                                cAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGtt                                                                                                                    <  4:488730/149‑1 (MQ=255)
                                  gATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTAc                                                                                                                  <  3:162182/149‑1 (MQ=255)
                                              ttatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGCGCAGGCTGCAGCAGGGGGGGCGCTTTTCCCATCTCCCTCGCCCTGTCCCCAGAGGGGCGAGGGGACCGTACTGTGCAAATATTGTTACCCCCGGaacaaa                                                                                                      <  2:421548/149‑1 (MQ=255)
                                               tataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaac                                                                                                     >  3:381773/1‑149 (MQ=255)
                                                 taATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAg                                                                                                   <  1:24374/149‑1 (MQ=255)
                                                 taATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTATTGTGCGAATGTTCTTATCCCAGCCACAAACAg                                                                                                   >  1:692953/1‑149 (MQ=38)
                                                  aaTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCACATATTGTTACCCCAGCAACAAACAgg                                                                                                  >  2:252555/1‑149 (MQ=255)
                                                  aaTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAgg                                                                                                  >  4:399669/1‑149 (MQ=255)
                                                   aTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGc                                                                                                 <  1:377388/149‑1 (MQ=255)
                                                   aTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGc                                                                                                 >  1:654193/1‑149 (MQ=255)
                                                    tCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGCGGCGAGGGGACCGTATTGTGCAAATCTTGTTACCCCAGGCACAAACAGGCt                                                                                                >  3:807130/1‑149 (MQ=255)
                                                              tgtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGccc                                                                                      >  3:608135/1‑149 (MQ=255)
                                                                   gTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAAc                                                                                 <  2:231156/149‑1 (MQ=255)
                                                                    tAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                                                                                                                             <  3:557482/88‑1 (MQ=255)
                                                                      gCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAggg                                                                                                                                           <  1:763218/88‑1 (MQ=38)
                                                                      gCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAggg                                                                                                                                           >  2:763218/1‑88 (MQ=38)
                                                                       cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaa                                                                                                      >  4:591367/1‑124 (MQ=255)
                                                                       cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaa                                                                                                      <  3:591367/124‑1 (MQ=255)
                                                                       cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCtt                                                                             >  3:171147/1‑149 (MQ=255)
                                                                        tataGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCttt                                                                            >  2:886330/1‑149 (MQ=255)
                                                                        tataGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCttt                                                                            <  2:209225/149‑1 (MQ=255)
                                                                            gTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATg                                                                        <  1:42022/149‑1 (MQ=255)
                                                                              tGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGc                                                                      <  4:38824/149‑1 (MQ=255)
                                                                              tGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGc                                                                      <  4:381773/149‑1 (MQ=255)
                                                                              tGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGc                                                                      <  3:656733/149‑1 (MQ=255)
                                                                               gCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCg                                                                     >  3:159317/1‑149 (MQ=255)
                                                                                ccAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGa                                                                    <  2:183629/149‑1 (MQ=255)
                                                                                ccAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGa                                                                    >  4:434180/1‑149 (MQ=255)
                                                                                  aGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTTTTACCCCAGCAACAAACAg                                                                                                   >  4:176620/1‑116 (MQ=255)
                                                                                  aGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAg                                                                                                   <  3:176620/116‑1 (MQ=255)
                                                                                   ggTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATgg                                                                 <  4:837760/149‑1 (MQ=255)
                                                                                    gTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGAAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGc                                                                >  4:581178/1‑149 (MQ=255)
                                                                                        gACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTggg                                                            >  4:420595/1‑149 (MQ=255)
                                                                                        gACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTggg                                                            >  2:362024/1‑149 (MQ=255)
                                                                                            ggAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTgggttag                                                        <  1:656014/149‑5 (MQ=255)
                                                                                                cAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCa                                                    <  3:505989/149‑1 (MQ=255)
                                                                                                  ggCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGa                                                  <  1:69106/149‑1 (MQ=255)
                                                                                                  ggCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAACCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGa                                                  >  3:365474/1‑149 (MQ=255)
                                                                                                   gCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGAACGGTTCAGAc                                                 <  3:581178/149‑1 (MQ=255)
                                                                                                      gCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCtt                                              <  3:228836/149‑1 (MQ=255)
                                                                                                      gCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCtt                                              <  3:399669/149‑1 (MQ=255)
                                                                                                                   cgcTTTTCCCCTCACCCTAACCCTCTCCCCAGAggg                                                                                                                                                  >  4:626329/1‑36 (MQ=37)
                                                                                                                   cgcTTTTCCCCTCACCCTAACCCTCTCCCCAGAggg                                                                                                                                                  <  3:626329/36‑1 (MQ=37)
                                                                                                                   cgcTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATtctc                                 <  1:154724/149‑1 (MQ=255)
                                                                                                                     cTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCa                               <  3:509966/149‑1 (MQ=255)
                                                                                                                      ttttCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAg                              >  3:285603/1‑149 (MQ=255)
                                                                                                                          ccccTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCAcc                          <  3:106519/149‑1 (MQ=255)
                                                                                                                             cTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtc                       <  4:618171/149‑1 (MQ=255)
                                                                                                                               cACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtctc                     <  1:711961/149‑1 (MQ=255)
                                                                                                                                  ccTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATtct                                  <  4:762633/133‑1 (MQ=255)
                                                                                                                                  ccTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATtct                                  >  3:762633/1‑133 (MQ=255)
                                                                                                                                  ccTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCAt                  >  1:480880/1‑149 (MQ=255)
                                                                                                                                      aCCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGttt              >  3:389889/1‑149 (MQ=255)
                                                                                                                                       cccTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGtttt             <  1:239859/149‑1 (MQ=255)
                                                                                                                                        ccTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATg                                                                        >  2:401442/1‑89 (MQ=255)
                                                                                                                                        ccTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATg                                                                        <  1:401442/89‑1 (MQ=255)
                                                                                                                                         ctctCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTca           <  2:594744/149‑1 (MQ=255)
                                                                                                                                         ctctCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTca           <  3:138559/149‑1 (MQ=255)
                                                                                                                                         ctctCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTca           <  1:47168/149‑1 (MQ=255)
                                                                                                                                         cgctcCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCTCCTAACCCTTTCATGGCGATGGCTGCGACGGTTCAGACCTTGCCGAATATTCTCCTGCACCGCCTCCATGGTTTca           <  2:124215/147‑1 (MQ=255)
                                                                                                                                           ctcCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcacc         >  2:146553/1‑149 (MQ=255)
                                                                                                                                              cccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAata                                      >  1:616407/1‑117 (MQ=255)
                                                                                                                                              cccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAata                                      <  2:616407/117‑1 (MQ=255)
                                                                                                                                              cccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTca           <  3:616747/144‑1 (MQ=255)
                                                                                                                                              cccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTca           >  4:616747/1‑144 (MQ=255)
                                                                                                                                              cccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACcaca      >  4:20422/1‑149 (MQ=255)
                                                                                                                                              cccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACcaca      >  2:619683/1‑149 (MQ=255)
                                                                                                                                              cccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACcaca      >  4:731244/1‑149 (MQ=255)
                                                                                                                                               ccAGAGGTGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaa     <  3:567201/149‑1 (MQ=255)
                                                                                                                                               ccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCttt                                                                            <  1:640586/78‑1 (MQ=255)
                                                                                                                                               ccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCttt                                                                            >  2:640586/1‑78 (MQ=255)
                                                                                                                                               ccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcac          <  4:333093/144‑1 (MQ=255)
                                                                                                                                               ccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcac          >  3:333093/1‑144 (MQ=255)
                                                                                                                                               ccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaa     >  4:339681/1‑149 (MQ=255)
                                                                                                                                                 agagGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaaca   <  2:654193/149‑1 (MQ=255)
                                                                                                                                                 agagGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaaca   >  1:254806/1‑149 (MQ=255)
                                                                                                                                                 agagGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaaca   >  3:512152/1‑149 (MQ=255)
                                                                                                                                                  gagGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAg  <  3:648778/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAG  >  NC_000913/1351082‑1351376

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