Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,420,793 (G)6→7 coding (110/423 nt) ydaT → Rac prophage; protein YdaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,420,7871.G98.2% 196.6 / ‑3.4 55coding (104/423 nt)ydaTRac prophage; protein YdaT
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (23/31);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGG  >  NC_000913/1420652‑1420918
                                                                                                                                        |                                                                                                                                      
tCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAgggg                                                                                                                                   <  2:1587/140‑1 (MQ=255)
tCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAgggg                                                                                                                                   >  1:887076/1‑140 (MQ=255)
    aCCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAgggggggc                                                                                                                               >  2:1024801/1‑138 (MQ=255)
    aCCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAgggggggc                                                                                                                               >  2:1069971/1‑138 (MQ=255)
     ccTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAgggggggcg                                                                                                                              <  2:18798/140‑4 (MQ=255)
      cTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAgggggggcgg                                                                                                                             >  2:726360/1‑136 (MQ=255)
         aaaCATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGa                                                                                                                           >  2:304511/1‑140 (MQ=255)
           aCATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAAt                                                                                                                         <  2:887076/140‑1 (MQ=255)
            cATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAAtt                                                                                                                        <  2:551222/140‑1 (MQ=255)
              tCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAgggggggcgg                                                                                                                             >  2:523507/1‑128 (MQ=255)
              tCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAgggggggcgg                                                                                                                             <  1:523507/132‑5 (MQ=255)
              tCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAgggggggc                                                                                                                               <  2:48921/130‑3 (MQ=255)
              tCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAgggggggc                                                                                                                               >  1:48921/1‑128 (MQ=255)
              tCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGTGGGGGCGGCGAATTGc                                                                                                                      >  2:1253779/1‑140 (MQ=255)
                 agagGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGcccc                                                                                                                   <  1:1070519/140‑1 (MQ=255)
                 agagGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGcccc                                                                                                                   <  2:1190495/140‑1 (MQ=255)
                 agagGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGcccc                                                                                                                   >  2:1264063/1‑140 (MQ=255)
                 agagGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGcccc                                                                                                                   <  1:304511/140‑1 (MQ=255)
                  gagGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCt                                                                                                                  <  1:1006119/140‑1 (MQ=255)
                       aaaTGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTc                                                                                                             <  2:1285221/140‑1 (MQ=255)
                       aaaTGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTc                                                                                                             >  2:978668/1‑140 (MQ=255)
                         aTGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCa                                                                                                           >  1:1322585/1‑140 (MQ=255)
                             ttCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCa                                                                                                       <  1:1264063/140‑1 (MQ=255)
                             ttCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCa                                                                                                       <  2:198138/140‑1 (MQ=255)
                             ttCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCa                                                                                                       >  1:690394/1‑140 (MQ=255)
                               cATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAg                                                                                                     <  2:824617/140‑1 (MQ=255)
                               cATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAg                                                                                                     <  2:690394/140‑1 (MQ=255)
                                aTGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAgg                                                                                                    >  2:1074395/1‑140 (MQ=255)
                                          aGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACaa                                                                                          >  2:846700/1‑140 (MQ=255)
                                              tGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcag                                                                                      <  1:1069971/140‑1 (MQ=255)
                                              tGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcag                                                                                      <  1:978668/140‑1 (MQ=255)
                                               gAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAg                                                                                                     >  1:102497/1‑124 (MQ=255)
                                               gAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAg                                                                                                     <  2:102497/124‑1 (MQ=255)
                                               gAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcagc                                                                                     <  2:1322585/140‑1 (MQ=255)
                                                aGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcagca                                                                                    <  2:547478/140‑1 (MQ=255)
                                                  cacaTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGa                                                                                 <  2:350196/140‑1 (MQ=255)
                                                  cacaTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGa                                                                                 <  2:624372/140‑1 (MQ=255)
                                                     aTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGa                                                                                 >  1:1213972/1‑137 (MQ=255)
                                                     aTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGa                                                                                 <  2:1213972/137‑1 (MQ=255)
                                                     aTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAata                                                                              >  1:188391/1‑140 (MQ=255)
                                                     aTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAat                                                                               <  2:953865/139‑1 (MQ=255)
                                                     aTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAat                                                                               >  1:953865/1‑139 (MQ=255)
                                                             aGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCAc                                                                      <  1:846700/140‑1 (MQ=255)
                                                                gttgttTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCgt                                                                   <  1:1254117/140‑1 (MQ=255)
                                                                gttgttTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCgt                                                                   <  1:79546/140‑1 (MQ=255)
                                                                       gCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGa                                                            <  2:1001009/140‑1 (MQ=255)
                                                                       gCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGa                                                            <  1:1160969/140‑1 (MQ=255)
                                                                        cccTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGaa                                                           >  1:429121/1‑140 (MQ=255)
                                                                                gcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGaaa                                                   >  2:649252/1‑140 (MQ=255)
                                                                                                caTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCa                                        >  2:597256/1‑134 (MQ=255)
                                                                                                caTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCa                                        <  1:597256/134‑1 (MQ=255)
                                                                                                         ccAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCa                         >  2:395807/1‑140 (MQ=255)
                                                                                                            aTGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGaa                      >  2:1319385/1‑140 (MQ=255)
                                                                                                             tGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAg                     >  2:390824/1‑140 (MQ=255)
                                                                                                             tGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAg                     >  2:242898/1‑140 (MQ=255)
                                                                                                              gCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGc                    >  1:711279/1‑140 (MQ=255)
                                                                                                              gCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAAGATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGc                    >  2:1223232/1‑140 (MQ=255)
                                                                                                                  ttACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGAt                >  2:512759/1‑140 (MQ=255)
                                                                                                                       gaagaaTATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAg           <  2:1103889/140‑1 (MQ=255)
                                                                                                                        aagaaTATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCTAAGAGAAAAAATTCAGAAGCt                   <  1:1165140/131‑1 (MQ=255)
                                                                                                                        aagaaTATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCTAAGAGAAAAAATTCAGAAGCt                   >  2:1165140/1‑131 (MQ=255)
                                                                                                                        aagaaTATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGc          <  2:865720/140‑1 (MQ=255)
                                                                                                                        aagaaTATCTGCGCCAGGGGGGGCGCCGATTTGCCCCTGTTCCCAGCTAAGACCTGGACCATCAGCACGAATATCAACAACAGTTGATTGAACGGTAAACCGAAACGCTAAAGGGAAAAAATTAAAAGGCTAATCCTAGc          >  2:364292/1‑140 (MQ=255)
                                                                                                                           aaTATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAAt       <  1:512759/140‑1 (MQ=255)
                                                                                                                             tatCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTc     >  1:767181/1‑140 (MQ=255)
                                                                                                                                cTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTgg  <  1:395807/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
TCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGG  >  NC_000913/1420652‑1420918

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: