Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,457,031 G→A A710V (GCA→GTA)  fadJ ← 3‑hydroxyacyl‑CoA dehydrogenase FadJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,457,0310GA90.5% 177.3 / 8.1 63A710V (GCA→GTA) fadJ3‑hydroxyacyl‑CoA dehydrogenase FadJ
Reads supporting (aligned to +/- strand):  ref base G (5/1);  new base A (34/23);  total (39/24)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

TCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGA  >  NC_000913/2456898‑2457163
                                                                                                                                     |                                                                                                                                    
tCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTcagttacagttg                                                                                                                                >  2:454329/1‑140 (MQ=255)
 ccATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAgttgt                                                                                                                               <  2:506953/140‑1 (MQ=255)
   aTTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGttt                                                                                                                             >  1:265701/1‑140 (MQ=255)
   aTTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGttt                                                                                                                             >  1:69549/1‑140 (MQ=255)
        ccGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCaa                                                                                                                        <  1:454329/140‑1 (MQ=255)
           aaTGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCaaaaa                                                                                                                     >  2:471845/1‑140 (MQ=255)
           aaTGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCaaaaa                                                                                                                     >  1:331253/1‑140 (MQ=255)
           aaTGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCaaaaa                                                                                                                     >  2:421340/1‑140 (MQ=255)
            aTGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAAc                                                                                                                    >  2:145708/1‑140 (MQ=255)
                      gTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCAcg                                                                                                          <  2:727902/140‑1 (MQ=255)
                          tttACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCAcgcgcg                                                                                                      <  2:69549/140‑1 (MQ=255)
                          tttACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCAcgcgcg                                                                                                      <  2:265701/140‑1 (MQ=255)
                          tttACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCAcgcgcg                                                                                                      <  1:471845/140‑1 (MQ=255)
                               cTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCAt                                                                                                 >  2:500390/1‑140 (MQ=255)
                                         tAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCaaa                                                                                       >  1:149663/1‑140 (MQ=255)
                                                aGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGc                                                                                >  2:85688/1‑140 (MQ=255)
                                                aGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGc                                                                                >  2:159159/1‑140 (MQ=255)
                                                aGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGc                                                                                >  1:263208/1‑140 (MQ=255)
                                                  tCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCaa                                                                              >  2:748118/1‑140 (MQ=255)
                                                  tCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCa                                                                               >  1:709067/1‑139 (MQ=255)
                                                  tCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCa                                                                               <  2:709067/139‑1 (MQ=255)
                                                            tttACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAAc                                                                    >  1:233798/1‑140 (MQ=255)
                                                            tttACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAAc                                                                    <  1:748118/140‑1 (MQ=255)
                                                              tACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACgg                                                                  <  2:233798/140‑1 (MQ=255)
                                                                cATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGa                                                                <  2:149663/140‑1 (MQ=255)
                                                                    ttGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCa                                                            >  2:166280/1‑140 (MQ=255)
                                                                            cacTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGt                                                    >  1:346509/1‑140 (MQ=255)
                                                                                 cGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCaa                                               >  1:189894/1‑140 (MQ=255)
                                                                                   ccATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGt                                             <  2:346509/140‑1 (MQ=255)
                                                                                   ccATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGt                                             >  1:822911/1‑140 (MQ=255)
                                                                                   ccATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGt                                             <  2:822911/140‑1 (MQ=255)
                                                                                      ttCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTg                                                       >  2:633780/1‑127 (MQ=255)
                                                                                      ttCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTg                                                       <  1:633780/127‑1 (MQ=255)
                                                                                      ttCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGt                                          >  1:597880/1‑140 (MQ=255)
                                                                                      ttCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGt                                          >  1:159438/1‑140 (MQ=255)
                                                                                         aGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGc                                       >  1:484373/1‑140 (MQ=255)
                                                                                          gcgcGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCa                                      <  1:145708/140‑1 (MQ=255)
                                                                                          gcgcGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCa                                      <  2:263208/140‑1 (MQ=255)
                                                                                          gcgcGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCa                                      <  2:159438/140‑1 (MQ=255)
                                                                                           cgcgGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCat                                     >  1:419482/1‑140 (MQ=255)
                                                                                           cgcgGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCat                                     >  2:122593/1‑140 (MQ=255)
                                                                                            gcgGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCatt                                    <  1:421340/140‑1 (MQ=255)
                                                                                                          cTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTc                      >  1:333701/1‑140 (MQ=255)
                                                                                                          cTTTGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTc                      >  1:313444/1‑140 (MQ=255)
                                                                                                            ttGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGc                    >  2:378808/1‑140 (MQ=255)
                                                                                                            ttGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGc                    >  1:700200/1‑140 (MQ=255)
                                                                                                             tGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGcc                   >  2:340461/1‑140 (MQ=255)
                                                                                                             tGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGcc                   <  2:700200/140‑1 (MQ=255)
                                                                                                             tGACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGcc                   >  2:763750/1‑140 (MQ=255)
                                                                                                              gACCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGccc                  >  1:399625/1‑140 (MQ=255)
                                                                                                               aCCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAACCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCcg                 >  2:279166/1‑140 (MQ=255)
                                                                                                               aCCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCcg                 >  1:431524/1‑140 (MQ=255)
                                                                                                               aCCTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCcg                 >  2:146899/1‑140 (MQ=255)
                                                                                                                ccTTCTTATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCcgc                <  2:399625/140‑1 (MQ=255)
                                                                                                                  ttcttATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCgcgc              <  2:313444/140‑1 (MQ=255)
                                                                                                                    cttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCg            >  1:171354/1‑140 (MQ=255)
                                                                                                                    cttATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCg            >  2:535085/1‑140 (MQ=255)
                                                                                                                      tATTGCAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgag          <  1:763750/140‑1 (MQ=255)
                                                                                                                           cAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAAt     <  1:340461/140‑1 (MQ=255)
                                                                                                                           cAGGTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAAt     <  1:4905/140‑1 (MQ=255)
                                                                                                                              gtaagttgcagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGa  >  2:502265/4‑140 (MQ=255)
                                                                                                                              gTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGcc                   >  2:414018/1‑123 (MQ=255)
                                                                                                                              gTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGcc                   <  1:414018/123‑1 (MQ=255)
                                                                                                                              gTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGa  >  2:444660/1‑140 (MQ=255)
                                                                                                                              gTCAGTTACAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGa  >  2:208923/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                    
TCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGA  >  NC_000913/2456898‑2457163

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