Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,320,735 G→A 5.1% G171G (GGC→GGT phnH ← ribophosphonate triphosphate synthase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,320,7350GA 5.1% 260.3 / 2.8 78G171G (GGC→GGTphnHribophosphonate triphosphate synthase subunit
Reads supporting (aligned to +/- strand):  ref base G (37/36);  new base A (2/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAA  >  NC_000913/4320589‑4320876
                                                                                                                                                  |                                                                                                                                             
gCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGc                                                                                                                                              <  2:320231/148‑1 (MQ=255)
          cGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGaa                                                                                                                                    <  4:179267/148‑1 (MQ=255)
          cGGCGCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGaa                                                                                                                                    <  2:3414/148‑1 (MQ=255)
             ctctctTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCg                                                                                                                                        <  4:106166/141‑1 (MQ=255)
             ctctctTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCg                                                                                                                                        >  3:106166/1‑141 (MQ=255)
                 ctTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATACCGAGCGGGAACGGATGc                                                                                                                             >  3:282462/1‑148 (MQ=255)
                           gtgttcggctTCTATCGCCTTCTCGCCCACTTTCACTGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATACCGATCGGGAACGGATCCGGGCGCTCgt                                                                                                                   <  4:282462/143‑2 (MQ=255)
                           gTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCgg                                                                                                                   <  1:203520/148‑1 (MQ=255)
                                               tctcGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCAc                                                                                               >  1:386842/1‑148 (MQ=255)
                                                          ttCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGa                                                                                                     >  5:148330/1‑131 (MQ=255)
                                                          ttCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGa                                                                                                     <  6:148330/131‑1 (MQ=255)
                                                          ttCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTg                                                                                    >  2:63532/1‑148 (MQ=255)
                                                          ttCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTg                                                                                    >  5:3340/1‑148 (MQ=255)
                                                          ttCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATACCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTg                                                                                    >  8:380877/1‑148 (MQ=255)
                                                           tCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGc                                                                                   >  4:242493/1‑148 (MQ=255)
                                                            cACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCg                                                                                  >  6:384275/1‑148 (MQ=255)
                                                             aCGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCg                                                                                  >  8:353281/1‑147 (MQ=255)
                                                              cGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGc                                                                                                 >  1:53606/1‑131 (MQ=255)
                                                              cGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGc                                                                                                 <  2:53606/131‑1 (MQ=255)
                                                                gCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGc                                                                              >  4:166577/1‑148 (MQ=255)
                                                                    cGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATc                                                                          <  1:63532/148‑1 (MQ=255)
                                                                    cGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATc                                                                          <  7:353281/148‑1 (MQ=255)
                                                                     gTACATCAGCCCACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCGCGCCACAGGTCTGGATCAGGTCGATTCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGa                                                                                                     <  5:64598/120‑1 (MQ=255)
                                                                     gTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGa                                                                            <  2:235379/145‑1 (MQ=255)
                                                                     gTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGa                                                                            >  1:235379/1‑145 (MQ=255)
                                                                      tACATCAGCACACATCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCGCGGTGAGTTCGTGCAGAATGCCCTCCGGCAGCTGCGGAGCGATCag                                                                        >  2:302266/1‑147 (MQ=255)
                                                                       aCATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCa                                                                                        <  2:235203/131‑1 (MQ=255)
                                                                       aCATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCa                                                                                        >  1:235203/1‑131 (MQ=255)
                                                                        cATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCAtt                                                                       <  5:384275/147‑1 (MQ=255)
                                                                          tCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGt                                                                    >  8:240065/1‑148 (MQ=255)
                                                                           cAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGtt                                                                                                             >  5:195647/1‑106 (MQ=255)
                                                                           cAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGtt                                                                                                             <  6:195647/106‑1 (MQ=255)
                                                                                cacCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCg                                                              <  2:386842/148‑1 (MQ=255)
                                                                                 acCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCgg                                                             >  1:126767/1‑148 (MQ=255)
                                                                                 acCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGACAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTAGGTGAg                                                                                                               >  5:135420/1‑98 (MQ=255)
                                                                                  ccTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGc                                                            <  2:126767/148‑1 (MQ=255)
                                                                                   cTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCa                                                           <  6:3340/148‑1 (MQ=255)
                                                                                     ccACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCg                                                                     >  4:125725/1‑136 (MQ=255)
                                                                                     ccACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCg                                                                     <  3:125725/136‑1 (MQ=255)
                                                                                     ccACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCaa                                                          >  3:339768/1‑147 (MQ=255)
                                                                                      cacaTGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCa                                                           <  4:339768/145‑1 (MQ=255)
                                                                                         aTGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGa                                                                                                     <  1:121288/100‑1 (MQ=255)
                                                                                         aTGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGa                                                                                                     >  2:121288/1‑100 (MQ=255)
                                                                                          tGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCgg                                                             >  5:255504/1‑139 (MQ=255)
                                                                                          tGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCgg                                                             <  6:255504/139‑1 (MQ=255)
                                                                                                 ttCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATACCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGt                                             <  7:380877/148‑1 (MQ=255)
                                                                                                  tCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACgg                                                                                                                                 <  1:268399/63‑1 (MQ=255)
                                                                                                  tCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACgg                                                                                                                                 >  2:268399/1‑63 (MQ=255)
                                                                                                    gcgGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAg                                          >  1:398900/1‑148 (MQ=255)
                                                                                                            gCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCa                                               <  7:67932/135‑1 (MQ=255)
                                                                                                            gCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCa                                               >  8:67932/1‑135 (MQ=255)
                                                                                                            gCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAAc                                   <  2:398900/147‑1 (MQ=255)
                                                                                                            gCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGGCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAAc                                   <  4:108817/147‑1 (MQ=255)
                                                                                                              cagcagGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGgcgc                                      >  2:182947/1‑142 (MQ=255)
                                                                                                              cagcagGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGgcgc                                      <  1:182947/142‑1 (MQ=255)
                                                                                                                gcagGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCg                              <  3:242493/148‑1 (MQ=255)
                                                                                                                gcagGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCg                              >  6:258996/1‑148 (MQ=255)
                                                                                                                     cgcTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGa                                                                                                     >  3:302854/1‑72 (MQ=255)
                                                                                                                     cgcTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGa                                                                                                     <  4:302854/72‑1 (MQ=255)
                                                                                                                     cgcTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTccgc                         <  7:240065/148‑1 (MQ=255)
                                                                                                                               aGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCt               >  6:360813/1‑148 (MQ=255)
                                                                                                                                  tCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGc            >  7:303503/1‑148 (MQ=255)
                                                                                                                                   cAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGc            >  1:127547/1‑147 (MQ=255)
                                                                                                                                      gATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTc                                                                      >  2:290417/1‑86 (MQ=255)
                                                                                                                                      gATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTc                                                                      <  1:290417/86‑1 (MQ=255)
                                                                                                                                      gATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAAc                                   >  3:194553/1‑121 (MQ=255)
                                                                                                                                      gATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAAc                                   <  4:194553/121‑1 (MQ=255)
                                                                                                                                       aTCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAAc                                   <  5:154004/120‑1 (MQ=255)
                                                                                                                                       aTCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAAc                                   >  6:154004/1‑120 (MQ=255)
                                                                                                                                       aTCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGa          <  4:188643/145‑1 (MQ=255)
                                                                                                                                       aTCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGa          >  3:188643/1‑145 (MQ=255)
                                                                                                                                         cAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGt                                             <  4:36595/108‑1 (MQ=255)
                                                                                                                                         cAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGt                                             >  3:36595/1‑108 (MQ=255)
                                                                                                                                         cAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTCAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTg      >  8:126988/1‑147 (MQ=255)
                                                                                                                                            gTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTaaa  <  3:166577/148‑1 (MQ=255)
                                                                                                                                            gTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTaaa  <  5:360813/148‑1 (MQ=255)
                                                                                                                                             tCGATGCCGCGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCgc                                                                            >  1:65343/1‑72 (MQ=255)
                                                                                                                                             tCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGa                                                                            <  2:65343/73‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAA  >  NC_000913/4320589‑4320876

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: