Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,235,687 A→T 5.4% intergenic (+30/+22) fadR → / ← ycgB fatty acid metabolism regulon transcriptional regulator/SpoVR family stationary phase protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,235,6870AT 5.4% 265.8 / 3.4 75intergenic (+30/+22)fadR/ycgBfatty acid metabolism regulon transcriptional regulator/SpoVR family stationary phase protein
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base T (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCG  >  NC_000913/1235542‑1235831
                                                                                                                                                 |                                                                                                                                                
aaGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACg                                                                                                                                               >  4:491948/1‑149 (MQ=255)
aaGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACg                                                                                                                                               >  1:15508/1‑149 (MQ=255)
  ggCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAg                                                                                                                                             >  1:175736/1‑149 (MQ=255)
   gcgcgcACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAgg                                                                                                                                            <  4:1280539/149‑1 (MQ=255)
   gcgcgcACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAgg                                                                                                                                            <  3:491948/149‑1 (MQ=255)
    cgcgcACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAggg                                                                                                                                           <  2:175736/149‑1 (MQ=255)
     gcgcACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAgggg                                                                                                                                          <  3:698682/149‑1 (MQ=255)
         aCGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAg                                                                                                                                             <  1:736916/142‑1 (MQ=255)
         aCGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAg                                                                                                                                             >  2:736916/1‑142 (MQ=255)
          cGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGttttt                                                                                                                                     >  4:444383/1‑149 (MQ=255)
             tCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAAc                                                                                                                                                <  4:522633/135‑1 (MQ=255)
             tCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAAc                                                                                                                                                >  3:522633/1‑135 (MQ=255)
                  tgtACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGttt                                                                                                                             >  1:1151807/1‑149 (MQ=255)
                   gtACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGtttt                                                                                                                            <  2:1151807/149‑1 (MQ=255)
                    tACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTttg                                                                                                                                   <  3:404666/141‑1 (MQ=255)
                    tACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTttg                                                                                                                                   >  4:404666/1‑141 (MQ=255)
                    tACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGttttt                                                                                                                           >  3:603306/1‑149 (MQ=255)
                       gAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACa                                                                                                                        <  3:310955/149‑1 (MQ=255)
                         aaCAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGa                                                                                                                      >  2:387095/1‑149 (MQ=255)
                              tGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTcc                                                                                                                 <  4:603306/149‑1 (MQ=255)
                              tGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTcc                                                                                                                 <  1:387095/149‑1 (MQ=255)
                                cGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCa                                                                                                               <  3:444383/149‑1 (MQ=255)
                                  tCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGtttttt                                                                                                                                    >  3:167539/1‑126 (MQ=255)
                                  tCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGtttttt                                                                                                                                    <  4:167539/126‑1 (MQ=255)
                                           gCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTTATTTAGCCATTCAGGGGCGATAAACCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTttgt                                                                                                                                  <  3:596065/119‑1 (MQ=255)
                                           gCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTttgt                                                                                                                                  >  4:596065/1‑119 (MQ=255)
                                           gCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCggg                                                                                                    >  1:965209/1‑149 (MQ=255)
                                                        gagaTTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGttttt                                                                                                                                     >  1:435562/1‑103 (MQ=255)
                                                        gagaTTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGttttt                                                                                                                                     <  2:435562/103‑1 (MQ=255)
                                                        gagaTTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCa                                                                                       >  3:459145/1‑149 (MQ=255)
                                                        gagaTTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTCTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCa                                                                                       >  2:805712/1‑149 (MQ=255)
                                                                   ccGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTCTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGAtgctg                                                                            <  1:805712/149‑1 (MQ=255)
                                                                           aGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTc                                                                    >  1:395904/1‑149 (MQ=255)
                                                                            gAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCGTTTGAGGGGtt                                                                                                                                        >  1:837824/1‑80 (MQ=255)
                                                                            gAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCGTTTGAGGGGtt                                                                                                                                        <  2:837824/80‑1 (MQ=255)
                                                                             aaaaaTCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGt                                                                  >  1:457199/1‑149 (MQ=255)
                                                                                 aTCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTg                                                              <  4:459145/149‑1 (MQ=255)
                                                                                         ggTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGa                                                      <  2:965209/149‑1 (MQ=255)
                                                                                                  gCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGt                                                          <  3:587549/136‑1 (MQ=255)
                                                                                                  gCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGt                                                          >  4:587549/1‑136 (MQ=255)
                                                                                                     aTTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCaa                                          >  1:1233071/1‑149 (MQ=255)
                                                                                                       tCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAAt                                        <  2:395904/149‑1 (MQ=255)
                                                                                                         aGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTttgt                                                                                                                                  <  4:455932/57‑1 (MQ=255)
                                                                                                         aGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTttgt                                                                                                                                  >  3:455932/1‑57 (MQ=255)
                                                                                                                  aaTCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCg                             >  2:581279/1‑149 (MQ=255)
                                                                                                                      ccTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGc                         <  2:1233071/149‑1 (MQ=255)
                                                                                                                        ttCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCaaa                                         <  2:1420121/131‑1 (MQ=255)
                                                                                                                        ttCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCaaa                                         >  1:1420121/1‑131 (MQ=255)
                                                                                                                        ttCCGTCTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGcac                       >  2:1040319/1‑149 (MQ=255)
                                                                                                                         tCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGcaca                      <  2:457199/149‑1 (MQ=255)
                                                                                                                            gTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATccc                                     >  1:907221/1‑131 (MQ=255)
                                                                                                                            gTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATccc                                     <  2:907221/131‑1 (MQ=255)
                                                                                                                             tttAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCt                  <  1:581279/149‑1 (MQ=255)
                                                                                                                                  agagCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACAt                     <  1:508623/141‑1 (MQ=255)
                                                                                                                                  agagCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACAt                     >  2:508623/1‑141 (MQ=255)
                                                                                                                                   gagCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGa            >  3:751730/1‑149 (MQ=255)
                                                                                                                                      aaaaCCCCTCGTTTGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAg              >  2:513391/2‑144 (MQ=255)
                                                                                                                                       aaaCCCCTCGTTTGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGc                                                             <  1:926805/96‑1 (MQ=255)
                                                                                                                                       aaaCCCCTCGGTTGAGGGTATTTGTGGTGTTTTTAGAGATTTACCATGCGTGACGGGAGACGATCCAGCAACTCGATGCTTCTGTCTTCTTTTTGCTGTTCGAGCATCACATCTAATTCCCCCAGCCTGTCTACACGCTTCAGGAcctc        >  4:264877/1‑149 (MQ=255)
                                                                                                                                       aaaCCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCca                                   <  1:1380806/122‑1 (MQ=255)
                                                                                                                                       aaaCCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCca                                   >  2:1380806/1‑122 (MQ=255)
                                                                                                                                        aaCCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGActtct       >  1:544597/1‑149 (MQ=255)
                                                                                                                                          ccccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCttt     <  2:15508/149‑1 (MQ=255)
                                                                                                                                           cccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGa                            >  2:1190048/1‑125 (MQ=255)
                                                                                                                                           cccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTg    >  1:968472/1‑149 (MQ=255)
                                                                                                                                           cccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTg    >  4:206608/1‑149 (MQ=255)
                                                                                                                                           cccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTg    >  4:355959/1‑149 (MQ=255)
                                                                                                                                           cccTCAAACGAGGGGTTTCTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCCACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGa                            <  1:1190048/125‑1 (MQ=255)
                                                                                                                                            ccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTgc   <  1:1040319/149‑1 (MQ=255)
                                                                                                                                            ccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTgc   <  3:206608/149‑1 (MQ=255)
                                                                                                                                             cTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGc                         >  2:1331513/1‑126 (MQ=255)
                                                                                                                                             cTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGc                         <  1:1331513/126‑1 (MQ=255)
                                                                                                                                             cTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTgcg  >  3:1221431/1‑149 (MQ=255)
                                                                                                                                              tCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGAtg                                                                               <  2:1416597/71‑1 (MQ=255)
                                                                                                                                              tCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGAtg                                                                               >  1:1416597/1‑71 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCG  >  NC_000913/1235542‑1235831

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