Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,950,342 Δ1 bp coding (21/99 nt) ilvL → ilvG operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,3420G.97.0% 276.1 / 1.5 67coding (21/99 nt)ilvLilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (32/32);  minor base C (2/0);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGAC  >  NC_000913/3950211‑3950479
                                                                                                                                   |                                                                                                                                         
ttagattattAAACGCATCATAAAAAGCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACATCCCTTCTACGAGTATTAGCCTTGTc                                                                                                                               <  2:584092/139‑1 (MQ=255)
tGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTc                                                                                                                               <  2:717941/143‑1 (MQ=255)
  aaaTTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGAATTAGCCTGGTCGt                                                                                                                             >  1:176401/1‑143 (MQ=255)
    attattAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGa                                                                                                                           <  2:621138/143‑1 (MQ=255)
    attattAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGa                                                                                                                           >  1:920449/1‑143 (MQ=255)
       attaAACGCATCATAAAAATCGGCCAAAAAATATCTTCTACTATTGACAAAACATAAGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTa                                                                                                                        <  2:811809/143‑1 (MQ=255)
         taAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGc                                                                                                                      <  2:569679/143‑1 (MQ=255)
         taAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTATACGAGTATTAGCCTGGTCGTGCTTAGc                                                                                                                      >  2:388637/1‑143 (MQ=255)
            aCGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTGGTGATTAGCgtg                                                                                                                   >  1:446665/1‑143 (MQ=255)
            aCGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGCATTAGCCTGGGCGTGATTAGCgtg                                                                                                                   >  2:158289/1‑143 (MQ=255)
             cGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCGACGAGTATTAGCCTGGTCGTGATTAGCgtgg                                                                                                                  >  1:1047772/1‑143 (MQ=255)
              gCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGACCTTCTACGACTATTAGCCTGGTCGTGATTAGCgtggt                                                                                                                 >  2:499438/1‑143 (MQ=255)
                atcatAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCgtggtgg                                                                                                               >  1:815849/1‑143 (MQ=255)
                   agaaacaTCGGCCAAAAAATATCTTGTACTATTTACAAAACCTCTGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGa                                                                                                            <  1:567986/141‑1 (MQ=255)
                          tCGGCCAAAAAATAGCTTGTACTATTTACAAACTAAATGAAAACTCGCTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGGCGTGATTAGCGTGGTGGTGATtattat                                                                                                     <  2:595701/143‑1 (MQ=255)
                              ccAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCa                                                                                                 >  2:1133947/1‑143 (MQ=255)
                                aaaaaaTATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGAAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGGGATTATTATaccacc                                                                                               >  2:18436/1‑143 (MQ=255)
                                      tatCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCgg                                                                                         <  1:868575/143‑1 (MQ=255)
                                      gatCTTGTACTATTTACGAAACATATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGGCAGAATGTCAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCgg                                                                                         <  1:162784/142‑1 (MQ=255)
                                            gTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTCTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTACTATACCACCGTGCGGGGCTGc                                                                                   >  2:584656/1‑143 (MQ=255)
                                            gTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGc                                                                                   >  2:1160780/1‑143 (MQ=255)
                                             tactacGTACAAAACCTACGGTAACTCATTAGGCATCCCGTCGAACAAGATGCAACAAAATACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCa                                                                                  <  1:768781/143‑1 (MQ=255)
                                               cTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTGTTATCCCACCGTGCGGCGCTGCAcc                                                                                >  2:1343553/1‑142 (MQ=255)
                                                tATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACtt                                                                               >  1:461760/1‑143 (MQ=255)
                                                  tttACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCAcattg                                                                             >  1:1122186/1‑139 (MQ=255)
                                                  tttACAAAACCTATGGGAACTCTTTAGGCATTCCTTCGAACACGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTgg                                                                             <  2:76272/143‑1 (MQ=255)
                                                    tACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGAc                                                                           <  2:558547/143‑1 (MQ=255)
                                                      cAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTCCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGa                                                                         <  1:147948/143‑1 (MQ=255)
                                                      cAAAACATACGGGAACGCTATAGGCAGGCCTTCGAAGAAGATGCAAGAAAAGACAAACTGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGa                                                                         <  1:319364/143‑1 (MQ=255)
                                                       aaaaCCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCGCTTCTACGAGTATTAGCCTGGTAGTGATTAGCGTGGTGGAGAGTATTATCCCACCGTGCGGGGCTGCACTTGGACGAg                                                                        >  1:411702/1‑143 (MQ=255)
                                                       aaaaCCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAg                                                                        >  1:1183459/1‑143 (MQ=255)
                                                       aaaaCCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAg                                                                        >  2:30616/1‑143 (MQ=255)
                                                       aaaaCCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATACCACCGTGCGCGGCTGCACTTGGACGac                                                                        >  2:874109/1‑142 (MQ=255)
                                                         aaCCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCACACCGTGCGGGGCTGCACTTGGACGAGGa                                                                      >  1:683605/1‑143 (MQ=255)
                                                           cagagggTAACCCGTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGaaa                                                                    <  2:1159513/138‑1 (MQ=255)
                                                                ggTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAGAAGACAAAATGACAGCCCTTCTACGAGTCATACCCAGGCGTTGTTTACCGTCGTGGTGATGATTTTCCCCCCGTGCTGGGCTGCTCTTGGCCGCGGAAAAGcac                                                               >  2:1097413/1‑141 (MQ=255)
                                                                   aaCTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTaga                                                            <  2:1215959/143‑1 (MQ=255)
                                                                   aaCTCTCTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTaga                                                            <  1:1065088/143‑1 (MQ=255)
                                                                    aCTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGCAAGGCTTagag                                                           >  1:122442/1‑143 (MQ=255)
                                                                        tttAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCa                                                       >  1:424926/1‑143 (MQ=255)
                                                                          tAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAg                                                     <  2:713595/143‑1 (MQ=255)
                                                                           aGGCATTCCTTCGAACAAGATGCAAGAAAAGACAACATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGc                                                    <  1:117178/143‑1 (MQ=255)
                                                                                ttCCTTGGAACAAGATGCAAGAAAAGACAGAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTaa                                               <  1:1311826/143‑1 (MQ=255)
                                                                                    ttCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGaa                                           >  1:1078917/1‑143 (MQ=255)
                                                                                    ttCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTCAACGaa                                           >  2:1239717/1‑143 (MQ=255)
                                                                                       gAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGTTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGTCTTAGAGATCAAGCCTTAACGAACTa                                        >  1:776825/1‑143 (MQ=255)
                                                                                       gAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTa                                        <  1:1066539/143‑1 (MQ=255)
                                                                                        aatcaGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTaa                                       <  1:830688/139‑1 (MQ=255)
                                                                                        aaCAAGCTGCAACGAAAGACACAATGACAGCCCTTCTACTCATATTAGCCTGGTCGTGATTAGCGTGGTGGTTATTATTACCCCACCGTGCGGGGCTGCACTTGGACCAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTaa                                       <  2:1192598/143‑1 (MQ=255)
                                                                                        aaCAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTaa                                       >  1:1011270/1‑143 (MQ=255)
                                                                                          cAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGACCAAGCCTTAACGAACTAAGa                                     >  2:808612/1‑143 (MQ=255)
                                                                                           aaGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGAc                                    <  1:343495/143‑1 (MQ=255)
                                                                                            aGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGAcc                                   <  1:175814/143‑1 (MQ=255)
                                                                                               tGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGAccccc                                >  1:1074828/1‑143 (MQ=255)
                                                                                                   agagAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCAc                            <  1:1021407/143‑1 (MQ=255)
                                                                                                    gAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCgtggtggt                                                                                                              >  1:1351612/1‑60 (MQ=255)
                                                                                                    gAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCgtggtggt                                                                                                              <  2:1351612/60‑1 (MQ=255)
                                                                                                    gAAAAGACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAGGGCGTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCAcc                           >  1:1220577/1‑143 (MQ=255)
                                                                                                         gACAAAATGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGACAg                      >  1:1148925/1‑143 (MQ=255)
                                                                                                             aaaTGACAGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTCAGACCCCCGCACCGAAAGGGcc                  >  1:278797/1‑143 (MQ=255)
                                                                                                               gtGACAGCCCTTCCACGAGTATGAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCgg                <  1:554310/142‑1 (MQ=255)
                                                                                                                tGACAGCCCGTCTACGAGTCTTCGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCCCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCggg               <  1:529730/143‑1 (MQ=255)
                                                                                                                    aGCCCTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGCCCTAAGACCCCCGCACCGAAAGCTCCGGGGGtt           >  2:1245802/1‑143 (MQ=255)
                                                                                                                      cccTTCTACGAGTATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCACACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCgggggagtt         >  1:885120/1‑139 (MQ=255)
                                                                                                                        tcgcTACGAGTATTAGCCTGGGTGTGATTAGCGTGGTGGTGATTATTGTCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGtttttt       <  2:176401/140‑1 (MQ=255)
                                                                                                                           cTACGAGTATTAGCCTGGTCGTGGTTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTg    <  2:815849/143‑1 (MQ=255)
                                                                                                                             aCGAGTATTAGCCTGGTCGTGATTAGCGGGGGGGTGATTATCATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGAc  <  1:219272/143‑1 (MQ=255)
                                                                                                                              acggggATTAGCCTGGTCGTGCTTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGAc  <  2:237425/138‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                         
TGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGAC  >  NC_000913/3950211‑3950479

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: