Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,460,624 A→G 9.6% intergenic (‑155/+26) fadI ← / ← yfcZ 3‑ketoacyl‑CoA thiolase FadI/DUF406 domain‑containing protein YfcZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,460,6240AG 9.6% 194.6 / 10.9 73intergenic (‑155/+26)fadI/yfcZ3‑ketoacyl‑CoA thiolase FadI/DUF406 domain‑containing protein YfcZ
Reads supporting (aligned to +/- strand):  ref base A (37/29);  new base G (4/3);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-01

CCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTC  >  NC_000913/2460476‑2460766
                                                                                                                                                    |                                                                                                                                              
ccccTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGAc                                                                                                                                               >  1:1030990/1‑150 (MQ=255)
ccccTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGAc                                                                                                                                               >  1:262232/1‑150 (MQ=255)
     gTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGTCAGTTGCTGGCAgggg                                                                                                                                          <  2:173228/150‑1 (MQ=255)
     gTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAgggg                                                                                                                                          >  2:1061681/1‑150 (MQ=255)
     gTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGACAGAAACTGACAgggg                                                                                                                                          >  2:989359/1‑150 (MQ=255)
      tAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGc                                                                                                                                         >  1:358440/1‑150 (MQ=255)
      tAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGc                                                                                                                                         >  2:265045/1‑150 (MQ=255)
       aaaaaaaTAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcg                                                                                                                                        >  1:1127348/1‑150 (MQ=255)
        aaaaaaTAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgt                                                                                                                                       >  2:1343559/1‑150 (MQ=255)
         aaaaaTAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgta                                                                                                                                      <  2:1273981/150‑1 (MQ=255)
           aaataaatAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgtacg                                                                                                                                    >  2:360019/1‑150 (MQ=255)
              taaataaaGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAACCGCCCCTGTCAGTTGCTGGCAGGGGcg                                                                                                                                        <  1:985056/143‑1 (MQ=255)
              taaataaaGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGTCAGTTGCTGGCAGGGGcg                                                                                                                                        >  2:1337471/1‑143 (MQ=255)
              taaataaaGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGTCAGTTGCTGGCAGGGGcg                                                                                                                                        <  1:1337471/143‑1 (MQ=255)
              taaataaaGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGTCAGTTGCTGGCAGGGGcg                                                                                                                                        >  2:985056/1‑143 (MQ=255)
                   aaaGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCgat                                                                                                                            >  1:751197/1‑149 (MQ=255)
                    aaGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCgaga                                                                                                                           >  1:1240462/1‑150 (MQ=255)
                      gTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAgg                                                                                                                         <  2:112803/150‑1 (MQ=255)
                      gTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAgg                                                                                                                         <  1:778907/150‑1 (MQ=255)
                          tCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAAtt                                                                                                                     >  2:287070/1‑150 (MQ=255)
                           cTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTa                                                                                                                    >  1:1361295/1‑150 (MQ=255)
                            tGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgta                                                                                                                                      <  1:267261/131‑1 (MQ=255)
                            tGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGcgta                                                                                                                                      >  2:267261/1‑131 (MQ=255)
                             gACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAAc                                                                                                                  <  2:751197/150‑1 (MQ=255)
                                  gATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGa                                                                                                             >  2:721122/1‑150 (MQ=255)
                                     cATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGAccc                                                                                                          >  2:486882/1‑150 (MQ=255)
                                      aTAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCa                                                                                                         >  1:1127720/1‑150 (MQ=255)
                                       tAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAg                                                                                                        >  1:743631/1‑150 (MQ=255)
                                         gTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCTGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCAAGCt                                                                                                      >  1:407548/1‑150 (MQ=255)
                                             tAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGaa                                                                                                  >  1:335732/1‑150 (MQ=255)
                                              aaCCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGaaa                                                                                                 <  2:356377/150‑1 (MQ=255)
                                               aCCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGaaaa                                                                                                >  2:10757/1‑150 (MQ=255)
                                                 cATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATc                                                                                              <  2:1127348/150‑1 (MQ=255)
                                                     tttttACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCa                                                                                          <  2:358440/150‑1 (MQ=255)
                                                      ttttACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCAt                                                                                         <  1:287070/150‑1 (MQ=255)
                                                           cATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAg                                                                                    >  1:1389339/1‑150 (MQ=255)
                                                            aTTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGc                                                                                   >  1:1414428/1‑150 (MQ=255)
                                                             tttAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCt                                                                                  >  1:841020/1‑150 (MQ=255)
                                                               tAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTc                                                                                >  2:1267619/1‑150 (MQ=255)
                                                                aGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCg                                                                               <  2:266423/150‑1 (MQ=255)
                                                                  ccAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCa                                                                             <  1:486882/150‑1 (MQ=255)
                                                                   cAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAg                                                                            <  2:262232/150‑1 (MQ=255)
                                                                        ggAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCaaa                                                                       >  2:259801/1‑150 (MQ=255)
                                                                          agaAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAg                                                                     <  1:1343559/150‑1 (MQ=255)
                                                                          agaAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAg                                                                     <  2:259589/150‑1 (MQ=255)
                                                                               aGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGaaa                                                                >  2:433243/1‑150 (MQ=255)
                                                                               aGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGaaa                                                                >  1:825438/1‑150 (MQ=255)
                                                                                 ggAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATc                                                              <  1:552973/150‑1 (MQ=255)
                                                                                    aaGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGat                                                           <  2:335732/150‑1 (MQ=255)
                                                                                     aGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGata                                                          <  1:563903/150‑1 (MQ=255)
                                                                                          gAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAg                                                     <  2:243684/150‑1 (MQ=255)
                                                                                           aGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAgg                                                    <  2:1030990/150‑1 (MQ=255)
                                                                                               aTGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTAAGTACCCGAGAGGAATTAACGCAGACACAGCTGGAAAATCAGCATTTCAGCTTCGAAGGCAAAAGTGAAATCGATATCCAGGCGca                                                >  2:78653/1‑150 (MQ=255)
                                                                                                  acacAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCacac                                             <  2:1389339/150‑1 (MQ=255)
                                                                                                  acacAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTAAGTAACAGAGAGGAAATAAACAAGAACAATATATAAAATAAGCATAACAAAATATAAAACAAAAGTGAAATCGATATCCAGGCGCacac                                             >  2:60949/1‑150 (MQ=255)
                                                                                                       gagagaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCa                                        >  1:1151790/1‑150 (MQ=255)
                                                                                                        agagaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAg                                       >  2:483109/1‑150 (MQ=255)
                                                                                                           gaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGAct                                    >  2:327656/1‑150 (MQ=255)
                                                                                                               agaagaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTc                                <  1:265045/150‑1 (MQ=255)
                                                                                                               agaagaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTc                                <  2:46382/150‑1 (MQ=255)
                                                                                                                 aagaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAg                              <  1:955536/150‑1 (MQ=255)
                                                                                                                   gaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACTGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCATGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTg                            <  2:481543/150‑1 (MQ=255)
                                                                                                                    aagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGa                           <  2:1240462/150‑1 (MQ=255)
                                                                                                                     agagAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGaa                          >  2:821158/1‑150 (MQ=255)
                                                                                                                        gaaaagaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGt                       >  1:1064373/1‑150 (MQ=255)
                                                                                                                             gaaaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCgggg                  <  1:1061681/150‑1 (MQ=255)
                                                                                                                               aaaCGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTg                >  2:582610/1‑150 (MQ=255)
                                                                                                                                  cGCCCCTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAtt             >  1:996614/1‑150 (MQ=255)
                                                                                                                                   gCCCCTGTCAGTTGCTGGCAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAttt            >  2:564792/1‑150 (MQ=255)
                                                                                                                                   gCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGAttt            >  2:681360/1‑150 (MQ=255)
                                                                                                                                          ccAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAgg     <  1:828098/150‑1 (MQ=255)
                                                                                                                                          ccAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAgg     <  2:1361295/150‑1 (MQ=255)
                                                                                                                                             gCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTc  <  1:721122/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CCCCTGTAAAAAAATAAATAAAGTGGTCTGACCTGATCATAGTCTTAACCATTTTTTTACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTC  >  NC_000913/2460476‑2460766

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