Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,082,060 A→G pseudogene (601/726 nt) efeU → ferrous iron permease (pseudogene)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,082,0600AG97.5% 282.8 / ‑0.9 80pseudogene (601/726 nt)efeUferrous iron permease (pseudogene)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (37/41);  minor base . (1/1);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

GTCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCG  >  NC_000913/1081925‑1082180
                                                                                                                                       |                                                                                                                        
gtcgCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGttt                                                                                                                       <  2:778536/139‑1 (MQ=255)
     cgcAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGc                                                                                                                  >  2:1534111/1‑139 (MQ=255)
       cAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTa                                                                                                                <  2:1463826/139‑1 (MQ=255)
          ggCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCa                                                                                                             >  1:823326/1‑139 (MQ=255)
             tGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGa                                                                                                          <  2:826159/139‑1 (MQ=255)
              ggCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGaa                                                                                                         <  2:1497236/139‑1 (MQ=255)
                   cTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGcc                                                                                                    <  2:568526/139‑1 (MQ=255)
                      gTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAg                                                                                                 >  2:49495/1‑139 (MQ=255)
                        gCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCg                                                                                               <  1:1661000/139‑1 (MQ=255)
                          cATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTc                                                                                             <  2:721767/139‑1 (MQ=255)
                            ttCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTa                                                                                                                >  1:403351/1‑118 (MQ=255)
                            ttCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTa                                                                                                                <  2:403351/118‑1 (MQ=255)
                                         tGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTc                                                                              <  2:782382/139‑1 (MQ=255)
                                           aaGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTg                                                                            <  1:582895/139‑1 (MQ=255)
                                             gCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGt                                                                          <  1:625164/139‑1 (MQ=255)
                                               cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTTATGGAAGGGGttttt                                                                                                                     <  2:1155766/94‑1 (MQ=255)
                                               cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGttttt                                                                                                                     >  1:1155766/1‑94 (MQ=255)
                                               cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGttt                                                                        >  1:998064/1‑139 (MQ=255)
                                               cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGt                                                                          <  1:130767/137‑1 (MQ=255)
                                               cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGt                                                                          >  2:130767/1‑137 (MQ=255)
                                                 gATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGtttat                                                                      >  2:1324284/1‑139 (MQ=255)
                                                  aTTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTGAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCg                                                                                <  1:1644218/128‑1 (MQ=255)
                                                  aTTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCg                                                                                >  2:1644218/1‑128 (MQ=255)
                                                           ccACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATc                                                            >  1:400750/1‑139 (MQ=255)
                                                              cTTTCAGGGAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAgcg                                                                                                      >  1:1027066/1‑94 (MQ=255)
                                                              cTTTCAGGGAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAgcg                                                                                                      <  2:1027066/94‑1 (MQ=255)
                                                                ttCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGgc                                                       >  1:731298/1‑139 (MQ=255)
                                                                        aTCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGc                                                                                                                  >  2:1058770/1‑72 (MQ=255)
                                                                        aTCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATccc                                                          <  2:1217840/128‑1 (MQ=255)
                                                                        aTCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATccc                                                          >  1:1217840/1‑128 (MQ=255)
                                                                        aTCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGCTTTTGGc                                                                                                                  <  1:1058770/72‑1 (MQ=255)
                                                                         tCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGc                                              >  1:800499/1‑139 (MQ=255)
                                                                            ccTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCAtt                                           <  2:823326/139‑1 (MQ=255)
                                                                              ttCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTg                                         <  2:800499/139‑1 (MQ=255)
                                                                                         tGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGccacca                              >  2:653408/1‑139 (MQ=255)
                                                                                          gCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATcc                                                           <  1:59014/109‑1 (MQ=255)
                                                                                          gCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATcc                                                           >  2:59014/1‑109 (MQ=255)
                                                                                           cGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCa                                             <  1:748672/122‑1 (MQ=255)
                                                                                           cGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCa                                             >  2:748672/1‑122 (MQ=255)
                                                                                           cGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCAcg                            >  1:1651785/1‑139 (MQ=255)
                                                                                            ggTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                               <  2:881666/87‑1 (MQ=255)
                                                                                            ggTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                               >  1:881666/1‑87 (MQ=255)
                                                                                             gTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTa                                                                                                                <  1:63635/53‑1 (MQ=255)
                                                                                             gTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTa                                                                                                                >  2:63635/1‑53 (MQ=255)
                                                                                             gTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTg                                    >  2:9108/1‑129 (MQ=255)
                                                                                             gTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTg                                    <  1:9108/129‑1 (MQ=255)
                                                                                             gTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCAcgcg                          >  1:215837/1‑139 (MQ=255)
                                                                                             gTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCAcgcg                          >  1:1343666/1‑139 (MQ=255)
                                                                                              tGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGtttat                                                                      <  1:595999/94‑1 (MQ=255)
                                                                                              tGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGtttat                                                                      >  2:595999/1‑94 (MQ=255)
                                                                                              tGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACgcgc                         <  1:653408/139‑1 (MQ=255)
                                                                                              tGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACgcgc                         <  1:49495/139‑1 (MQ=255)
                                                                                                ctctCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                               >  1:1099689/1‑83 (MQ=255)
                                                                                                ctctCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGt                                                                               <  2:1099689/83‑1 (MQ=255)
                                                                                                 tctcAACTCACTCGCTGTTTGGCACGCTGATGGGAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCg                                                                                >  2:228442/1‑81 (MQ=255)
                                                                                                 tctcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTgg                                                                                                                   >  1:1528249/1‑46 (MQ=255)
                                                                                                 tctcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTgg                                                                                                                   <  2:1528249/46‑1 (MQ=255)
                                                                                                 tctcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCg                                                                                <  1:228442/81‑1 (MQ=255)
                                                                                                   tcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAgggg                    <  1:1534111/139‑1 (MQ=255)
                                                                                                   tcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTGGCTTTCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGc                   <  2:317527/139‑1 (MQ=255)
                                                                                                   tcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGGTTTTGGCTTTCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGc                   >  1:317527/1‑139 (MQ=255)
                                                                                                          actcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAgggg                    >  1:1280183/1‑132 (MQ=255)
                                                                                                          actcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAgggg                    <  2:1280183/132‑1 (MQ=255)
                                                                                                           ctcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGt                                                                          <  1:1597265/77‑1 (MQ=255)
                                                                                                           ctcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGt                                                                          >  2:1597265/1‑77 (MQ=255)
                                                                                                            tcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTgg                                                                                                                   >  2:771919/1‑35 (MQ=255)
                                                                                                            tcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTgg                                                                                                                   <  1:771919/35‑1 (MQ=255)
                                                                                                            tcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCg                                                                                         >  2:1623661/1‑61 (MQ=255)
                                                                                                            tcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCg                                                                                         <  1:1623661/61‑1 (MQ=255)
                                                                                                            tcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGt           >  1:27230/1‑139 (MQ=255)
                                                                                                            tcGCTGTTTGGCACGCTGATGGAAGGGGTTTTTAGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCNGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGt           >  1:1410873/1‑139 (MQ=255)
                                                                                                             cGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCtggt                                                 <  2:241074/100‑1 (MQ=255)
                                                                                                             cGCTGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCtggt                                                 >  1:241074/1‑100 (MQ=255)
                                                                                                                tGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCg       <  2:731298/139‑1 (MQ=255)
                                                                                                                tGTTTGGCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCg       <  2:1651785/139‑1 (MQ=255)
                                                                                                                     ggCACGCTGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCcg  <  2:27230/139‑1 (MQ=255)
                                                                                                                            tGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGtc          >  1:847515/1‑124 (MQ=255)
                                                                                                                            tGATGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGtc          <  2:847515/124‑1 (MQ=255)
                                                                                                                              aTGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGttta                                                                       <  1:1319274/61‑1 (MQ=255)
                                                                                                                              aTGGAAGGGGTTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGttta                                                                       >  2:1319274/1‑61 (MQ=255)
                                                                                                                                       |                                                                                                                        
GTCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCG  >  NC_000913/1081925‑1082180

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