Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,002,813 C→T *127* (TAG→TAA yigG ← PRK11371 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,002,8130CT100.0% 250.8 / NA 76*127* (TAG→TAAyigGPRK11371 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/38);  total (38/38)

CTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGATATA  >  NC_000913/4002685‑4002947
                                                                                                                                |                                                                                                                                      
ctTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTtg                                                                                                                              <  2:239087/139‑1 (MQ=255)
          tGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGt                                                                                                                    <  2:141303/139‑1 (MQ=255)
          tGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGt                                                                                                                    >  1:908439/1‑139 (MQ=255)
               ttGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTg                                                                                                               <  1:403009/139‑1 (MQ=255)
               ttGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTg                                                                                                               <  2:1196167/139‑1 (MQ=255)
               ttGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTg                                                                                                               >  1:1196167/1‑139 (MQ=255)
                tGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGc                                                                                                              <  2:908439/139‑1 (MQ=255)
                 gCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCTCTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAAt                                                                                                                        <  2:869496/128‑1 (MQ=255)
                 gCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCTCTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAAt                                                                                                                        >  1:869496/1‑128 (MQ=255)
                 gCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCa                                                                                                             <  2:1515489/139‑1 (MQ=255)
                 gCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCa                                                                                                             <  1:1433922/139‑1 (MQ=255)
                            tCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGt                                                                                                                                  >  2:1225105/1‑107 (MQ=255)
                            tCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGt                                                                                                                                  <  1:1225105/107‑1 (MQ=255)
                               aGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACt                                                                                               >  2:243366/1‑139 (MQ=255)
                                    aaCCATGTTGATCACAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAg                                                                                          >  1:201727/1‑139 (MQ=255)
                                      ccATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGt                                                                                                                    >  1:477216/1‑111 (MQ=255)
                                      ccATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGt                                                                                                                    <  2:477216/111‑1 (MQ=255)
                                                cATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAAc                                                                              >  2:620619/1‑139 (MQ=255)
                                                 aTAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTg                                                                                                               >  2:1153601/1‑105 (MQ=255)
                                                 aTAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTg                                                                                                               <  1:1153601/105‑1 (MQ=255)
                                                   aaaaaTTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGc                                                                                                              <  1:855534/104‑1 (MQ=255)
                                                   aaaaaTTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGc                                                                                                              >  2:855534/1‑104 (MQ=255)
                                                    aaaaTTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCa                                                                                                         <  2:303021/108‑1 (MQ=255)
                                                    aaaaTTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCa                                                                                                         >  1:303021/1‑108 (MQ=255)
                                                        ttAAATCGGTATTTCCATTTTTCACTTAAAGAACCATCGTTCATAGATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATGCTGAAAGTGGAAAAATAACAAATCaaa                                                                      >  1:929908/1‑139 (MQ=255)
                                                                      ccATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAAcc                                                        <  1:243366/139‑1 (MQ=255)
                                                                                           aTCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGc                                                   <  2:850632/123‑1 (MQ=255)
                                                                                           aTCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGc                                                   >  1:850632/1‑123 (MQ=255)
                                                                                             cGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCtt                                 >  1:1181123/1‑139 (MQ=255)
                                                                                               ttCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTc                                >  1:1559939/1‑138 (MQ=255)
                                                                                               ttCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTc                                <  2:1559939/138‑1 (MQ=255)
                                                                                               ttCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCg                               >  2:932239/1‑139 (MQ=255)
                                                                                               ttCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCg                               >  2:1074700/1‑139 (MQ=255)
                                                                                                tCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAATTAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCgg                              <  1:932239/139‑1 (MQ=255)
                                                                                                tCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCgg                              <  1:1074700/139‑1 (MQ=255)
                                                                                                       ttCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTg                                                                                                               >  2:1371571/1‑51 (MQ=255)
                                                                                                       ttCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTg                                                                                                               <  1:1371571/51‑1 (MQ=255)
                                                                                                        tCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTa                           <  2:1496536/134‑1 (MQ=255)
                                                                                                        tCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTa                           >  1:1496536/1‑134 (MQ=255)
                                                                                                          cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtgt                    >  1:450094/1‑139 (MQ=255)
                                                                                                          cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtgt                    >  1:1246122/1‑139 (MQ=255)
                                                                                                          cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtgt                    >  1:754928/1‑139 (MQ=255)
                                                                                                           ttACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtgtg                   <  1:620619/139‑1 (MQ=255)
                                                                                                           ttACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtgtg                   >  1:1081992/1‑139 (MQ=255)
                                                                                                            tACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGtga                  >  1:941972/1‑139 (MQ=255)
                                                                                                            tACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGtga                  >  1:627747/1‑139 (MQ=255)
                                                                                                              cTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGtgatg                >  1:1124483/1‑139 (MQ=255)
                                                                                                                cATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTgatgat              <  2:1081992/139‑1 (MQ=255)
                                                                                                                  tAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATag            <  2:941972/139‑1 (MQ=255)
                                                                                                                  tAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATag            <  2:627747/139‑1 (MQ=255)
                                                                                                                   aaaTACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTg                                                >  1:1055086/1‑102 (MQ=255)
                                                                                                                   aaaTACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTg                                                <  2:1055086/102‑1 (MQ=255)
                                                                                                                   aaaTACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATaga           >  1:742255/1‑139 (MQ=255)
                                                                                                                    aaTACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag          <  2:754928/139‑1 (MQ=255)
                                                                                                                    aaTACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag          >  2:1217199/1‑139 (MQ=255)
                                                                                                                       aCTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAgatg       <  2:929908/139‑1 (MQ=255)
                                                                                                                       aCTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAgatg       <  2:742255/139‑1 (MQ=255)
                                                                                                                       aCTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAgatg       <  1:122434/139‑1 (MQ=255)
                                                                                                                         tCCATAATTATTGTTTTGATGAATCAGTAGGTTCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAatcatc                                                             >  1:957340/1‑83 (MQ=255)
                                                                                                                         tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAatcatc                                                             <  2:957340/83‑1 (MQ=255)
                                                                                                                         tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAAcc                                                        >  1:851362/1‑88 (MQ=255)
                                                                                                                         tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAAcc                                                        <  2:851362/88‑1 (MQ=255)
                                                                                                                         tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAGTCTTCGGCTATCATTGtga                  <  2:1076194/126‑1 (MQ=255)
                                                                                                                         tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCa                        <  2:309402/120‑1 (MQ=255)
                                                                                                                         tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCa                        >  1:309402/1‑120 (MQ=255)
                                                                                                                         tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag          >  2:1269572/1‑134 (MQ=255)
                                                                                                                         tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag          <  1:1269572/134‑1 (MQ=255)
                                                                                                                         tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAgatgat     >  1:1071849/1‑139 (MQ=255)
                                                                                                                         tCCATAATTATTGTTCTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACTATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGtga                  >  1:1076194/1‑126 (MQ=255)
                                                                                                                           cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGATCCATTGCCTGaa                                              <  1:131390/96‑1 (MQ=255)
                                                                                                                           cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGaa                                              >  2:131390/1‑96 (MQ=255)
                                                                                                                           cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTc                                <  2:815850/110‑1 (MQ=255)
                                                                                                                           cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTc                                >  1:815850/1‑110 (MQ=255)
                                                                                                                           cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGatat   >  2:1574097/1‑139 (MQ=255)
                                                                                                                            aTAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGatata  <  2:1246122/139‑1 (MQ=255)
                                                                                                                             tAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag          <  1:335391/130‑1 (MQ=255)
                                                                                                                                |                                                                                                                                      
CTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGATATA  >  NC_000913/4002685‑4002947

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: