Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,757,874 T→C H225R (CAT→CGT)  b1680 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,757,8740TC99.6% 792.7 / ‑4.8 243H225R (CAT→CGT) b1680hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (143/99);  total (143/100)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

Only 100 of 247 total aligned reads displayed.

ATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTT  >  NC_000913/1757784‑1757962
                                                                                          |                                                                                        
 tAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGcc                                                                                       >  2:2846107/1‑93 (MQ=255)
 tAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGcc                                                                                       >  1:4529952/1‑93 (MQ=255)
   gAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCg                                                                                     >  1:839781/1‑93 (MQ=255)
     aCCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGga                                                                                   <  2:3111208/93‑1 (MQ=255)
     aCCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTCCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGga                                                                                   >  2:1017781/1‑93 (MQ=255)
         cccccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAac                                                                               <  1:5507350/93‑1 (MQ=255)
         cccccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAac                                                                               >  1:1844057/1‑93 (MQ=255)
         cccccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAac                                                                               <  1:2855748/93‑1 (MQ=255)
          ccccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAaca                                                                              >  1:2548393/1‑93 (MQ=255)
           cccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAacac                                                                             >  2:6448967/1‑93 (MQ=255)
           cccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAacac                                                                             >  1:4841092/1‑93 (MQ=255)
           cccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAacac                                                                             <  1:474811/93‑1 (MQ=255)
            ccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACg                                                                            >  2:861128/1‑93 (MQ=255)
              ttCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTa                                                                          >  1:2644329/1‑93 (MQ=255)
              ttCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTa                                                                          >  1:4597351/1‑93 (MQ=255)
               tCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTaa                                                                         <  2:4538718/93‑1 (MQ=255)
                cccACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTaaa                                                                        <  2:1444777/93‑1 (MQ=255)
                 ccACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTaaaa                                                                       <  1:2593711/93‑1 (MQ=255)
                 ccACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTaaaa                                                                       <  2:4442234/93‑1 (MQ=255)
                 ccACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTaaaa                                                                       <  2:756026/93‑1 (MQ=255)
                   aCGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGt                                                                     <  1:3697279/93‑1 (MQ=255)
                    cggcggCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTc                                                                    >  2:1279305/1‑93 (MQ=255)
                      gcggcATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGc                                                                  >  1:2786056/1‑93 (MQ=255)
                      gcggcATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGc                                                                  >  2:1145912/1‑93 (MQ=255)
                        ggcATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCaa                                                                <  1:1390556/93‑1 (MQ=255)
                         gcATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAAt                                                               >  2:5690683/1‑93 (MQ=255)
                         gcATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAAt                                                               >  2:6724552/1‑93 (MQ=255)
                          cATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATc                                                              >  2:1890872/1‑93 (MQ=255)
                           aTCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATcc                                                             <  2:6200364/93‑1 (MQ=255)
                            tctcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCa                                                            >  2:4589906/1‑93 (MQ=255)
                            tctcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCa                                                            >  2:5186074/1‑93 (MQ=255)
                            tctcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCa                                                            >  2:4018516/1‑93 (MQ=255)
                            tctcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCa                                                            <  2:2390465/93‑1 (MQ=255)
                            tctcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATAAAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCa                                                            <  1:3164152/93‑1 (MQ=255)
                             ctcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAg                                                           >  1:4300461/1‑93 (MQ=255)
                                cTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAgcgc                                                        >  2:3548722/1‑93 (MQ=255)
                                cTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAgcgc                                                        >  2:4140367/1‑93 (MQ=255)
                                 tGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGcc                                                       >  2:1331582/1‑93 (MQ=255)
                                    aacaaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGa                                                    <  2:822811/93‑1 (MQ=255)
                                     acaaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGaa                                                   >  2:2028578/1‑93 (MQ=255)
                                     acaaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGaa                                                   >  2:4668667/1‑93 (MQ=255)
                                     acaaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGaa                                                   <  1:6306677/93‑1 (MQ=255)
                                      caaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAAc                                                  >  1:2178314/1‑93 (MQ=255)
                                      caaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAAc                                                  <  1:5012428/93‑1 (MQ=255)
                                      aaaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGTCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAAc                                                  <  2:2534576/92‑1 (MQ=255)
                                       aaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACa                                                 <  1:4001790/93‑1 (MQ=255)
                                       aaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACa                                                 >  1:1915571/1‑93 (MQ=255)
                                       aaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACa                                                 >  1:128964/1‑93 (MQ=255)
                                        aGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACAt                                                >  2:6360422/1‑93 (MQ=255)
                                         ggCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATc                                               >  1:1524111/1‑93 (MQ=255)
                                          gCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATcc                                              >  1:1669316/1‑93 (MQ=255)
                                          gCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATcc                                              <  1:6475138/93‑1 (MQ=255)
                                           cttcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCa                                             >  2:3371093/1‑93 (MQ=255)
                                            ttcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAc                                            >  2:3284742/1‑93 (MQ=255)
                                             tcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAcc                                           >  2:1365092/1‑93 (MQ=255)
                                               tttCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCgg                                         >  2:3834804/1‑93 (MQ=255)
                                               tttCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCgg                                         >  2:4774949/1‑93 (MQ=255)
                                               tttCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCgg                                         >  1:397037/1‑93 (MQ=255)
                                               tttCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCAACCgg                                         >  2:145578/1‑93 (MQ=255)
                                                  cacaTAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGgatg                                      >  2:4312302/1‑93 (MQ=255)
                                                    caTAAAGAATGCCAATTCCGGTGGTGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGgatgat                                    <  1:2661937/93‑1 (MQ=255)
                                                    caTAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGgatgat                                    >  1:3408734/1‑93 (MQ=255)
                                                     aTAAAGAATGCCAATTCCGGTGTGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGGAATCCAGCGCCTGATCATACATCTgctgatt                                   >  1:5868056/1‑92 (MQ=255)
                                                     aTAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGatgatg                                   >  1:22647/1‑93 (MQ=255)
                                                     aTAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGatgatg                                   <  2:5321309/93‑1 (MQ=255)
                                                     aTAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGatgatg                                   >  1:318677/1‑93 (MQ=255)
                                                      tAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGc                                  >  1:150482/1‑93 (MQ=255)
                                                       aaaGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCa                                 <  1:6211129/93‑1 (MQ=255)
                                                        aaGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCAt                                >  1:3801426/1‑93 (MQ=255)
                                                         aGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATc                               >  2:2949620/1‑93 (MQ=255)
                                                         aGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATc                               >  2:2056303/1‑93 (MQ=255)
                                                         aGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATc                               >  1:6593180/1‑93 (MQ=255)
                                                          gAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCa                              <  1:248345/93‑1 (MQ=255)
                                                          gAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCa                              >  2:3705100/1‑93 (MQ=255)
                                                                 aaTTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTg                       >  2:5527853/1‑93 (MQ=255)
                                                                 aaTTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTg                       <  1:1638998/93‑1 (MQ=255)
                                                                  aTTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGa                      >  1:6024638/1‑93 (MQ=255)
                                                                  aTTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGa                      >  2:605764/1‑93 (MQ=255)
                                                                   ttCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAg                     >  1:5786556/1‑93 (MQ=255)
                                                                   ttCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAg                     <  1:5758131/93‑1 (MQ=255)
                                                                      cGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAgcgc                  <  2:3662025/93‑1 (MQ=255)
                                                                      cGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAgcgc                  <  1:1862289/93‑1 (MQ=255)
                                                                       ggTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGcc                 <  1:5054618/93‑1 (MQ=255)
                                                                        gTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCa                >  1:5473457/1‑93 (MQ=255)
                                                                        gTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCa                <  1:4592140/93‑1 (MQ=255)
                                                                         tGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCAt               <  2:4218199/93‑1 (MQ=255)
                                                                          ggggCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATc              >  1:111180/1‑93 (MQ=255)
                                                                              ccATATAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATccacc          <  1:4236851/93‑1 (MQ=255)
                                                                              ccATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATccacc          >  2:625340/1‑93 (MQ=255)
                                                                              caaTACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATccacc          <  2:1512005/91‑1 (MQ=255)
                                                                               cATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATccacca         <  2:1524909/93‑1 (MQ=255)
                                                                               cATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATccacca         <  2:1717149/93‑1 (MQ=255)
                                                                                agaCAGTTTAGGCCCGGAGAACACGTAAAACTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCTCCATCCACCAg        <  1:2036451/91‑1 (MQ=255)
                                                                                  aCAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCa      <  1:3338757/93‑1 (MQ=255)
                                                                                   cAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCAc     <  1:71444/93‑1 (MQ=255)
                                                                                     gTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACtt   <  2:5016235/93‑1 (MQ=255)
                                                                                     gTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACtt   <  2:3461602/93‑1 (MQ=255)
                                                                                     gTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACtt   <  1:293798/93‑1 (MQ=255)
                                                                                      tttACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACttt  >  2:6175609/1‑93 (MQ=255)
                                                                                      tttACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACttt  <  1:2931929/93‑1 (MQ=255)
                                                                                          |                                                                                        
ATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTT  >  NC_000913/1757784‑1757962

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