Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,792,826 C→T *286* (TAG→TAA yibB ← YibB family protein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,792,8260CT100.0% 247.2 / NA 75*286* (TAG→TAAyibBYibB family protein, function unknown
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (42/33);  total (42/33)

TACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATA  >  NC_000913/3792699‑3792960
                                                                                                                               |                                                                                                                                      
tACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGa                                                                                                                             <  1:1176966/139‑1 (MQ=255)
tACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGa                                                                                                                             <  2:518192/139‑1 (MQ=255)
     ttttGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGa                                                                                                                        >  1:503194/1‑139 (MQ=255)
               gCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAg                                                                                                              >  2:1369685/1‑139 (MQ=255)
                          attgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttt                                                                                                   >  1:281686/1‑139 (MQ=255)
                             gattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg                                                                                                <  1:1369685/139‑1 (MQ=255)
                                   aCTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCAGCATTTTTGCCCCTa                                                                                          <  1:1729935/139‑1 (MQ=255)
                                        aTTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtg                                                                                     >  1:189745/1‑139 (MQ=255)
                                         ttGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtgt                                                                                    >  2:457082/1‑139 (MQ=255)
                                                 aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGccc                                                                                             <  2:1407883/122‑1 (MQ=255)
                                                 aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGccc                                                                                             >  1:1407883/1‑122 (MQ=255)
                                                      aaaaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAagca                                                                       >  1:1521303/1‑139 (MQ=255)
                                                      aaaaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAagca                                                                       >  2:607342/1‑139 (MQ=255)
                                                      aaaaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGACCCTAAAGTGTTACTCCTGAagca                                                                       >  2:2050146/1‑139 (MQ=255)
                                                        aaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtgt                                                                                    <  2:890212/124‑1 (MQ=255)
                                                        aaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtgt                                                                                    >  1:890212/1‑124 (MQ=255)
                                                                gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACaaa                                                             <  1:2050146/139‑1 (MQ=255)
                                                                        atcatTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCa                                                     >  2:990334/1‑139 (MQ=255)
                                                                             ttttCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTc                                                <  1:457082/139‑1 (MQ=255)
                                                                               ttCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCct                                              <  1:1280419/139‑1 (MQ=255)
                                                                                tCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             <  1:607342/139‑1 (MQ=255)
                                                                                   tATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTcc                                          <  1:1333444/139‑1 (MQ=255)
                                                                                   tATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTcc                                          <  1:990334/139‑1 (MQ=255)
                                                                                    aTCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTccc                                         >  1:104424/1‑139 (MQ=255)
                                                                                     tctcAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGa                                                                                                                             >  2:518937/1‑54 (MQ=255)
                                                                                     tctcAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGa                                                                                                                             <  1:518937/54‑1 (MQ=255)
                                                                                      ctcAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCa                                       <  2:1521303/139‑1 (MQ=255)
                                                                                         aGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCatat                                    >  2:430587/1‑136 (MQ=255)
                                                                                           ttGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt                                                                                                 <  2:652881/76‑1 (MQ=255)
                                                                                           ttGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt                                                                                                 >  1:652881/1‑76 (MQ=255)
                                                                                           ttGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAg                                  >  2:1448091/1‑139 (MQ=255)
                                                                                                aTCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             >  2:1542958/1‑123 (MQ=255)
                                                                                                aTCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             <  1:1542958/123‑1 (MQ=255)
                                                                                                aTCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAg                             >  1:499749/1‑139 (MQ=255)
                                                                                                 tCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGa                                                                           <  2:390676/92‑1 (MQ=255)
                                                                                                 tCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGa                                                                           >  1:390676/1‑92 (MQ=255)
                                                                                                     aaCGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGa                        >  2:1462577/1‑139 (MQ=255)
                                                                                                        gCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAaga                     >  2:929902/1‑139 (MQ=255)
                                                                                                          tACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGccc                                                                                             <  2:1954515/65‑1 (MQ=255)
                                                                                                          tACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGccc                                                                                             >  1:1954515/1‑65 (MQ=255)
                                                                                                               aCTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCggg              >  2:338393/1‑139 (MQ=255)
                                                                                                                 tAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGc            <  2:189745/139‑1 (MQ=255)
                                                                                                                 tAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGc            >  1:1925166/1‑139 (MQ=255)
                                                                                                                 tAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGc            <  1:1462577/139‑1 (MQ=255)
                                                                                                                   aTGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGt          <  2:499749/139‑1 (MQ=255)
                                                                                                                    tGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGtt         >  2:1078646/1‑139 (MQ=255)
                                                                                                                    tGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGtt         >  1:312617/1‑139 (MQ=255)
                                                                                                                    tGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGtt         <  2:312617/139‑1 (MQ=255)
                                                                                                                      gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCg                                                          >  2:1670317/1‑88 (MQ=255)
                                                                                                                      gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCg                                                          <  1:1670317/88‑1 (MQ=255)
                                                                                                                      gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTcc                                          >  1:1001859/1‑104 (MQ=255)
                                                                                                                      gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTcc                                          <  2:1001859/104‑1 (MQ=255)
                                                                                                                      gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAg                                      <  2:723191/108‑1 (MQ=255)
                                                                                                                      gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAg                                      >  1:723191/1‑108 (MQ=255)
                                                                                                                      gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTCCTCTCCCCAGATAGt                                 <  1:1818531/112‑1 (MQ=255)
                                                                                                                      gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTCCTCTCCCCAGATAGt                                 >  2:1818531/1‑112 (MQ=255)
                                                                                                                       aGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTaa                                                                                         <  1:819753/56‑1 (MQ=255)
                                                                                                                       aGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTaa                                                                                         >  2:819753/1‑56 (MQ=255)
                                                                                                                        gggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg                                                                                                <  1:434795/48‑1 (MQ=255)
                                                                                                                        gggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg                                                                                                >  2:434795/1‑48 (MQ=255)
                                                                                                                        gggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTa     >  1:2051617/1‑139 (MQ=255)
                                                                                                                         ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCg                                                          <  1:1042299/85‑1 (MQ=255)
                                                                                                                         ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCg                                                          >  2:1042299/1‑85 (MQ=255)
                                                                                                                         ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTaa    >  1:891205/1‑139 (MQ=255)
                                                                                                                         ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTaa    >  2:886263/1‑139 (MQ=255)
                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTcc                                               <  2:1636543/94‑1 (MQ=255)
                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTcc                                               >  1:1636543/1‑94 (MQ=255)
                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             >  2:1276302/1‑96 (MQ=255)
                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             >  2:71364/1‑96 (MQ=255)
                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             <  1:71364/96‑1 (MQ=255)
                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             <  1:1276302/96‑1 (MQ=255)
                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTcc                                          >  1:733657/1‑99 (MQ=255)
                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTcc                                          <  2:733657/99‑1 (MQ=255)
                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGTCGTTTGTAATa  >  2:1342302/1‑139 (MQ=255)
                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATa  >  1:2037014/1‑139 (MQ=255)
                                                                                                                               |                                                                                                                                      
TACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATA  >  NC_000913/3792699‑3792960

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