Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,460,614 C→T 8.7% intergenic (‑145/+36) fadI ← / ← yfcZ beta‑ketoacyl‑CoA thiolase, anaerobic, subunit/UPF0381 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,460,6140CT 8.7% 224.2 / 9.0 81intergenic (‑145/+36)fadI/yfcZbeta‑ketoacyl‑CoA thiolase, anaerobic, subunit/UPF0381 family protein
Reads supporting (aligned to +/- strand):  ref base C (43/30);  new base T (4/3);  total (47/34)
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.28e-01

CCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGG  >  NC_000913/2460466‑2460763
                                                                                                                                                    |                                                                                                                                                     
ccATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGcc                                                                                                                                                      <  2:94079/150‑1 (MQ=255)
ccATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGACAAGAGGAGAAAAGGGAAAGTGGGAGATATGACACAGAGAGAAAAGAAGATGAGAAAAGAAAACGCCCCTGcc                                                                                                                                                      >  2:788109/1‑150 (MQ=255)
ccATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTAAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGcc                                                                                                                                                      >  2:980392/1‑150 (MQ=255)
 cATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCa                                                                                                                                                     >  1:323072/1‑150 (MQ=255)
       acTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACt                                                                                                                                               >  1:1090677/1‑150 (MQ=255)
       acTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACt                                                                                                                                               <  2:910294/150‑1 (MQ=255)
         tCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGa                                                                                                                                             >  2:793078/1‑150 (MQ=255)
          ccccTGTAAAAATATAAATAAAGGGGTCTGACCTGATCATAGGCTTAACCAGTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGTCAGCAACTGAc                                                                                                                                            <  2:38800/150‑1 (MQ=255)
           cccTTTAAAAAAATAAAAAAAGTGGTCTGAACTCATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGACCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACa                                                                                                                                           <  1:619701/150‑1 (MQ=255)
             cTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAgg                                                                                                                                         >  1:928917/1‑150 (MQ=255)
              tGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAggg                                                                                                                                        >  2:557708/1‑150 (MQ=255)
               gTAAAAAAATAAATAAAGTGGTCTGACCGGAGCATAGTCTAAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAgggg                                                                                                                                       <  2:1048307/150‑1 (MQ=255)
                    aaaaTAAATAAAGTGGTATGACCTGATAATAGTCTGTACCATTTTTTTACATTTAGCAAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGATAAAAGAAGAAGCGAAAAGAAAACCCCCCTGCCAGCAACTGACAGGGGcgtcc                                                                                                                                  >  1:667385/1‑148 (MQ=255)
                     aaataaatAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgtacg                                                                                                                                 <  2:115799/150‑1 (MQ=255)
                     aaataaatAAAGTGGGCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCTAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgtacg                                                                                                                                 <  2:777232/150‑1 (MQ=255)
                      aataaataAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgtacgt                                                                                                                                >  1:520674/1‑150 (MQ=255)
                       ataaataaAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGTCAGTTGCTGGCAGGGGcgttttct                                                                                                                               <  2:1214423/150‑6 (MQ=255)
                       ataaataaAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGTCAGTTGCTGGCAGGGGcgttttct                                                                                                                               <  2:478285/150‑6 (MQ=255)
                       ataaataaAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgtacgta                                                                                                                               >  2:415473/1‑150 (MQ=255)
                        taaataaaGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAAGTGACAGGGGCgtacgtac                                                                                                                              >  1:647844/1‑150 (MQ=255)
                        taaataaaGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCgtacgtac                                                                                                                              <  2:221022/150‑1 (MQ=255)
                           ttacaGTGGTATGACCTGATCAAAGTCTTAACCATTTTTTTACATTAAGCCAATAGGAGAAAAGTGAAAGTGGGAGCTATGACACAGAGAGAAAACAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCg                                                                                                                           <  2:173615/146‑1 (MQ=255)
                              aaGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCgaga                                                                                                                        >  2:501234/1‑150 (MQ=255)
                              aaGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCgaga                                                                                                                        >  1:1251050/1‑150 (MQ=255)
                              aaGTGGTCAGACCTGATCCTCGGCTTAACCATTATTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTCGCAGCTATGACACAGATAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCgaga                                                                                                                        <  1:215531/150‑1 (MQ=255)
                               aGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTAAGCCAAGTGGAGAAAAGGGAAAGTGGGCGCTAGGACACAGAGAGAAAAGACGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGagag                                                                                                                       <  1:1055710/150‑1 (MQ=255)
                                   gTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAAt                                                                                                                   >  1:159145/1‑150 (MQ=255)
                                      tGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTaa                                                                                                                >  2:860847/1‑150 (MQ=255)
                                       gACCTGATCATAGTCTTAACCCTTTTTTTACATTTAGCCAAGTGCAGAAAAGGGCAAGTGGGAGCTATGACACAGAGAGAACAGTAGAAGAGAAAAGAACACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAAc                                                                                                               <  2:169311/150‑1 (MQ=255)
                                       gACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGTGAAAAGGGAAAGTGGGTGCTATGACACAGAGAGAAATGTAGAAGAGAAAAGAAAACGCCCCTTCCAGCCACTGAAAGGGGCGTAACTACCCGAGAGGAATTAAc                                                                                                               >  2:715132/1‑150 (MQ=255)
                                       gACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGCCAGGGGCGTACGTACCCGAGAGGAATTAAc                                                                                                               >  2:596549/1‑150 (MQ=255)
                                         ccTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGc                                                                                                             >  1:1178997/1‑150 (MQ=255)
                                           tGATCATAGTCTTAACCATTTTTTTACATTTAGACAAGTGGAGAAATGGGAAAGTGGGAGCTATTACAAAGAGAGAAAAGAAGAAGAGAAAAGAAAACTCCCCTGCAAGCAACTCACAGGGGCGTACTTAACCGAGAGGAATAACCGCac                                                                                                           >  1:206201/1‑149 (MQ=255)
                                                 tAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAg                                                                                                     <  1:18651/150‑1 (MQ=255)
                                                  aGTCTTACCCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGc                                                                                                    <  2:683677/150‑1 (MQ=255)
                                                   gTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCt                                                                                                   >  2:580217/1‑150 (MQ=255)
                                                     cTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTgg                                                                                                 <  2:557345/150‑1 (MQ=255)
                                                      ttAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATAAACGCAGACCCAGCTGGa                                                                                                <  2:515061/150‑1 (MQ=255)
                                                       tAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGATAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAATAGAAAACTCCCCTGCCAGCTACTGACAGGCCCGTCCGTACCCGACAGCAATAAACGCAGAACCAGATGTaa                                                                                               >  2:227506/1‑150 (MQ=255)
                                                          ccATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAAt                                                                                            >  1:199949/1‑150 (MQ=255)
                                                          ccATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAAt                                                                                            >  1:446135/1‑150 (MQ=255)
                                                          ccATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACACGAGAGGAATTAACGCAGACCCAGCTGGAAAAt                                                                                            >  1:592141/1‑150 (MQ=255)
                                                          ccATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCACTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCAGGAAAAt                                                                                            >  1:248431/1‑150 (MQ=255)
                                                                 tttACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCAtt                                                                                     >  1:1066473/1‑150 (MQ=255)
                                                                      aTTTAGCCAAGTGGATAAATGGGAATGTGGGAGCTTTTACACAGAGAGAAATGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGc                                                                                <  2:1217239/150‑1 (MQ=255)
                                                                           gCCAAGTGGAGAAAAGGGAAAGTGGAAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGc                                                                           >  1:446581/1‑150 (MQ=255)
                                                                           cccATGTGCATAAAAGGTAAAGTGGTTGCTATTACACAGAGATAAAAGAATAAGAGAAAATAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGc                                                                           <  1:1155135/149‑1 (MQ=255)
                                                                                    ataaaaGGTAAAGTGGGCGCTTTGACACAGACAGAAAAGAAGAAGAGAAAAGAACACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAg                                                                  <  1:1183418/148‑1 (MQ=255)
                                                                                    agaAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAg                                                                  <  2:323072/150‑1 (MQ=255)
                                                                                    agaAAAGGCAAACTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAg                                                                  <  1:112000/150‑1 (MQ=255)
                                                                                       aaaGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGa                                                               >  2:685245/1‑150 (MQ=255)
                                                                                               aGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTCGAAAATCAGCATTTCAGCTTCCCAGGCAAAAGAGAAATCGata                                                       >  2:1166751/1‑150 (MQ=255)
                                                                                                gTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGatat                                                      >  2:361141/1‑150 (MQ=255)
                                                                                                 tGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATc                                                     <  1:879278/150‑1 (MQ=255)
                                                                                                     aGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAACACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAgg                                                 <  1:1251944/150‑1 (MQ=255)
                                                                                                     aGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAgg                                                 <  2:135175/150‑1 (MQ=255)
                                                                                                      gCTATGACACTGAGAGAAAAGAAGAAGTCAAAACATCACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGgc                                                <  2:239145/150‑1 (MQ=255)
                                                                                                      gCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGgc                                                >  2:1090085/1‑150 (MQ=255)
                                                                                                       cTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGgcg                                               <  2:386362/150‑1 (MQ=255)
                                                                                                       cTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCAGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGgcg                                               >  1:34418/1‑150 (MQ=255)
                                                                                                             cacaGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACAcc                                         >  1:953708/1‑150 (MQ=255)
                                                                                                                    agaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGAc                                  <  1:501234/150‑1 (MQ=255)
                                                                                                                     gaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAAGCGATATCCAGGCGCACACCGGCAGAct                                 >  1:448432/1‑150 (MQ=255)
                                                                                                                         agaagaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTc                             <  2:1254067/150‑1 (MQ=255)
                                                                                                                         agaagaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTGCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGCCTCTTc                             >  1:761289/1‑150 (MQ=255)
                                                                                                                            agaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCCCACCGTCAGACTCTTCAGt                          >  1:265951/1‑150 (MQ=255)
                                                                                                                            agaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGt                          >  1:609060/1‑150 (MQ=255)
                                                                                                                                  gaaaagaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCTGACTCTTCAGTGtaatt                    >  2:225477/1‑148 (MQ=255)
                                                                                                                                  gaaaagaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGAAGGCAAAAGTGAAATCGATATCCAGGCGCACCCCGTCAGCCTCTTCAGTGAAAGt                    >  2:228864/1‑150 (MQ=255)
                                                                                                                                    aaagaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCg                  >  2:1275897/1‑150 (MQ=255)
                                                                                                                                      agaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCggg                >  1:1285958/1‑150 (MQ=255)
                                                                                                                                      agaaaaCGCCCATGACAGCAACTGACATGTGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTTGAAAATCAGCATTTCAGCTTCGCAGTCAAAAGTGAATTCGATATCCAGGCGCACACCGTCAGACTCTTCAGTTAAAGTCggg                <  2:1223861/150‑1 (MQ=255)
                                                                                                                                       gaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAAGGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGATATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCgggg               <  1:710667/150‑1 (MQ=255)
                                                                                                                                       gaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCCGGCGCACCCCGTCAGCCTCTTCAGTGAAAGTCgggg               >  1:558339/1‑150 (MQ=255)
                                                                                                                                        aaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGt              <  1:793078/150‑1 (MQ=255)
                                                                                                                                           aCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAt           <  1:1266177/150‑1 (MQ=255)
                                                                                                                                            cGCCCCTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAATTCGATATCCAGGCGCCCACCGTCAGACTCTTCAGTGAAAGTCGGGGGGAtt          >  2:74077/1‑150 (MQ=255)
                                                                                                                                            cGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGCCCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtt          >  1:660339/1‑150 (MQ=255)
                                                                                                                                             gCCCCTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAttt         >  2:677775/1‑150 (MQ=255)
                                                                                                                                              ccccTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtttt        >  1:1111683/1‑150 (MQ=255)
                                                                                                                                              ccccTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCGTCAGTGAAAGTCGGGGTGAtttt        >  2:140472/1‑150 (MQ=255)
                                                                                                                                              ccccTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTACCGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGGGAAAGCGATCTCCAGGCGCACACCGGCAGACTCTTCAGGGCAAGTCGGGGTGAtttt        >  1:764712/1‑150 (MQ=255)
                                                                                                                                              ccccTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtttt        >  1:579906/1‑150 (MQ=255)
                                                                                                                                                  tGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCaa    <  1:580217/150‑1 (MQ=255)
                                                                                                                                                    ccAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAgg  <  1:446326/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCATAACACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGG  >  NC_000913/2460466‑2460763

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: