Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,460,619 A→T 8.4% intergenic (‑150/+31) 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,6190AT 8.4% 207.0 / 7.0 81intergenic (‑150/+31)fadI/yfcZbeta‑ketoacyl‑CoA thiolase, anaerobic, subunit/UPF0381 family protein
Reads supporting (aligned to +/- strand):  ref base A (40/33);  new base T (5/2);  total (46/35)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

ACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGG  >  NC_000913/2460473‑2460763
                                                                                                                                                  |                                                                                                                                                
acTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACt                                                                                                                                               <  2:910294/150‑1 (MQ=255)
acTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACt                                                                                                                                               >  1:1090677/1‑150 (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                                                                                                                        >  1:1251050/1‑150 (MQ=255)
                       aaGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCgaga                                                                                                                        >  2:501234/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:446135/1‑150 (MQ=255)
                                                   ccATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAAt                                                                                            >  1:199949/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)
                                                                                                                                                  |                                                                                                                                                
ACTCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGG  >  NC_000913/2460473‑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: