Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,460,624 A→G 8.8% intergenic (‑155/+26) 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,6240AG 8.8% 217.1 / 11.1 80intergenic (‑155/+26)fadI/yfcZbeta‑ketoacyl‑CoA thiolase, anaerobic, subunit/UPF0381 family protein
Reads supporting (aligned to +/- strand):  ref base A (39/33);  new base G (4/3);  total (44/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 = 5.03e-01

TCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGAT  >  NC_000913/2460475‑2460770
                                                                                                                                                     |                                                                                                                                                  
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:478285/150‑6 (MQ=255)
              ataaataaAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGTCAGTTGCTGGCAGGGGcgttttct                                                                                                                                      <  2:1214423/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)
                                                                                                                                                tgcTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAg    <  2:1111683/148‑1 (MQ=255)
                                                                                                                                                  cTGACAGGGTCGTACGTACCCGAGAGGTATTAACGCAGACCCAGCTTGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACTGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGAt  <  1:166918/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TCCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGAT  >  NC_000913/2460475‑2460770

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: